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ABSTRACT 


Tracks produced by 5-7 Bev protons in nuclear emulsions have been followed 
till they interact with emulsion nuclei, and a mean free path of 35-6+2-4 em 
has been obtained. About 24% of the shower tracks are found to be protons, 
all of which emerge at angles less than 20° with the direction of the primary. 
The mean multiplicity of pion production is found to be 2-6+0-3, the 
coefficient of inelasticity 0-49+0-07, and the transverse momentum of pions 
220+ 30 mMev/e. 

The angular distribution of the shower tracks has been analysed on the 
assumption that the pions are produced in a single collision system with a 
forward—backward symmetry in the centre of mass reference frame. The 
mass of the particle with which the incident proton forms such a collision, 
is shown to have an upper limit of 1:08+0-15 proton mass units. Hence 
the meson production on an average may be considered as arising in a nucleon— 
nucleon collision even with a complex target nucleus. The angular distri- 
bution of the pions in the C-system is found to be isotropic, within the 
limits of experimental errors. 


§ 1. INTRODUCTION 


Hiaeu energy protons and «-particles from cosmic radiation have provided 
extensive data on multiple meson production. Such observations have 
been made both with cloud chambers and photographic emulsions. In 
those cases where few or no heavily ionizing prongs are associated with the 
event, it has been customary to assume that the secondary particles are 
created in a single collision of the incident proton with an individual 
nucleon of the target nucleus. The energy of the incident proton is then 
estimated from the median angle of the shower tracks, i.e. the angle which | 
contains half the number of shower tracks (Dilworth et al. 1953). 

During recent years, beams of high energy protons have become available, 
and diffusion chambers and bubble chambers filled with liquid hydrogen 
have been used to investigate proton—proton interactions at various energies. 
These energies, though small as compared with those available in the cosmic 
radiation, are much above the threshold for multiple meson production. 
While photographic emulsion is not the ideal medium for studying such 
nucleon-nucleon collisions, it has nevertheless often been used for this 
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purpose (Duke et al. 1955, Cester et al. 1956, Kalbach et al. 1959). Moreover, 
proton—nucleus collisions in the emulsion are capable of providing a wealth 
of useful information. It is intended to give here an analysis of such 
collisions, initiated by 5-7 Bev protons, in relation to various aspects of 
meson production, such as multiplicity, angular distribution, inelasticity, 
etc. Although heavy mesons, hyperons and anti-nucleons, can be pro- 
duced at this energy, the cross section for their production is very much 
smaller than for 7 mesons; hence it is assumed that all the newly created 
particles are z mesons. The primary energy being known, the angular 
distribution has been used to examine the hypothesis of nucleon—nucleon 
collision. 
Fig. 1 
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§ 2. EXPERIMENTAL 


A small stack of Ilford G5 emulsions (10 em x 4em x 400) was exposed 
toa single pulse of the 5-7 Bev proton beam from the Berkeley bevatron. 
The tracks were parallel to the shorter edge and had negligible dip. 
The beam had a high flux density of ~ 5 x 10° particles per em2. Only two 
plates of the stack were used in this experiment. 
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Adopting a procedure of scanning along the track, 266 interactions have 
been obtained. An interaction has been defined as an event where the 
beam track seems to terminate, with or without the association of secondary 
track(s) of any degree of ionization. Events that look like a mere deflec- 
tion of the incident track were ignored, as they would require a detailed 
examination, which is often inconclusive, to decide between meson pro- 
duction and elastic scatter. Hach interaction was classified according to 
the number of NV» of the heavily ionizing prongs (black and grey) and the 
number, ns, of the shower tracks, i.e. tracks with a blob density less than 
1-4 times that of the beam proton tracks. The projected and dip angles 
of each shower track were measured and the space angle calculated. 


Fig. 2 
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§ 3. PRELIMINARY RESULTS 


3.1. ns and Ny Distributions 
The ng distribution is shown in fig. 1. Our mean value of Ms is 2-3, as: 
compared with 2-1, of Schein et al. (1956). Events in which ng is 1, 2 or 3 
are equally abundant and these together constitute 60% of the total. 


2N2 
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The Ny distribution is shown in fig. 2. There is a high proportion 
(viz. 23%) of events in which Ny»<1. Half the number of events are 
included in the interval Ny < 4. 

The observations do not show any correlation between Ny and ns. The 
mean multiplicity ns does not seem to vary with Np, as may be seen from the 
table below. 

Table 1. Variation of ng with Ny 


Nn OSL 2,3 4-6 7-10 SSlal 
Mean ng 2-4. 2°2 2°5 2°3 2-2 


3.2. Mean Free Path 


A total track length of 75-9 m was followed and 213 interactions obtained, 
thus giving a mean free path of 35-6+2-4cm. (53 additional interactions 


Fig. 3 
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Differential angular distribution of all the shower tracks of the composite star. 
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were found, but the path length for these was not recorded. All the 266 
interactions have been used for the further analysis.) The mean free path 
isin good agreement with other observations at the same energy. Cavanaugh 
et al. (1955), obtain 37-6+5-3cem from 50 events and the Bombay group 
(Appa Rao et al. 1956) find 34-7 + 3-4 em from 103 events. 


§ 4. Mepian ANGLE OF THE 7 MESONS 


Since the average number of shower tracks per star is small, all the events 
have been grouped together to form a ‘ composite star ’, to which it is possible 
to apply statistical considerations. The angular distribution of this 
composite star is shown in fig. 3. 

Measurements on the blob density, b, and multiple Coulomb scattering, @, 
were made on all secondary tracks longer than 5mm with a view to identi- 
fying the particles producing them. The blob density was normalized to 
that of nearby proton tracks at the appropriate depth. 1000 blobs were 
counted on the shower tracks and 3000 on the proton tracks. 


Fig. 4 
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The plateau value of the blob density has been obtained from tracks of 


electron pairs. 


It is necessary to ascertain that the scattering measurements even after 
correction for the stage noise and reading errors, are not seriously affected 
by the existence of distortion, which commonly occurs on account of the 
differential shear of the emulsion during processing. The effect of such | 
distortion is to impart a curvature to straight tracks. Observations were 
made on a number of steep grey tracks to determine the magnitude of this 
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effect in several regions. It was then concluded that the curvature due to 
this cause, imposed on tracks longer than 5mm, would introduce an 
error in the scattering measurement much less than that due to statistical 
fluctuations arising from the use of a finite number of independent cells. 

The second effect to be considered is the spurious scattering (Biswas e¢ al. 
1955) which becomes an important source of noise in measurements using 
large cell sizes. Spurious scattering is believed to be caused by local 
dislocations in the emulsion, completely unrelated to the differential shear 
mentioned above. From scattering measurements on a large number of 
beam proton tracks it was found that the spurious scattering in these plates 
was fairly small, and thus the scattering measurements may be considered 
to give a satisfactory estimate of the energy of the particle, subject only to 
the statistical errors associated with the measurements. 

The results of the measurements of blob density and pf from scattering on 
64 tracks which satisfied the length criterion are shown in fig. 4. The 
plateau value of the blob density was determined from the tracks of electron 
pairs materialising from the y-rays created in 7°-decay. The plateau value 
of ionization and the shape of the ionization curve are believed to be 
sensitive to the development of the emulsion. The curves shown in the 
diagram are similar to those of Shapiro and Stiller (1952), where the rise in 
the blob density is more gradual than that obtained by Daniel ef al. (1952). 
However, the uncertainty in the shape of the curves, and the statistical 
errors associated with the measurements of blob density and scattering, do 
not allow an unambiguous identification of all the tracks into those due to 
protons or pions. Cester ef al. (1956) have arrived at a similar conclusion, 
while investigating the interactions of 3 Bev protons in emulsion. In spite 
of these difficulties, it is possible to obtain upper and lower limits to the 
median angle of the 7 mesons. 

(a2) In the sample of measured tracks, all those with pf < 1200 Mev/c 
can be classified with confidence as 7 mesons. If the remaining, (i.e. 
tracks with p68 > 1200 mev/c) are all assumed to be protons, one obtains an 
upper limit to the number of protons in the sample. The angular distri- 
bution of the measured tracks is given below. 


Tabie 2. Angular Distribution of the Measured Shower Tracks 


—— eee EEE eee 


Angle 0-10° | 10°-20° | 20°-30°| >30° Total 
| Number of tracks with é : 
pB > 1200 Mev /c 2s - ° si ve 
Total number of tracks 43 6 4 ll 64 


Assuming this upper limit to the fraction of protons also holds in the 
corresponding angular intervals of all the shower tracks (of the composite 
star), one obtains an upper limit to the number of protons as 157+. The 


} Total number of shower particles is 606. 
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median angle of the remaining, considered to be 7 mesons, would then be 
28-9°+ 1-6°, since none of the tracks that could be ascribed to protons, 
are observed at angles greater than 20°. This median angle is an upper 
limit to the true median angle of the pions alone, since under the assumption, 
the number of pions at small angles is underestimated. 

(6) If all the shower tracks of the composite star are considered to be 
7 mesons, the median angle would be 20-7° + 0-8°. This represents a lower 
limit to the true median angle of the pions since the number of pions at 
small angles is over estimated in this case. 

An independent estimate of the median angle may be obtained in the 
following way: 

(c) From statistical considerations, the incident proton, will re-emerge 
as a proton in about half the number of collisions, and as a neutron in the 
other half. If in all cases, the secondary protons are relativistic, approxi- 
mately 133 shower tracks emitted at small angles could be ascribed to such 
protons. The median angle of the remaining, all assumed to be pions, 
would be 27-7° + 1-6°. 

(d) Notwithstanding the remarks made above, about the use of (b*, pf) 
curves, it would be interesting to see what proportions of protons among 
the shower tracks is predicted if the tracks are identified (as protons or 
pions) with their aid in the best possible way. When this is done, the 
following angular distribution is obtained. 


Table 3 
Interval 0-10° | 10°-20° | 20°-30° | >30° Total 
Number of protons 21 2 0 0 23 
Number of pions 22 4 4 1] 41 
Total 43 6 4 ll 64 


Assuming the same proportions of protons and pions to hold for all the 
shower tracks in the corresponding angular intervals, the total number of 
protons among the shower tracks is estimated (all of these are within 20°) 
to be 121; the median angle of the pions would then be 26-9° + 1:6°. 

The median angle in the neighbourhood of this value does not seem to be 
very sensitive to the errors in the estimate of the number of protons. In 
fact the results under assumptions (a), (c) and (d) are equal within limits of 
experimental errors. 


§ 5. Kinematics or Muson PRopuction 


Let it be assumed that the incident proton, colliding with a particle or 
group of particles within the target nucleus (hereafter designated as the 
collision partner) produces all the mesons in a single encounter. We can, 
then define the C-system, i.e. the coordinate system in which the centre 
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of mass of the colliding particles is at rest. This system has a velocity pe 
(in units of velocity of light) with respect to the laboratory system in which 
the target nucleus is at rest. , is given by the relation, 
B= ee 
E+, 
where p and # are the momentum and total energy of the incident proton 
and H, the rest energy of the collision partner. 
A meson emitted with a velocity B* in the C-system at an angle ¢ to the 
shower axis, i.e. the direction of motion of the incident proton, will appear 
in the laboratory system to be ejected at an angle 6. The relation between 


6 and ¢ is sind 2) 
Ve tan §= B,/B* +. cos 
where y= 1//(1—B,”)- 2 Rhee ee 


Assuming forward—backward symmetry in the C-system, half the mesons 
would be emitted within an angle 7 in the laboratory system (which 
corresponds to 7/2 in the C-system), given by 


7,6, tan n =f*. erie ree te 


The total energy of the meson in the C-system, y*, expressed in terms 
of the energy (y) and momentum (fy) in the L-system provides another 
useful relation, viz: 


yr=yy(l—PBoos8) ww se 4) 


where yf,,= yf cos @ is the component along the shower axis, of the momen- 
tum of the meson in the L-system. 

The measurements of the angles of emission, and multiple Coulomb 
scattering of the identified pions (see table 3) provide the values of y, 8 and 
B,. It will be recalled that the tracks available for the measurements are 
biased in favour of those emitted at small angles, because of the selection 
criterion. The geometrical loss for large angles has to be corrected for by 
means of an appropriate weight factor, and the true (weighted) mean 
values of y, 8 and 8, determined. 

It has been found that these values are very little affected by the errors 
of identification of the shower tracks. The ambiguity of identification 
occurs only for tracks of high energy and these are mostly emitted at small 
angles. ‘The more serious errors in the mean values of 8 and 8, arise 
mostly due to the small numbers of tracks at wide angles, and the corre- 
sponding large weight factors, apart from the large statistical errors in 
the measurement of Coulomb scattering. 

Equations (8) and (4) may now be solved for the two unknowns, f, and 
B* with the appropriate value of the median angle y. For this purpose 
it was found convenient to start with an assumed value of 6, (chosen. 
appropriate to a proton—nucleon collision) and work by the method of 
successive approximations. The mass of the collision partner is then 
obtained by substituting in (1) the final convergent value of £.. 
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From the lower and upper limits of the median angle (viz. 20-7° and 
28-9°) the limits for the mass of the collision partner M/ in units of proton 
mass are found to be: 

M yin = 0°46 + 0-05, 
M nax = 1:08 + 0-15. 

The median angles obtained under assumptions (c) and (d) give for M/ 
the values of 0-98+ 0-15 and 0-91+0-15 proton mass units respectively. 
Hence one may conclude that the meson production, on an average, occurs 
as if the incident proton suffers a collision with a particle of protonic mass, 
within the complex target nucleus. 


§ 6. IsorRopy oF EMISSION 
Conversion from the measured angles to the C-system angles may now 
be done by means of the transformation eqn. (2) of the previous section, 
with £,, y, chosen appropriate to a proton—nucleon collision. $* may be 
obtained in several ways, e.g. using eqn. (3), and the median angle under 
assumption (c) which, as we have seen earlier, gives a value of M very 
nearly equal to unity. 


Fig. 5 
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C-system angular distribution. The shaded portion represents the contribution 
due to protons. The number of these is estimated to be 145 + 20. 


This allows the C-system angular distribution to be determined for all 
the shower tracks of the composite star. This is shown in fig. 5. When 
the peak in the forward direction is interpreted as due to fast protons, the 


546 V. Y. Rajopadhye on 


isotropic character of the distribution for the pions becomes immediately 
obvious. The number of protons in the forward peak is then estimated to 
145 + 20, all of which are confined to the interval 12 cos ¢ > 0-4, which 
corresponds to 0-19° in the laboratory system. Both these conclusions 
are consistent with the observations made earlier. 

The values of 8* obtained in various ways do not seem to differ from one 
another by more than a few per cent. To allow for a possible small uncertainty 
in the value of 6*, the C-system angular distribution has also been deter- 
mined for two more values of £,/8*, viz. 1 and 0-92, and in all cases it 
is found to be closely similar. Thus it may be concluded that isotropy 
obtained is a genuine effect. 


Fig. 6 


Corrected number of secondary tracks 


Transverse momentum 


P, histogram for all shower tracks, corrected for the geometric loss. The dotted 
line represents the actual number of tracks measured. 


Isotropy of emission for pions produced in 6 Bev proton interactions 
has been obtained by King (1958) and Kalbach e¢ al. (1959). However, 
Daniel e¢ al. (1959) have obtained a forward—backward peaking, an effect 
probably due to their small statistics, involving the use of only those 
tracks, which were identified as mesons by (b*, p8) curves, and the large 
weight factors to correct the geometrical loss. 


§ 7. TRANSVERSE MoMENTUM 
The importance of this dynamical quantity for the understanding of 
high energy nuclear interactions was first realized by Nishimura (Koba 
1956). Extensive measurements on the transverse momentum P » on the 
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charged and neutral pions emitted from jets were made by Edwards et al. 
(1958). They find P, to be 500 Mev/c, over a very wide range of primary 
energies. Recent observations by the Japanese Emulsion Group (1959) 
on the P, of neutral pions in jets with primary energies 10!2+ 10!ev show 
a narrow distribution with a peak at about 400 mev/c. They also find no 
significant variation with either primary energy or pion energy. 

In the present experiment, the transverse momentum of all the secon- 
dary particles has been determined from the measurements of the Coulomb 
scattering and the angle of emission. It does not seem to depend upon 
the identity (7 or p) assigned to the particle when pf, > 1200 mev/c and thus 
it has the same value for all secondary tracks. The histogram of P,, for all 
shower tracks, in which the observations are duly weighted to compensate 
for the geometrical loss is shown in fig. 6 from which the mean value of P, 
is found to be 220 + 30 Mev/c, somewhat smaller than that observed in ‘jets’. 
It should be noted that a large number of tracks available for measure- 
ments are emitted at small angles. Constancy of P, would require very 
high energies for the mesons emitted at small angles. This is clearly 
impossible in view of the finite energy of the incident proton. With the 
small statistics available, it isnot possible to draw any conclusions regarding 
the variation, if any, of the transverse momentum with the angles of 
emission. 


§ 8. ENERGY BALANCE 


It is possible to see how the kinetic energy of the incident proton is 
accounted for in the interaction. It is convenient to consider all the inter- 
actions together and utilize the results of the scattering measurements on 
the sample of 64 tracks. 

As a starting point, the average kinetic energy of the protons and the 
total energy of the pions, in the different angular intervals may be deter- 
mined under assumption (a), §4. The relative numbers of protons and 
pions shown in table 2 then allows the energy of all the shower tracks of the 
composite star to be calculated. A correction has to be made for the 
energy carried away by the neutrons, and 7° mesons emitted. About 
66°% of the collisions would be with the heavy nuclei Ag and Br, in which - 
the neutron proton ratio is 13:10 and thus collisions with neutrons would 
be more numerous in that ratio. In the case of light nuclei C, N and O, 
probability of a collision with a proton or a neutron is the same. From 
all these considerations, the energy released per star, that accounts only 
for all shower particles (charged and neutral) is found to be 3-5 + 0-8 Mev. 

Thus a lower limit is set for the energy of the shower tracks, since the 
assumption (a) causes the average energy of both pion and proton groups 
to be underestimated. Assumption (d), which gives a higher proportion 
of pions, may now be used as comparison, to calculate the average energy 
of protons and pions. From this, one obtains the total energy of the 
shower tracks per star to be 4-1 + 1-1 Bev. 
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To this must be added the energy of the black and grey tracks. Froma 
study of 9 Bev proton stars in emulsion Beliakov et al. (1958) have found 
the energy (H’) of black and grey tracks in light and heavy nuclei to 
be respectively 470+70 and 1400+ 120mev. Considering the propor- 
tion of light and heavy nuclei in emulsion, the average energy per star 
is ~lBev. This is in good agreement with that deduced from cosmic 
ray data (Camerini ef al. (1951), Powell et al. (1959)). The latter give 
HE’ =124N), +30 Mev which corresponds to 860 Mev in the present experl- 
ment. Thus the total energy released per star agrees well with the primary 
energy. 


Fig. 7 
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Histogram of the C-system kinetic energy for the mesons, corrected for the 
geometric loss. The dotted line shows the actual number of tracks. 


§9. Pron Muurrenictry anpD INELASTICITY OF COLLISION 


The number of tracks considered to be due to protons among the 606 
shower tracks varies according to the assumptions made, but it seems to 
lie between fairly definite limits of 121 and 157. The number of mesons 
then would be uncertain to that extent. Mesons emitted backwards in 
the C-system would appear slow in the laboratory system, and would 
produce grey tracks. These have not been considered in the present 
analysis, as their number is presumably small. Assuming the number of 
protons to be 145 + 20 as obtained from fig. 5, the number of charged mesons 
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per star is found to be 1-7. According to the principle of charge indepen- 
dence in meson production, the ratio of neutral to charged mesons is 1: 2. 
Thus the total pion multiplicity would be 2-6+0-3. This may be com- 
pared with 2-5 as predicted by Fermi’s theory. 

The inelasticity of collision is defined as the fraction of the available 
energy in the C-system, that is radiated off as mesons. The average 
energy of the mesons may be computed in several ways. When combined 
with the mean multiplicity, this would allow the inelasticity to-be deter- 
mined. 

The average velocity of the mesons, £*, has been found to be 0-91. 
This corresponds to a kinetic energy 210 Mev. 

The eqns. (3) and (4) of §5, when solved for 8* and 8,, with the median 
angle 26-9° under assumption (d) and 27-7° under assumption (c), give B* 
as 0-935 and 0-931 respectively. These correspond to energies of 285 
and 250 Mev. 

On the basis of the identification of mesons made under assumption (d), 
their energies and momenta have been calculated from the measurements 
of multiple Coulomb scattermg. ‘Since the mesons are assumed to be 
produced in a nucleon-nucleon collision, their C-system energy can be 
calculated by means of appropriate Lorentz transformation equations. 
Figure 7 shows the C-system energy spectrum of the mesons, thus obtained, 
with due corrections for the geometrical loss. The weighted mean kinetic 
energy is found to be 280 + 80 Mev. 

The errors of identification of the tracks do not seem to affect the mean 
value obtained above. Using assumption (a), for discriminating pions 
from protons we get nearly the same spectrum, and the same mean energy. 
The small weight factors associated with the tracks at small angle do not 
alter noticeably the spectrum by their inclusion or exclusion. 

All these results may not be considered inconsistent with one another, 
since the statistical errors are rather large. Taking the mean energy as 
280 + 80 Mev, the inelasticity would be 0-52 + 0-15. 

Since the statistical errors in the mean energy derived above are too 
high, we have preferred to use the results of Kalbach e¢ al. (1959) who find 
the mean energy to be 236+20Mev, from measurements on 143 pion 
tracks. This gives the inelasticity as 0-49 + 0-07. : 


§ 10. Discussion 


The median angle method for the analysis of meson showers is based 
on the assumption of an on-the-average forward—backward symmetry in 
the C-system. With events of low multiplicity as in the present experi- 
ment, the median angle of the shower tracks would be considerably smaller 
than that for the mesons only, because the proton often emerges from the 
interaction with only a small change of direction. The primary energy 
deduced from the median angle of the shower tracks would then be over- 
estimated (Haber-Schaim 1956). At 6Bev primary energy, the mean 
velocity of the mesons f*, in the C-system is found to be somewhat higher 
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than the velocity of the C-system itself. Ignoring this would have the 
effect of under-estimating the primary energy. The error involved, 
however, is of the same order as that due to the statistical error in the 
median. angle. 

The fluctuations in the median angle with Ny observed by Haber- 
Schaim, (1956) are probably due to the small statistical weight of his data. 
In the present experiment, the variation of the median angle with Ny is 
found to be as follows: 


Table 4. Variation of the Median Angle of Shower Tracks with Vp 


Na 0 1 <4 >5 
Median angle 16-3420 | 14041-7 | 17-5410 | 25-041-4 


A difference between, the median angles for small and large Ny groups 
has been also found by Beliakov e¢ al. (1958), in a study of 9 Bev proton 
stars. For collision with protons, light and heavy nuclei they find the 
median angles for shower tracks to be respectively 18°, 25°, and 28°. 

Events with small Ny are generally considered to be good approxi- 
mations to nucleon-nucleon collisions. The number of such events is 
not very large in this experiment. However, adopting the same pro- 
cedure as before, the C-system angular distribution for pions, in these 
events has been determined. A plot similar to fig. 5 shows isotropy of 
emission for pions, allows the number of protons in the forward peak and 
hence the median angle of the pions to be determined. Although the 
median angle of the shower tracks is observed to depend upon the Vy 
group, there is no evidence that the median angle of the pions shows a 
similar behaviour. The total multiplicity of the pions for the small Ny 
group is found to be 2-04+0-36 which, although lower than Fermi’s 
prediction, is in agreement with Daniel et al. (1959), who made a more 
detailed examination of events with Ny<4. The uncertainty in the 
number of protons causes large errors, and a comparison with theory 
loses much of the significance. 

Considering all events, an upper limit to the number of fast protons 
per star is found to be 0-:59+0-06. For events of small Ny, the 
corresponding figure seems to be significantly higher, although the 
uncertainty in the number of protons estimated introduces a relatively 
greater error. 

It is to be noted that the conclusions of the experiment are based on 
the analysis of shower tracks only, and from events of all Ny grouped 
together. Slow mesons, of which the number is presumably small, are 
not considered. Since an upper limit to the proton contamination of the 
shower tracks could be set up, it was possible to deduce that the meson 
production occurs as if the incident proton suffered a collision with a 
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particle of mass no greater than 1-08 + 0-15 proton mass units. Thus the 
hypothesis of a nucleon-nucleon collision seems to be justified. More- 
over all the fast protons are found to be emitted at angles less than 20° 
and the majority of them at very small angles. 
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ABSTRACT 

An attempt is made to interpret the horizontal ‘ drift velocities ’ observed 
in ionospheric fading-wave experiments by examining statistically the 
fluctuations in space and time of the density of ions and electrons. The 
ionization is assumed to be convected by a turbulent flow in the neutral air, but 
the magnetic field and the dynamo electric field are also taken into account. 
The drift velocity to be calculated is defined by means of the space-time 
correlation function. 

For the D region and lower E region, the estimation of this drift velocity 
is found to be mathematically practicable, using certain assumptions. 
Results are obtained which show the relative importance of the mean aero- 
dynamie velocity and of the electric field in determining the drift of irregu- 
larities in electron density. 

For higher regions, the theoretical problem is much less tractable, for 
reasons discussed here; but it is pointed out that in this case it is unlikely 
that the irregular behaviour can be ascribed to turbulence anyway, and that 
the problem of the drift velocity cannot properly be investigated until a 
satisfactory explanation of the cause of the irregularities is available. 

Some remarks are offered on the relation of this theory to the observations 
and to the dynamo theory. 


§ 1. INTRODUCTION 


OBSERVATIONS of the fading of radio waves reflected from the ionosphere 
reveal that the density of ionization has random, or ‘irregular’ fluctuations 
superposed on the gradual variation with height. If two or more spaced 
receivers are used, comparison of the records reveals a noticeable impression 
of large-scale movement of the irregularities. The procedure for evaluating 
the velocity from the fading records was given by Briggs e¢ al. (1950). 
Further random movements are superposed on this general ‘drift’. 

In the E region, since the gas is lightly ionized and collisions are frequent, © 
the movements of the ionization are related to those of the neutral air. It 
is thought that in this region the air may often be in turbulent motion; if 
so its movement is described in terms of a mean ‘wind’ on which are 
superposed irregular turbulent velocities. 1t might then be expected (and 
indeed has sometimes been assumed) that the mean ‘wind’ in the neutral 
air is equal to the mean ‘drift’ of the ionospheric irregularities. This 
would be the case if the density of the plasma behaved as a simple convected 
function, like temperature, but unfortunately the Earth’s magnetic field 
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makes the ionospheric problem much more complicated than that. The 
purpose of this paper is to investigate the relation between the two mean 
velocities. It will be shown that in general they are not the same, and 
that the relationship between them depends on the height. 

Clemmow et al. (1955) and Clemmow and Johnson (1959) gave idealized 
solutions of the relevant equations which showed clearly that results of this 
kind must be expected, and this investigation began as an attempt to 
generalize their work. The present approach, however, differs from theirs 
in several respects. In their formulation, only a uniform motion of the 
air was permitted ; a model ‘irregularity’ of cylindrical or sinusoidal form 
was introduced at time t=0 (in the plasma only), and its subsequent 
behaviour in imposed electric and magnetic fields was found. The present 
paper attempts instead to discuss the problem statistically. The situation 
is supposed to be one of equilibrium: the turbulent flow continually 
creates and distorts irregularities, but these, if left to themselves, would 
gradually be smoothed out by diffusion. It is therefore necessary to make 
some assumptions about the nature of the air turbulence, and this in itself 
poses some difficult problems arising from the effect of the buoyancy forces. 
By restricting attention to the problem of the drift velocity, these assump- 
tions are kept to a minimum, and we shall only consider the simplest which 
can reasonably be made. This will at least make clear the nature of the 
problem, and it will also become apparent that to discuss further features 
of the radio fading, i.e. to discuss the random movements in addition to 
the statistically defined ‘mean drift’, would require much deeper analysis. 
This has recently been initiated by Miinch and Wheelon (1958), but so far 
their discussion of temporal correlations concerns only the components of 


turbulent velocity and not the density of a plasma convected in a magnetic 
field. 


§ 2. OUTLINE oF THE METHOD 


To carry out this programme, we have first to decide upon the assumed 
form of motion of the neutral air (assuming that the reaction of the plasma 
on it is negligible). §3 is taken up with a discussion of this, together with 
some points concerning the interpretation of the fading pattern on the 
ground. 

The standard dynamical equations for ions and electrons together with 
Maxwell’s equations then determine, in theory, the motion of plasma in 
external electric and magnetic fields. It is convenient to write all field 
quantities as a mean part together with arandom part whose mean is zero. 
The mean electric field is then the field of the ‘atmospheric dynamo’ 
(Baker and Martyn 1953). This is not genuinely external of course, but 
it is arbitrary so far as our equations are concerned since we suppose all 
mean quantities to vary so slowly (with height, with horizontal distance, 
and with time) that they may be regarded as constant in a sample of space- 
time large enough to be statistically significant. The random parts are 
supposed statistically stationary in such a sample, and the equations 
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governing them are developed in §5. As in ordinary turbulence theory, 
it is useful to work with the Fourier transforms of the random quantities, 
so we derive in §6 an ‘equation of motion’ for each Fourier component 
N(k,t) of the density fluctuation n(x,¢) of the plasma. This exhibits the 
temporal changes in N(k,t) as being due to a ‘convection’ dependent on 
the mean wind and electric field, and on the direction of k; 2 ‘decay’ 
due to diffusion, also directionally dependent; an ‘intermodulation’ with 
other Fourier components via their electrostatic fields; and a complicated 
process of ‘production’ by the turbulent flow, which includes amongst 
other things the mechanism first described by Dungey (1956). 

Before entering into these dynamical considerations, however, we give in 
§ 4 the natural definition for the ‘drift’ velocity of the pattern of a random 
scalar function in terms of its space-time correlation function. This is a 
straightforward generalization of that given by Briggs e¢ al. (1950) for an 
isotropic pattern in the amplitude of the signal on the ground. Since the 
statistical definition of the drift is thereby the same as the one used in 
analysing the records, it seems reasonable to identify the two, assuming 
the fading to be due to small structure in the ionosphere, as discussed in § 3. 

The results of § 6 are applied in §§ 7 and 8 to the drift so defined ; a tentative 
explicit expression for it for heights up to about 120 km is found assuming 
an isotropic spectrum, as predicted by Howells (1960). Other literature on 
turbulent convection (Villars and Weisskopf 1955, Gallet 1955, 1957, 
Wheelon 1957, Batchelor 1959, and, with the magnetic field included, 
Howells 1960) concentrates on trying to estimate the spectrum function 
of k(=|k|), apparently without general agreement. So it is interesting to 
note that our result is independent of this. A graphical representation 
of the way the result depends on height, through the E region, is given, 
from which one can deduce the relative importance of wind and electric 
field in determining the drift. 

For greater heights, a peculiar difficulty enters and the whole situation. 
becomes rather obscure. ‘The reason for this is discussed in §9. This also 
provides an opportunity to make some remarks on the analogous problem 
for the F region. ; 

A general discussion of the result and its relation with both the radio 
observations and the dynamo theory is provided in the concluding section. 


§ 3. Mecuanics OF THE Non-IONIZED AIR IN THE LOWER LONOSPHERE 


As we have already hinted, no adequate theory of the motion of the 
neutral air in the E region exists—indeed if it were not for the kind of 
difficulty which it is the purpose of the present paper to try to overcome, 
one would hope to use the radio observations to help to elucidate some 
features of the air motion. To make progress we have to adopt some 
idealized model, so some explanation is due on the question of how likely 
this is to be near the truth. 

Present indications (mainly from data concerning the rate of deformation. 
and diffusion of meteor trails) are that the air motion is to be pictured 
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somewhat as follows. At the smallest scales, kinetic energy is converted 
into heat by the action of viscosity and we have a ‘dissipation range’ in 
the ordinary way, with buoyancy forces quite unimportant. At rather 
larger scales, dissipation gradually becomes less important, and buoyancy 
more important, but there may be an ‘equilibrium range’ in which neither 
is very effective, the main feature being the interaction among all pairs of 
Fourier components of the velocity field, arising from the non-linear 
inertial term in the Navier-Stokes equation. So far, the energy spectrum 
function would be roughly isotropic. But at still larger scales the effect 
of buoyancy becomes the dominant feature, and the spectrum becomes 
anisotropic. The large ‘ flat’ eddies of about 6 km depth and up to 150 km 
in horizontal extent, observed by Greenhow and Neufeld (1959) come under 
this heading. Most of the kinetic energy resides in these eddies. An 
alternative description of the large-scale motions is obtained by ignoring 
the non-linear term so that one obtains the linearized theory of ‘ gravity 
waves’ inastratified medium (Hines 1959). This gives a useful insight into 
the general effect of the buoyancy forces and one can easily re-examine 
the non-linear term later, to establish the error involved, which, if small, 
can be interpreted as coupling between the various modes. (There is of 
course another class of motions, associated with compressibility, which at 
small scales are sound waves, but which at scales comparable with a scale- 
height are also profoundly modified by buoyancy. This is why Hines has 
‘ two turbulent sequences’.) Finally, on a scale still larger than the large 
eddies or gravity waves, are the world-wide meteorological motions. 

If this description is correct, the usual notion of degradation of energy 
from larger to smaller scales and finally into the dissipation range, may be 
used, but care should be exercised in estimating e¢, the rate of dissipation 
per unit mass, from the characteristic size L and velocity V of the energy 
containing eddies. or isotropic turbulence in a uniform fluid one usually 
writes «= V#/L, but recently reported observations of meteor trails indicate 
that « is much smaller than V*/L if the latter refers to the large gravity- 
dominated eddies. A qualitative explanation appears to be that gravity 
and the density stratification enable the flow to take on a more ordered 
structure whereby the excitation of smaller-scale motion through the 
non-linear terms is less than would be the case for a flow of comparable 
V and Z in a uniform fluid (assuming no special boundary conditions). 
Thus in our case, to apply the above formula, V and L should refer to the 
rather smaller-scale (and less energetic) part of the spectrum where 
buoyancy is just beginning to be dominant, say V x 1m/sec, L somewhat 
less than 1 km, as a very rough estimate. Indeed it might be better not to 
follow the terminology of the authors just cited, and reserve the word 
‘turbulence’ for the approximately isotropic part of the spectrum, so 
excluding the gravity waves. 

The reader will find more discussion of many of these points in the 
Proceedings of the International Symposium on Fluid Mechanics in the 
Ionosphere (Journal of Geophysical Research, December 1959) (and the 


Statistical Theory of Ionospheric Drifts 557 


author is indebted to many of the pegakers at this meeting for making 
clearer these ideas). 

This description has been given in some detail partly to emphasize the 
great complexity of the problem even before introducing the plasma in a 
magnetic field, and partly to show the significance of a fundamental 
assumption we shall make. This is to suppose that the radio fading, and 
the associated measurements of drift velocity, are to be explained in terms 
of the small scale, roughly isotropic, part of the spectrum. The wavelength 
for reflection from the E layer is about 100m, and the size of a Fresnel zone 
is about 5km. So any quite small structure in the distribution of plasma 
density is certainly effective in producing small-scale structure in the 
radio-wave diffraction pattern on the ground. On the other hand it is 
quite true that even with large-scale variations only, small-scale structure 
can appear in the diffraction pattern, as for example in an interference 
pattern due to specular reflection from a small number of points in a 
“wavy ionosphere. The theoretical problem of interpreting the observed 
distribution of signal on the ground in terms of its cause in the ionosphere 
is discussed by Ratcliffe (1956), who points out that the fading pattern does 
not contain enough information to distinguish between these possibilities. 
Briggs (private communication) suggests to the author that the highly 
irregular appearance of records obtained by him (by the rapid frequency - 
change method) for E region contrast rather markedly with those for the 
F region, when displayed as a panoramic record. The F region ones are 
more regular suggesting specular reflection, thus leading one to think that 
in the E region it is the small-scale fluctuations which are responsible for 
the fading. If it is subsequently shown that large structure plays an 
important réle, it will be necessary to study the convection of plasma by 
gravity waves in a magnetic field instead of turbulent convection as in this 
paper. Such a study is in any case relevant to the part of the ionosphere 
where viscosity is so large that no equilibrium range exists at all, and work 
on this problem has already been published by Hines (1955, 1956). 
Estimates of the height at which this might begin to be the case vary; 
one, due to Stewart, is at about 120 km (remarks at International 
Symposium, loc. cit.) 

Our model, then, is that the flow of neutral air consists of roughly 
isotropic inertial and dissipation ranges of turbulence, superposed on some 
uniform mean flow; the latter may incorporate, besides meteorological 
‘winds’, flow due to any gravity waves having a scale larger than a region 
scanned in making an observation by radio and uniform throughout the 
time required to make a fading record. Although this is rather an idealiza- 
tion, it enables us to concentrate on the special effect of the magnetic field 
on the convective process. 


§ 4. Tar Drier Vetocrry or 4 RANDOM ScALaR FIELD 


Suppose we have an ensemble of realizations of a random scalar n(x, ¢), 
stationary in space and time, with mean zero. Although n will later refer 
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to the ionization density, it will be convenient here to suppose the variance 
has been made unity by normalization. Our standard notation for 
Fourier transforms will be to use the corresponding capital letter, as in — 


nxt) = | NU, t)exp (ik:x)d*k. (4, 3 
The space-time correlation function is defined by the ensemble average 


ale, 7) = n(x, NOCH Gee ee ee rr 


and is independent of x and ¢ (by the stationary property). By Schwartz’s 
inequality this has a maximum (namely unity) at §=0,7=0, so assuming 
a Taylor expansion exists, we have 

p(&,7) =1—4[a,,€,€, — 2b,E7 +7?) + «.. ee 
where the quadratic form in brackets must be positive definite. (Summa- 
tion over i, 7, being assumed as usual). The coefficients a,,, b;, c are given 
by the second partial derivatives of p. 

With a suitable ergodic assumption, it is immaterial whether the average 
in (2) is taken over space, time, or the ensemble of realizations. Imagine, 
for a particular realization, an observer calculating the ‘cross correlation ’ 
between the values of n(x, ¢) at pairs of points separated by & and at two 
instants ¢and¢+7, where 7 is small; the averaging process being performed 
only over values of x. He then selects the value of § which maximizes 
the correlation, and defines §/r= W to be the drift velocity during the 
interval. An instantaneous drift is clearly provided by lim W if the 


70 


limit exists, which it does if (3) is valid; in fact one easily finds 


Weather 


ij? 


where a is the matrix [a,,;] and b the column-vector [;]. The correlation 
for small 7 given by (3) when § = Wr reduces to 


p(Wr,7)=1—3[c—b.at.bl?+... 2 . . . (5) 


so the square bracket is a measure of the random changes which remain 
after the drift has been removed. For §=W7+y, say, one has 
pl, 7)=1 — $1a,;+ (c—b.a™.b)r7]4+..:. Sh oa eel Od 
The problem of defining a mean drift was considered by Briggs ef al. 
(1950) for the case of random diffraction patterns on the ground produced 
when radio waves are reflected from the ionosphere. Their approach was 
conditioned by the fact that one has only a small number of aerials to 
observe that pattern so p can only be obtained for two or three values of 
§. Ifthe assumption as to the nature of p which they were thereby forced 
to make is true, their V, and V, are the same as would be given by (4) in 
this two dimensional case. These authors also assumed a spatial isotropy, 
80 a,;=a9;;, say, and their V,, is defined as (ac — b?)1?/a in our notation, with 
(6) becoming 


p(&, 7) =1—ha(n?+ V272) +... 5 TS SRR ECD) 
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When n(x, ¢) is given as a Fourier transform, (1), the definition (2) of p 
becomes 


6 7)= | [exp [ik . (x + )]exp (ik’. x) N(k’, DN(k, f-+ 7) d3 kd? h! 


= | [exp tik. 5) exp(tk”.x) V(k” —k, t)N(k,¢+7) dka3 kh”. 


At this point we extend the averaging to include an integration over x in 


a large volume 87%c, remembering p(§,7) is actually independent of x, 
to obtain 


p(§,7) = (Satay {| exp (7k. §)exp(tk”. x) N*(k—k”", t)N(k,¢+7) 
x Bkd? kh" dex 


where we have made use of the reality condition N(—k,t)=N*(k, 6). 
The usual rule for the inversion of Fourier transforms then gives 


plE,7)=o [exp (ike. &) VGN 1G 4 7) oer (S) 


Comparing with (3) and (4), we have 


0p - a 
Oy | Se Be a een ck 
a; car) o [ae |V(k, t) |? dk, 


x Bei, k,t) [et d® b. 
b, (55), ps faa (k, NON (k, ford? k 


(9) 


These are of course actually independent of t, since |N(k,t)|? is so; this 
also shows incidentally that b as given by (9) is real, as it should be. 

The above work uses Fourier integrals formally, without regard for 
convergence. Rigorous exposition of manipulations of this type can be 
given, provided the Stieltjes integral notation is used, see for instance 
Moyal (1952). We shall recover finite expressions on evaluating a~!.b. 


§ 5. DYNAMICS OF THE IONS AND ELECTRONS 


Detailed considerations of the effects likely to be important in the type 
of disturbance to be investigated, have been given by Dungey (1956) and 
by Clemmow and Johnson (1959) so need not be repeated here. We shall ~ 
assume : 

(i) Fluctuations in the magnetic field are to be neglected. A convenient 
choice of axes makes B= (0,0, 6). 

(ii) The plasma consists of electrons, and ions of one kind only, singly 
charged, and the dynamical processes take place so rapidly that ionization 
and recombination processes can be neglected. Inertial and gravitational 
terms are also neglected. 

(iii) The number densities of ions and electrons remain almost exactly 
equal, the difference being much smaller even than the perturbations in 
number density. (Clemmow and Johnson do not assume this at the outset, 
but the result of their linearized calculation is that additional modes of 
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oscillation exhibiting non-neutrality, but decaying very rapidly, enter the 
theory ; they are neglected here.) This is the usual procedure, (see Spitzer 
1956, p. 16) which is valid provided the length-scales of the problem are all 
much greater than the Debye distance, which for the ionosphere is a few 
millimetres. The number of variables is thus reduced by one, but Poisson’s 
equation need not be used. This is, however, contrary to Booker’s (1956) 
treatment, which introduces both ‘ionic’ and ‘electronic’ irregularities. 

We shall write script letters for the following mean quantities: 
for the density of ions and electrons, W*,Y— for their velocities, and W@ 
for the velocity of the neutral air. The random parts are denoted by the 
corresponding small letters n, v+, v-, u. But the mean electric field is 
written E, and the random part — Vd¢. Capital letters V, @, are the 
Fourier transforms as in (1). We shall frequently write symbols or 
whole equations with alternative signs +, for brevity, to indicate ions and 
electrons respectively. Thus m+, +e are the mass and charge, while v= 
are the frequencies of collisions with neutral particles (ion—electron collisions 
being neglected). We assume partial pressures. £7 for ions and electrons 
at constant temperature 7’. The gyro-frequencies are w+=eB/m+, and 
as in Baker and Martyn (1953) the useful angles « are defined by 


tan a= = w+/v+. 2 ost Se eee ee eee 

The general equations of continuity and momentum are 
= = —div(W¥ +), eerie ee 

= 26 eR ae 
mide ¥*— We) Fey =XB — 0 — VA oe cane 


which together with curlE=0 determine the motion of the plasma. 
Equation (12) may be solved for Y= in the form 


aaa He EL Mm er) ae Aas 
m+p+ m=v>A 
where the matrices M+ are defined by 
cos? a+ +sina+cosa+ 0 
M== | Fsinatcosat COs? x= ae See) 
0 0 1 


Since the script letters and E were defined as the steady mean values, 
all the derivatives in (11)-(13) are actually zero; we retained them formally 
so that the equations in the disturbed state are then obtained by replacing 
U by U+u, etc., and E by E~V¢. One then subtracts the equations 
satisfied by the mean quantities, 


eG 
wae +®) oS ee ae 
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to obtain 
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so using (11), and going into suffix notation, 


on 0 
— = —— aff” wy, ‘ 9). =E 
BE 7 3g, EAN to 
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Except in the non-magnetic case, when M* become unit matrices, it is 
impossible to eliminate explicitly the internal potential ¢ which couples 
the two simultaneous diffusion equations (17) and makes the problem 
non-linear. We are, however, only concerned to find enough about the 
solutions of (17) to make an estimate of the drift velocity W of §4, and to 
do this the equation will be converted into one satisfied by the Fourier 
transforms. 


(17) 


§ 6. THE EQuaTIon oF Motion oF A FouRIER COMPONENT 


One simple approximation in eqns. (17) will be made at once: the frac- 
tional changes in density are known observationally to be quite small, 
though the density gradients involved are larger than the steady gradient 
with height which we have already neglected. Hence we can ignore n in 
comparison with ./ in factors where n is not differentiated, i.e. we replace 
the two brackets (4 +7) by WV in (17). 

After doing this and substituting expressions of type (1) into (17), we 
obtain, using some straightforward changes of variable, two coupled 
integral equations for N(k,t) and O(k, ¢): 


Mike yore i He maT neers 
I = HIVE MN (bt) + My | FS bak (kf) — Bhp (lef 
é , im) , ! , 
| (ic, — ky’) N(k’, 1) (kk —k’, t)d3 ke 
=i [kj U (kek ON) BR IHU kD} eR ieee io): 


The terms on the right-hand side of (18) have the following interpretation. 
The first represents a convection of the Fourier components of wave 
number k in the ions or electrons with the mean drifts W* or W— 
respectively ; the first in the large bracket is the effect of the Fourier 
component with the same k in the electrostatic field, set up to prevent 
charge separation. The term involving £7’ gives the effect of molecular 
diffusion in smoothing out the Fourier component k, while the last term 
of all can be identified as Dungey’s mechanism (Dungey 1956) for 
producing irregularities in the plasma by means of a magnetic field. We 
note that in the absence of the field, M;,+ reduce to 6,;, and the term 
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becomes small because k.U(k,t)+0 since the turbulent velocities are 
much smaller than the sound velocity; the term then represents the 
compressions in the neutral air being repeated in the plasma, an effect 
rejected by Villars and Weisskopf (1955) and by Gallet (1955) as being 
too small to explain the strength of irregularity observed. (Maeda et al. 
1959, suggest similarly that the compressions are repeated in the plasma 
via the ionizing process rather than mechanically: this seems to the 
present writer to be far too small an effect, for similar reasons). 

To complete our description of the right-hand side of (18), there are the 
two integral terms representing intermodulation between all pairs of 
Fourier components k and k’. In the first case it is through the internal 
electric field, while the second integral is a generalization of that arising in 
the non-magnetic case, whereby the turbulent velocity field acting on 
the local gradients in n(x,t) due to the existing irregularities, continually 
redistributes irregularity throughout the spectrum Vk, f). 

The gentle gradient in the mean ./ is also involved in this last mentioned 
process, but (as emphasized by Batchelor (1959, p. 118) but not appreciated 
by some of the earlier writers) the gradients due to the irregularities them- 
selves are much more important, except at the smallest wave numbers. 
So (18) is still valid in the dissipation and equilibrium ranges, but omits a 
term describing the ‘input’ at small wave number, which however is 
immaterial to us in what follows. 

It will be seen from (18) that the difficulty of eliminating 4 between the 
ion- and electron-equations is at least as great as before; we have removed 
spatial differentiations at the price of introducing integral equations. 
Nevertheless (18) supplies the equation which is sufficient for the present 
purpose, namely that obtained by eliminating the term in which ®(k, f) 
appears explicitly, without regard for the term in which ®(k — k’,#) appears 
under integration. If 6 is the angle between k and B, so that 


cos 0=k,/k MRC bei ho ef), 
and introducing 
St = M,j+k,k,/k? = cos a+sin?@+cos?6 . . . . (20) 
the appropriate linear combination of equations (18) is 
ON (k,t 
. ) = i W(6). KN (k,t) — Dy(6)k2N (ke, 


+ Pis(8) | kj (kj — kj’) N(k’, t)O(k — k’, t) 3 ke! 


i) 


d 


= Q,(0)45 | ky U(kK—k’,t)N(k’, t) Bk! +i.V KU (ke o} ; 


y . a 
The following abbreviations are used in (21): Se 
- 2kT 
mtvt+m-v—- er 
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This is the standard “‘ambipolar”’ diffusion coefficient, and 


_ Sttanat¥—+ S-tana-W* 


W (6) = : 2) 
uy St+tan at + S- tan a- rote 8) Coal) 


(6) = S*S7 (tan at + tan «) 
us St+tanat+ Stan ar ’ 
2BUS+M-— S-M+ 
P(@)= qo a jee n(25) 
m*yvtm-v-(S* tan at + S- tan «) 


_ S*tanatM- + S~ tan aM 
St+tanat+ S-tana- — 


Q(6) (26) 

Equation (21) can be regarded as differing from (18) in that a common 
velocity of convection, W(@) has been achieved to replace W+, and a 
common diffusive decay as if due to a diffusivity Dy(0). The terms with 
P and @ can still be regarded as describing Dungey’s process, together 
with the intermodulation of different Fourier components, but in practice, 
lack of knowledge of ®(k, ¢) would prevent any exact use of (21) alone. 

W (6) is clearly identical with the ‘wave velocity’ found by Clemmow 
and Johnson (1959) for the linearized theory in the absence of turbulence. 
It can be regarded as a ‘compromise’ drift, between W*+ and W—-; being 
in fact a weighted mean of them, the weight depending on @ and on a+. 
It should not be confused with the W (written without an argument) 
defined in §4, though as we shall see, they are closely related; if it were 
possible for all the contributions to N(k,t) to arise only from those values 
of k corresponding to a fixed value of 6, W would then be just W(@) 
evaluated for that 6. Note also that y(@) is unity for 0=0, and decreases 
monotonically to its value at = $7. 


§ 7. INTEGRAL EXPRESSIONS FOR THE Drirt VELOCITY 


Although we have now set up a complete formal scheme for tackling 
the problem, namely eqns. (4), (9) and (18), it is already clear that even 
if we had enough information about the behaviour of U(k,t), an exact use 
of the equations would be prohibitively complex. But besides this, little 
is known about the temporal behaviour of the velocity field in a turbulent 
flow, even in the absence of buoyancy forces. Because of the non-linearity 
of the Navier-Stokes equation, the equation satisfied by the U(k,?) for 
each realization is of a similar intractable nature, and has not been much 
used (but see Miinch and Wheelon 1958). Clearly, then, we shall have to 
make some approximations or statistical assumptions if we are to make any 
estimate of the overal drift W. For this reason our own estimate of it 
must be regarded as tentative. 

The integrals (9) which we have to evaluate, extend over the whole of 
k-space. Although we shall continue to write them conventionally in 
this way, it should be noted that our assumption discussed in §3 makes 
the contribution from the large scales (small values of k) negligible. This 
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is now seen to be even more plausible on account of the factors k? and 
ko/ot in the integrands. 

As pointed out by Wheelon (1957) the diffusivity D is about twice the 
kinematic viscosity in this context. Hence the convective equilibrium 
range terminates roughly where the inertial one does, and the cutoff 
beyond is almost certainly very rapid (Batchelor 1959). So the dissipation 
range may also be omitted from our integrations as a first approximation. 

By retaining the range of integration as the whole of k-space we are really 
supposing that the temporal behaviour of all Fourier components is of the 
same form as for those in the equilibrium range, on the grounds that those 
parts of the range for which this is not so are unimportant anyway. 

Consider now the first of the integrals (9). To evaluate this we need 
only the ‘spectrum density’ |NV(k,t)|*. This has been considered recently 
by Howells (1960), who concludes that despite the apparent anisotropy 
of the situation, the spectrum is roughly isotropic in the equilibrium range. 
Two cases arise : (i) below about 125 km height, where v* > w*,1.e. tanat <1; 
(ii) above that level, where the opposite holds. Case (i) has been worked 
out in more detail, and one can say with confidence that the spectrum is 
close to isotropic. Even for case (ii) the departure from isotropy appears 
to be fairly small. (The anisotropy of the diffusion is irrelevant here 
because we ignore the range in which any diffusion is important. The 
highly elongated irregularities suggested by Booker (1956) are therefore 
unlikely to occur inthe parts of the ionosphere wherethere is local 
turbulence, though they may occur where irregularities exist for some 
other reason). 

We can, then, write |N(k,é)|?=g(k), say. Using polar angles 0, ¥, 
relative to the 03 axis, (9) gives (where I = k/k) 


(ce) (ar (*Qar 
y= 07} { k4g(|t) die | | ijl; sin 6 d0 cb, 
0 6=0 .! ¥=0 
=48;;, (27) 
where 
Peel LE 
ll k4g(k) dk (28) 


Turning to the second of the integrals (9), we must use the work of § 6 
to estimate 0M (k, f)/0t, and we shall use the form (21) in preference to the 
simultaneous eqns. (18). It is helpful to picture N(k,t), for each fixed k, 
as a locus in the complex plane. Then the first term of (21) isa tendency 
for the ‘vector’ representing N to rotate at a definite angular velocity 
—k.W (4); the second term need not be considered. The two integral 
terms are of a stochastic nature and represent some kind of ‘random 
walk ’ process superposed on the rotation. The last term, if it is important, 
is of a more systematic nature than the integrals, because at fairly small 
scales U,,(k, t) varies nearly systematically—the usual concept of convection 


of the small eddies by the mean or large-scale flow makes U tend to rotate 
with angular frequency —k. %. 
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Using (18) to write eqn. (9) free of time derivatives, we obtain 


bao [k[k.W(? NK Dak, he eee, 015) 
~ Don | xKRIN(K ORR . 2... . . . . (80) 
Hit | ky (h, —by')Py NUON, HOH a? ke 
+o71 i | kky'Q; (0) N*(k, DN(k’, )U,(k— Kt) d? kek’, 


Fo [kk Qu) N0, ‘(k, t) dk. Se eee), 


Here @ refers to the polar angle for the wave number k, not k’, and it is 
convenient to number separately parts of the right-hand side. 

For the part (29), continuing our assumption of isotropy, integration 
over |k| in polar coordinates and making use of the definition (28) of G 
yields a contribution to W(=a-.b) which is 


G 2a 7 
=| | I[L.W(@)]sindody.  . . . . (33) 
Ai v=0/ 0=0 

The part (30) is zero since the integrand is an odd function of k. 

The part (31) shows explicitly the contribution to drift velocity of the 
intermodulation effects. Triple correlations of N, ® and U evaluated at 
three different wave numbers are involved here. We shall assume that 
these correlations are negligible; or in other words that the interactions 
between different wave numbers are entirely random, showing no preferred 
direction in space-time. This is not to say that the non-linear effects are 
unimportant—they are indeed important in determining the form of the 
spectrum, i.e. g(k), but can be ignored in the drift problem. 

Finally, we have to examine part (32). For case (i) of the classification 
mentioned above, i.e. for heights rather less than 125 km, this term ‘does 
provide a systematic correlation, but it is small, because Dungey’s effect~ 
(see Dungey 1956) is itself small in this region. This can be seen explicitly 
in our formulation by means of eqn. (26), which shows @ as a weighted 
mean of M+ and M-. Below about 90km, where v>w for the electrons 
also, both M+ are very nearly the unit matrix. Between 90 and 125km, 
where M* is still nearly unit, M~ is not so, however the factor tan a7 
increases to such an extent that except for @ very close to $7, the mean is 
weighted heavily in favour of M+. Thus it is only near and above 125km 
that Q differs appreciably from the unit matrix, and until that happens 
(32) gives zero since k.U(k,t)~0. Since U(k,t) is isotropic, no special 
weight is attached to the small range of @ in which this argument fails, so 
as a first approximation, we can justifiably neglect (32) for case (i). 
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For case (ii), the position with regard to (32) is much more difficult, but 
some remarks are offered in § 9. 


§ 8. EvaLuaTion oF THE Drirr VELOCITY IN THE SIMPLE CASE 


For case (i), the only surviving contribution to W, according to the 
foregoing arguments, is given by (33). This is no more than a mean over 
all directions of the velocities for the Fourier components, each resolved 
along the ‘wave normal’. Our statistical definition of drift has thus given 
a new significance to Clemmow and Johnson’s ‘ wave velocity ’, originally 
derived by a linearized theory in the absence of turbulence. Note also 
that the distribution among wave numbers, g(k), disappears from the drift 
problem once an isotropic spectrum has been assumed. ‘This is of course 
to be traced back to the fact that the convection for any Fourier component 
k depends solely on @ and not on the wave number k. Since there has been 
some controversy (summarized by Batchelor 1959) concerning g(k) even 
for a non-magnetic case, we can regard this as being rather fortunate. 

The last step is to write 1=(sin@cos%, sin@sing, cos?) and use (20) 
and (23) to carry out the integrations (33); this involves only elementary 
integrations, and a convenient form of the result is 


Wy 2=(1-R)4% 1 ot + RY,. and Wz=(1—Ry)V 47+ BV 5. (34) 
Thus each component of W is a weighted mean of the corresponding 
components of Y+, the weighting factors being R, for those perpendicular 


to the field, and R, for that parallel to it. R, and R, depend solely upon 
a+, i.e. upon height in the ionosphere, according to the formulae 


sin a+ cos aW 
(= tase {1+ tan a tan (a7 — a*)I,[cos at Cos a~cos (a~ — a*) ]} 
(35) 
(for 7=1,3) 
where 3x (72 sind 6dé 
a a ate ca 
2 Jo cos?é+xsin?6 


x 3 l—a |"? 
_— 2 — ——_—_—__. Fi = —— 
eat x—5 RAR emo e | = ] \. (36) 


and I(x) = 2a[1 —I,(x)]. erty or. ae 
For small x one uses 


37 Or 
T(x) x —awl2—4 fee 3/2 
1() a B+ 7 +080... 


The first factor.of (35) can also be written 


sin a+ Cos a— 
sin (at ++ a) = [1 + tan a—/tan at] = [1 +. mtytlm-p-|-! ~ m-v—|mtvt, 

(38) 
80 is approximately independent of height in the ionosphere since y+/p- 
is nearly constant—we shall take it to be 1/30 in the numerical work which 
follows, and m+/m- to be 5-5 x 104. Then Rk, and R, depend only on say 
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tan a~ = w/v, since w* is then fixed. The figure shows the values calculated 
from (35). This can be applied for any atmospheric model if the two 
numerical parameters are as just quoted, but estimated heights are also 
marked on the horizontal axis, taking w-=7-5 x 10sec! and values of 
v~ based on recent atmospheric models (kindly supplied by Dr. K. Weekes). 


-OS 
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-OOOS 


60 70 80 90 100 120 {50 Height 


Illustrating the dependence upon height of the weighting factors R, and R, 
which give the relative importance of the ion-drift and electron-drift 
in determining the drift of the irregularities. 


Consider first the lowest levels, around 70km, where both «+ are small. 
The large bracket factor in (35) is then ~1, (for both j = 1,3) and we have 
Rx Rx [1+mtyt/m-v-}4. But (15) has unit matrices in place of 
M= in this approximation, so that, as expected, equations (33) revert to 
W = &, irrespective of E. 

Slightly higher in the ionosphere this simple result changes, partly 
because (15) becomes more complicated and partly because the coefficients 
R, and R, increase, especially R,. One might have expected intuitively 
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that however V+ vary in this way, the drift of irregularities would always 
remain approximately equal to V+ on account of the much greater mass of 
theions. The coefficients k; measure just the extent to which this must be 
corrected due to the motion of the electrons, an effect reaching a maximum 
in the middle of the region characterized by w— > v~ but wt <v*, i.e. pre- 
cisely the dynamo region. The effect is, however, fairly small unless E and 
M ave so arranged that W~ is somewhat larger than ¥~. It can also be 
interpreted as the dynamical effect of the ‘ polarization ’ of the irregularities 
due to the mean current flow and the fluctuations in conductivity (cf. 
Clemmow and Johnson’s example of a cylinder, 1955). 

Above the dynamo region, the coefficients R, decrease again to their 
former value, but the further complications discussed in the next section 
make this simple theory unreliable. 


§ 9. Toe Upper E REGION AND THE F REGION 


Above the dynamo region, W+ tend to equality, in which case the 

velocity given by (33) is also the same, namely 

Ww t+2v—-x (E#,/B, —E,/B, Us) ee ict) 
the coefficients R, being practically irrelevant. However, as @ is no 
longer nearly the unit matrix, the contribution from part (32) of the 
expression for b must be considered, and it has not proved possible to 
derive any mathematical results. 

Some insight into this situation may be gained by linearizing (21), and 
omitting the diffusion term, for simplicity. We then have a ‘small 
signal’ theory, in which each Fourier component behaves independently 
and Dungey’s effect may be studied in isolation. We have, for each k, 


= +7 W(0). KN = —-1WQ,kU;,. Wey 
Taking U, to be proportional to exp(—ipt) where p=k.@%, which as 
mentioned above is a good approximation for the smaller eddies, the 
solution for V(t) consists of the particular integral proportional to exp (—ipt) 
and a complementary function proportional to exp [—ik.W(8)t]. The 
former gives a finite correlation in both (29) and (32), and the latter as 
before in (29). But the amplitude of the complementary function is 
arbitrary in a linearized theory; its mean equilibrium amplitude in any 
situation to which the linear theory may approximate is determined by 
the non-linear terms. Setting all the complimentary functions zero leads 
of course to W=%. Presumably when both turbulent mixing and 
Dungey’s process are important, W is in some sense intermediate between 
WM and the value found in §8. 
As the height of the maximum electron density in the E region generally 
occurs at 120-130km, no reflections can be observed from this level up to 
something like 180 km, and it seems most unlikely that turbulence occurs 


at or above this latter height. So in practice the difficulty discussed in 
this section is probably rarely serious. 
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The main point of interest in this discussion is that it has emerged that 
there is no unique ‘ velocity of irregularities’ at a given level in the iono- 
sphere. It is necessary to enquire into the details of the process producing 
the irregularity—in the present case we are unable to proceed without 
knowing in greater detail the réles of turbulent mixing and Dungey’s 
process. Similarly our theory for the lower E region assumes turbulent 
mixing as the cause, but the whole question would need to be reconsidered 
if, for example it were proposed that inhomogeneities in the ionizing 
radiation produced irregular electron density. 

In the case of the F region, the difficulty becomes acute, as there is no 
adequate theory of the origin of the observed irregularities. Dagg (1957 a) 
has reviewed several explanations and criticized all of them. His own 
theory (1957b), of small-scale electric fields originating in the dynamo 
region, is also open to criticism, however, as such electric fields would 
not produce compressions in the plasma (Dougherty 1959). Alfvén waves 
travelling along lines of force and originating in the outer ionosphere have 
also been suggested, (Dessler 1958) but the velocity would then be the 
wave velocity, which is too large to agree with the observations. 


§ 10. Discussion or RESULTS 


Let us now reconsider our objective, the interpretation of the drift 
velocities derived from fading records. To do this, we work in new 
axes Ov, Oy and 0z, taken along magnetic East, North, and vertically 
upwards, respectively. Then assuming that the fading is imposed at or 
near the height of reflection, the spaced-receiver experiments yield W, 
and W,, appropriate to that height. 

Below about 90 kms, we can with confidence say that this gives directly 
the wind velocity for the neutral air, Y, and %,. (W, can generally be 
assumed to be an order of magnitude smaller.) This is true even if a 
magnetic field is present. 

In the neighbourhood of the dynamo layer, the electric field plays an 
important part in determining W, and one cannot proceed without 
ambiguity from W,,W, to Y,, Y,. Taking %,=0 and assuming £, 
adjusted to make 7 ,* = VY ,- (no vertical current), we still have two further 
unknown quantities, H,, L,. 

The further relations between Y%,, Y, and H,,H, must of course be 
provided by means of a dynamo theory of their variation over the whole 
globe. So the present theory can be used to find whether a proposed 
dynamo theory successfully predicts the drifts, and for example, their 
diurnal behaviour. 

It is certainly possible to construct numerical examples in which the 
error involved in identifying the components of W as those of @ would 
be very serious. 

For much of the published data (see for example Briggs and Spencer. 
1954) the height to which the drift refers was not measured, except so far 
as to identify the region, E or F, and this makes satisfactory interpretation 
virtually impossible. Recently, howeyer, Jones (1958) has reported 
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results in which the height was estimated ; simultaneous measurements at 
different heights near the dynamo region revealed a gradient in the 
velocity, which Jones expressed as a change of phase in the diurnal com- 
ponent of drift. A possible explanation is that although @ and E are each 
independent of height (as assumed in the dynamo theory) the diurnal 
variations of Y,, U,, H,, H,, have different phases. Since W depends 
increasingly upon E with increasing height, a gradual change in phase of W 
would automatically take place. 

In conclusion, we may claim that a theoretical basis has been laid for 
relating the ionospheric drift observations to the winds and electric fields 
of the dynamo region. 
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ABSTRACT 


Single-crystal films of silver were made by evaporation of silver onto the 
cleavage surface of a heated rock-salt crystal. The crystalline defects 
present in the detached films were studied by transmission electron 
microscopy. In fresh areas of the foil, numerous simple stacking faults 
confined to a single plane were observed. The number varied between 
5x 10° and 1x10!°cm™. There was a strong tendency for the partial 
dislocations to lie in ¢110) directions. The morphology of the simple 
fault shapes resulting is considered in detail. Electron beam stressing 
resulted in the dissociation of whole dislocations and in the widening of 
faults already present. The interaction of the faults to form complex 
faults on several planes, and the effects of holes and dislocations were 
studied. Some areas showed minor defects consisting of dislocation loops 
and tetrahedral stacking faults. Preliminary annealing experiments 
indicate that the simple faults are remarkably stable at 300°c. The origin 
of the various defects is discussed. It is concluded that most simple faults 
probably form by dissociation of dislocations which may (1) grow in from 
the substrate, (2) nucleate at silver patch edges during growth, (3) form by 
impingement of patches on a {110} plane during growth. 


§ 1. InTRODUCTION 


Unttt recently the electron microscope had not been used appreciably to 
study the structure of metal films condensed onto single-crystal substrates. 
It is, of course, well known that single-crystal films of metals may be made 
in this way as a result of ‘epitaxy’. Pashley (1956) has reviewed the 
methods of achieving this. Asa result of the work of Bollmann (1956 a, b), 
Hirsch et al. (1956), Whelan and Hirsch (1957 a, b) and Whelan et al. (1957), 
techniques are available for directly observing crystalline defects such as — 
whole and partial dislocations and stacking faults. 

The present work is an application of the direct observation techniques 
to the study of defects present in evaporated films of silver, both as-grown 
and after stressing in the beam of the electron microscope. This work was. 
carried out independently of the related study of Matthews (1959) on 
evaporated films of gold, silver and copper, and does not overlap greatly. 
Pashley (1959) and Bassett and Pashley (1958-59) have investigated the 
structure of evaporated films, particularly of gold, but their results are 
complementary to the present work. A preliminary account of the present 
work has been given (Phillips 1959). 


+ Communicated by the Author. 
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§ 2. EXPERIMENTAL DETAILS 


Films were prepared by evaporation of silver of 99-999% purity from a 
conical tungsten filament onto a freshly cleaved rock-salt crystal main- 
tained at 300°c. The evaporation was carried out in a vacuum of about 
10-mm of mercury. The films were cooled on the substrate in vacuo to 
room temperature which took about 20min. The silver was floated off 
on a water surface, washed and pieces transferred to specimen grids. 

Single crystal foils with an {001} plane parallel to the surface were 
obtained. The thickness of most foils ranged from 400-7004, although 
in the immediate vicinity of holes they were, of course, thinner. 

Specimens were examined in a Siemens ElmiskopI electron microscope 
operated at 100kv. Double condenser illumination was used. The 
samples were tilted by means of a stereo-holder to bring any defects present 
into contrast. A ciné film was made of the effects observable at about 
x 40 000 on the fluorescent screen. The improved resolution obtained 
by recording the image directly on plates was necessary for detailed study. 


§ 3. Imagine METHODS 


Stacking faults in {001} films are readily seen by diffraction contrast and 
typically show extinction fringes parallel to the intersection of the {111} 
fault plane with the foil surfaces, that is, parallel to ¢110) directions in 
the plane of the film. The detailed theory of these fringe effects has been 
worked out by Whelan and Hirsch (1957a,b). Figure 1 (schematic) 
depicts the type of electron diffraction pattern given by an {001} silver film. 
The ‘crosses’ such as that enclosed by circle B are a well-known feature 
of such patterns. (See for example Watanabe 1956, Ogawa et al. 1955.) 

In the present work stacking faults were examined both in bright field 
and in dark field. The respective objective aperture positions for these 
two methods are shown by the circles A and B in fig. 1. The objective 
aperture size was 50. Dark field images were obtained using the portion 
of the diffraction pattern encircled at B comprising a {200} reflection and 
satellite spots. The gun was tilted to bring an {002}! onto the axis of the 
microscope. Under these conditions, individual faults contributing to 
this portion of the diffraction pattern are brightly imaged and individual 
fringes are not usually discernible. Individual dislocations may be seen 
in dark contrast where there is an overlapping stacking fault. Improved 
resolution of individual faults at moderate magnifications appears to be 
obtained in dark field. In bright field individual faults may be imaged 
in dark contrast without resolving individual fringes. This usually 
occurs in the vicinity of a diffuse tilt fringe as in fig. 2 and some portions 
of fig. 4(b). The stacking fault contrast here is probably dynamical in 
origin. More commonly, in bright field, individual fringes are resolved at 
stacking faults in the way predicted by theory (Whelan 1958-59). 


{ Figures 2 to 9 are shown as plates. 
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§ 4. STRUCTURE OF THE Fins 


Figures 2-9 are electron micrographs of the films. Stacking faults were 
observed in all areas. The number of stacking faults varied from one film 
to another but was usually between 5x 109 and 1x 10!percm?. This 
number is relatively independent of the thickness since most faults ran 
from surface to surface. There was a big difference in the length and 
complexity of the faults depending on whether or not a particular area of 
foil had been previously exposed to the electron beam. Exposure to the 
beam caused buckling, as evidenced by movement of tilt fringes across the 
field, and the stresses produced caused wide dissociation and frequent 
combination, and interaction of stacking faults both with each other and 
with dislocations and holes through the foil. These effects will be dis- 
cussed later in more detail. Figures 2 and 4 are typical of the structure 
observed in areas of foil not previously irradiated. Figure 5 is typical 
of the appearance after exposing to the beam for some time. 
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Typical electron diffraction pattern from a silver film (schematic). Circles A 
and B show the portions of the diffraction pattern used in the light and 
dark field imaging methods. 


Defects which appear to be typical dislocation loops were observed in 
fresh areas of the films as at A in fig. 2. These loops ranged in diameter - 
from about 100-3004. Smaller loops would not be recognizable but may 
well be present. The average number of loops in the region shown, which 
is about 6504 thick, is approximately 2 x 10‘4cm~*. Assuming an average 
loop diameter of 2004, in order to calculate the total line length per unit 
volume, p, this loop number is equivalent to a dislocation density of p/8, 
that is, 4x 108em-*. Loops of similar sizes have been observed in quenched 
aluminium by Hirsch ef al. (1958). If the loop is assumed to be a collapsed 
disk of vacancies, then for an average diameter of 200A, the corresponding 
initial concentration of vacancies is 15x 10-5. Dislocation loops have | 
been observed in quenched bulk silver by Smallman and Westmacott 
(1959). 
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A number of parallel line features will be noted in fig. 2, for example 
B,_,. These are not due to slip contrast since this appeared to fade in less 
than 1/10sec. They are instead attributed to movements of loops during 
the 3sec period of the exposure. The projected direction of movement is 
parallel to a (110) direction in the plane of the film. Loops B, and B, 
(fig. 2) appear to have slipped out to holes through the foil. B; appears 
to have moved in a jerky fashion and is still visibile. If this interpretation 
is correct, it implies that the loops are prismatic loops, formed by nucleating 
a Shockley partial in the disk of stacking fault associated with the Frank 
sessile produced by collapse of a disk of vacancies on a {111} plane 
(Kuhlmann-Wilsdorf 1958). These would glide on a cylinder whose axis 
is parallel to a (110) direction not lying in the plane of the fault. Opposite 
halves of the loop would presumably glide in opposite directions as observed 
by Hirsch et al. (1958) in aluminium. Some doubt is thrown on this inter- 
pretation by the fact that the stacking-fault energy in silver is estimated 
(Thornton and Hirsch 1958) to be only about 35 erg em? which is hardly 
sufficient to make the nucleation of a Shockley partial energetically favour- 
able unless stress-assisted. The glissile nature of the loops would explain 
why they are only observed on first examination of an area, and why in 
similar observations on evaporated gold (Bassett and Pashley 1958-59), 
they apparently vanished between exposures. 

Perfect dislocations were commonly seen running obliquely through the 
film as at D in fig. 3. Dislocations running perpendicularly to the film 
from surface to surface would appear as black dots in bright field as for 
example at C in fig. 3. Many dot-like features are also visible in fig. 6, 
some of which are superimposed on stacking faults. Since there is no 
interaction with the faults they cannot be line dislocations running right 
through the film but must be small defects away from the plane of the 
fault, possibly small loops of less than 504 diameter formed by vacancy 
condensation. 

Tetrahedral stacking faults were occasionally observed confirming the 
finding of Matthews (1959). Smallman and Westmacott (1959) claim to 
have observed tetrahedra in quenched and aged bulk silver. Examples 
of tetrahedra are visible at E in fig. 3: the largest has a side length of 
about 500 A, and the smallest about 1504. This agrees well with the 
range of sizes of tetrahedra observed by Silcox and Hirsh (1959) in 
quenched gold, which is to be expected since the stacking fault energy of 
gold and silver are about the same. 

Small triangular defects such as F in fig. 3 and G in fig. 4 (a) were quite 
commonly observed. These correspond to equilateral triangles in a 
{111} plane, the edges of which are parallel to (110) directions lying in 
that plane. The smallest defects recognizable as of this type have side 
lengths of about 100A. The largest would have sizes limited by the 
thickness of the film; for example, that at A’ in fig. 4 (a) has a side length 
of about 650 A. 

The evidence is compatible with the view that the small triangular 
defects which do not in general intersect the film surface are Frank sessile 
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dislocations (Frank 1949) formed by a vacancy condensation mechanism, 
whereas the large defects which generally intersect the film surfaces and 
of which A’ (fig. 4 (a)) is one sort are of a different type. These large defects 
will be discussed in §§5 and 6. The triangular shape of the small defects is 
expected, since Kuhlman-Wilsdorf (1958) has pointed out that a Frank 


sessile dislocation will prefer to split up into segments parallel to (110) 
directions. 


§ 5. DESCRIPTION OF SIMPLE STACKING Favu.Lts 


In areas of the films not previously exposed to the beam most of the 
faults were relatively short and confined to a single plane as in figs. 2 and 4 
Faults confined to a single plane will be called ‘simple faults’. Complex 
faults forming on two or more of the four possible {111} planes were usually 
only observed after irradiation in the beam. Matthews (1959) mainly 
observed complex faults and has described these rather fully. The nature 
of simple faults will now be described. 

Most of the simple faults observed, with the exception of those ending at 
a hole, were straight-sided with 3-6 sides and extended from surface to 
surface of the film. About 5°% of the partials assumed the shortest possible 
length parallel to a {112) direction so that the fault viewed in a direction 
perpendicular to the foil appeared ‘square-ended’. A strong tendency 
was observed for the partial dislocations to align themselves in 110) 
directions as predicted by Kuhlmann-Wilsdorf (1958). Analysis of the 
directions of about 50 partials in figs. 2 and 3 showed that about 60% of 
them coincided with a possible [110] direction within the accuracy of 
measurement, and about 90% were within 15°. The direction of signi- 
ficant deviations was always such as to shorten the length of the partial. 

The stacking faults lie on any of the four {111} planes making a 54° 44’ 
angle with the plane of the foil. The latter will be taken as (001), so that the 
(111) and (111) planes intersect the surface along [110] while the (111) and 
(111) planes intersect it along [110]. It is not possible from the present 
micrographs to determine which of the fault boundaries is in the top surface 
of the foil, so that the stacking fault at T’ in fig. 2, for example, is either on 
the (111) or on the (111) plane if the long edge of the figure is parallel to 
[110]. 
Line drawings of the main possible shapes of simple fault projected onto 
the plane of the film are shown in figs. 10-13. The shapes have been 
classified according to the number of sides and the side direction. The con- 
vention adopted by Matthews (1959) is followed, viz. the intersections of 
the fault with the top and bottom surfaces of the film are shown by thick 
full and dashed lines, respectively. Fine full lines represent the extinction 
fringes caused by the fault. Dashed fine lines represent partial dislocations. 
In each of the four {111} fault planes there are three possible (110) direc- 
tions, one of which lies in the plane of the foil. Thus in the (1i1) plane, 
the partial may be in the [011] or [101] directions or may take the shortest 
distance through the film in the [112] direction. The intercept of the 
(111) plane with the surfaces and also the extinction fringe direction are 
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parallel to [110]. The projection of either [011] or [101] perpendicular 
to the foil makes an angle of 45° with the [110] direction, while the projec- 
tion of [112] is at 90° to [110]. 

Faults of type A’ may occur in two ways in each of the four planes, that 
is in eight ways. Similar information for all the fault types is given in 
brackets in figs. 10-13. Since the top and bottom fault boundaries are 
indistinguishable in the micrographs only a half of these possibilities would 
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Possible shapes of simple fault with partials in (110) directions. 


be distinguishable. Examples of 19 types of fault are identified in 
fig. 4 (7) and may be seen in fig. 4 (b) in light field. Examples of three more 
types are shown in fig. 2. Examples of G’, M’, N’, 0’, P’, V’, W’, AB’, 
AC’ and AF’ are not shown. Of these, G’ and AB’ have been seen in other 
areas. ‘The remaining types would not be expected to be common because 
partials parallel to [112] or its equivalent are relatively infrequent. Since 
the partials tend to deviate somewhat from (110) directions, many of the 
shapes deviate somewhat from those described. 
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Although there was a strong tendency for the partials to be straight, they 
were occasionally curved in a manner more typical of stainless steel 
(Whelan ef al. 1957) as at fault T in fig. 6 (c). 

Many of these simple fault types are also to be seen in a micrograph of 
gold grown directly on rock-salt published by Bassett and Pashley 
(1958-59). Itis probable that they are typical of face centred cubic metals 
of low stacking fault energy grown epitaxially on rock-salt. 


Fig. 11 
48°10) 83°60' 4-48°10' nae 
(111) a Nea as 
y, = Py A Sa a, 
[iid <4 bi. Bee so 
4I°50' ;.! § 
fot] [ii2] (io!) 
J’ 3SIDES K' 3 SIDES L' 4 SIDES 
(16) (8) (16) 
A PN, east 
on CEE P| ‘eraes, 
SS 2 eS SS) cc ee aes 
M' 4 SIDES N' 4 SIDES 0'5 SIDES 
(16) (16) (8) 
aes aE (cates tte ae 
Qa Carrere (Soe ey 
— OF Sey 
pi y| oe a BS Ren 
P' 5 SIDES Q' 5 SIDES R' 5 SIDES 
(8) (32) (32) 
of Te apt A EE 
Nees sere ara pe) 
LTE, Eee ea ae BS nee ae | 
S' 5 SIDES T' 6 SIDES U' 6 SIDES 
(32) (16) (16) 
eda aN 
+———_7 . aa 
ae et ey a | 
Meee Ramaanya Eta 
v' 6 SIDES w' 6 SIDES 
(16) (16) 


Possible shapes of simple fault with one partial in <110> directions, 
other in <112> or <110> directions. 

Comparison of fig. 4(a) with 4(6) shows that most, but not all faults 
are in contrast in fig. 4(a). Some changes have occurred between the two 
exposures, for example at Land M. These will be discussed in §6. The 
object I in fig. 4 is probably a twin. This is one of the very few examples 
observed. 
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H, in fig. 4 (a) is an example of a group of faults on parallel planes spaced 
200-5004 apart. H, appears to be a similar group of three faults more 
closely spaced. J, in fig. 4(a) appears at first sight to be a similar group 
containing three faults, one of which is annular. However, in bright 
field (fig. 4 (b)) the contrast is such that this interpretation appears unlikely. 
J, in fig. 4(a) and (6) is a similar example. The nature of these defects is 
uncertain. 
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Possible shapes of simple faults with both partials in (112) 
or partly in {112) and partly in (110) directions. 


The smallest defects such as K imaged in fig. 4(a) are about 504 in 
diameter. Comparison with fig. 4(b) shows that most of them are located 
at the edges of holes. These are thought to be tiny faults as discussed in 
$6. 

It is possible that simple planar defects are common in bulk material 
but have not yet been recognized. Thus the triangular defects observed 
in quenched and aged bulk silver by Smallman and Westmacott (1959) 
and Smallman et al. (1959-60) appear to be of type A’ although interpreted 
by the authors as defect tetrahedra. Similarly, some of the defects 
observed by Silcox and Hirsch (1959) in quenched and aged gold could be 
simple planar faults of types A’ and ©’ rather than tetrahedra. 


§ 6. Dissociation AND ForMATION OF ComMPLEX Fauuts 


Whole dislocations lying in {111} planes tended to dissociate into two 
(Shockley) partials parallel to (110) directions and connected by stacking 
fault. Dislocations Q, and Q, in fig. 6 illustrate this and also show how 
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this dissociation tends to be inhibited by the stress field associated with the 
dislocations R,_, piled-up on a neighbouring parallel plane against a group 
of stacking faults. The maximum dissociation of Q, is about 200A in 
(a), 100A in (6) and 4504 in (c). The dissociation of Q, in fig. 6 (6) has 
produced a defect of type C’ (fig. 10). Other examples of the blocking of 
the movement of whole dislocations by stacking faults are seen at Z in 
fig. 6 and at Y where the dislocation has moved out of the field between 
the exposures in fig. 6(b) and (c). 

Some individual dislocations such as D in fig. 3 showed no resolvable 
dissociation either because they did not lie in a {111} plane or because 
dissociation was inhibited by the local stress field as at Ry in fig. 6 (bd). 
Dissociation to a width of less than about 100A would not be detectable. 

Stressing in the beam caused the short simple faults typically seen in 
figs. 2 and 4 to widen and interact. Holes, dislocations and other faults 
acted as barriers to dissociation. Eventually most faults have extended 
as far as possible as in fig. 5 and dislocation movements become less frequent. 
For an average (extended) fault width of say 0-3, and assuming 101° 
faults/em?, an order of magnitude calculation shows that this corresponds 
to a plastic strain of about 5x 10-8. Figure 6 illustrates the structure at 
an intermediate stage. 

Dislocation movements were very readily produced in the foil. This is 
attributed partly to the cleanliness of the surface, for it is probable that 
oxide or other contaminant layers tend to impede dislocation movement, 
and partly to the presence of holes which acted as stress concentrators. 
Holes also acted as pinning points for extended faults as at P in fig. 6 (a) 
where one partial has run out to the hole surface giving a fault of type AG’ 
in fig. 13. Other fault shapes such as AI’ in fig. 13 may be observed. 
Exceptionally a fault may run into two holes eliminating both partials 
giving a fault of type AH’ in fig. 13. Thin places in the vicinity of holes 
or elsewhere can similarly reduce the length of both partials and stair-rod 
dislocations. Movement of such partial dislocations would cause lengthen- 
ing and is thus made difficult. Effects of this kind may stabilize normally 
unstable configurations of faults. 

The widening of a fault under stress occurs when one partial is pinned 
and the other moved under the action of the stress. Since movement of a. 
partial produces a shear, the stress may be locally relieved and movement 
stops, however, the shear produces a redistribution of the complex local 
stresses so that another fault is affected and expansion or contraction of this 
fault occurs. Thus movements tended to be rapid but to occur in stages, 
for example, fault S in fig. 6 (a) has extended a short distance in fig. 6 (0) 
and in fig. 6(c) has extended to a hole where the partial will tend to be 
anchored due to the reduction in its length. 

Combination reactions leading to complex stacking faults have been 
discussed by a number of authors (Cottrell 1952, Thompson 1953, Frank 
and Nicholas 1953, Friedel 1955, Whelan 1959). The Burgers vectors of 
the stable stair-rod dislocations have been listed by Frank and Nicholas 
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(1953). While such considerations are probably applicable to fully annealed 
clean films of uniform thickness, they are of doubtful value here because 
it is clear that the reduction in line energy of both partial and stair-rod 
dislocations at holes and thin places plays an important role, as also does 
the blocking of the movement of partials by other faults and dislocations. 
The formation of complex faults as a result of local stressing in the electron 
beam will now be considered. 

A common reaction observed was the combination of two simple faults 
to form a one-bend fault. This has occurred at L and M between the 
consecutive exposures shown in fig. 4(a) and (6). The angle between the 
faulted planes at L is acute (70°32’) and at M is obtuse (109° 28’). The 
stair-rod dislocations (Thompson 1953) at the bends cannot here be resolved. 
Examples of stair-rods may be seen in fig. 6 (a) at the obtuse angled fault 
S and at the acute angled fault y. Partial dislocations are sometimes 
visible, thus both partials are visible in fig. 6 (6) at the single bend fault W. 

Another example of the formation of a single-bend acute-angled fault 
from two simple faults is shown in a ciné sequence (fig. 7 (a) to (c)). 
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Possible shapes of simple faults running into holes—one or no partials. 


The fact that the triangular fault at L in fig. 4(a@) combined with the 
other fault would suggest that it is bounded by Shockley rather than 
Frank partials lending some support to the distinction of two types of 
' triangular fault in § 5. 

The nucleation of dislocations at thin foil edges in stainless steel under 
the action of thermal stresses has been observed by Whelan et al. (1957). 
They calculate, following a private communication from Cottrell, that it is 
relatively easy for a dislocation to nucleate at a thin edge ifit is of edge type. 
It was commonly observed in silver that faults tended to nucleate from the 
edges of holes through the films. The tiny faults K (fig. 4(a)) which are 
located at hole edges (fig. 4(b)) are thought to represent a very early stage 
in this process. These faults are only 50-1004 in diameter. Figure 6 (a) 
to (c) shows the apparent nucleation of three successive faults from the 
hole edge at point W. The first and last faults to form are on parallel 
planes as indicated by the change in fringe contrast. A further example 
of simple overlapping faults is seen;at X in fig. 6(b); these are seen to have 
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lengthened in fig. 6(c). Another instance of fault nucleation at a hole is 
shown at N in fig. 9. The source of the dislocations Q, and Q, in fig. 6 is 
probably the hole edge to the left of Q,. 

The behaviour is frequently more complex. Comparison of fig. 8 (a) 
and (b) shows that the single-bend fault DEF has become a complex double- 
bend fault DEHGFI. 

Simple faults were frequently observed to become shorter as at PQ in 
fig. 9(b) to (c). Complex faults were frequently observed to contract 
forming simpler faults. Examples of this behaviour are seen at UT in 
fig. 6 (a) to (c) and at JK LM in fig. 9 (a) to (c). 

Preliminary experiments have been made in which films were annealed 
mm vacuo on grids for half an hour at 300°c followed by slow cooling. It 
appears that this treatment reduced the number of faults by a factor of 2-3. 
After annealing most of the faults observed are short simple faults. Whole 
dislocations are also present. Due to a loss in resolution occasioned by 
surface deterioration during annealing it was not possible to study the 
minor defects. Electron beam stressing caused the simple faults to widen 
and interact to form complex faults as before. This annealing experiment 
emphasizes the remarkable stability of the simple faults. 

Widening and interaction processes of the kind described create complex 
fault patterns which impede free dislocation movement. These processes 
may play an important role in work hardening in bulk materials such as 
silver and gold if simple faults are present which may be the case as suggested 
in §6. If such defects are peculiar to evaporated films then they may 
account for the claims that the strengths of such films may considerably 
exceed that of bulk material. 


§ 7. ORIGIN oF INITIAL DEFECTS 


The formation of the minor defects present namely dislocation loops, 
stacking-fault tetrahedra and small triangular defects is attributed to a 
vacancy collapse mechanism as discussed in §4. The dot-like features are 
believed to be some kind of vacancy aggregate, perhaps unresolved loops. 
A vacancy concentration of about 5 x 10-° is sufficient to account for the 
formation of all these defects. This seems a reasonable value for an 
evaporated foil. : 

The major defects present, namely stacking faults with their associated 
dislocations, could have formed by a number of alternative processes : 


I. Matthews (1959) growth impingement mechanism. 
II. Vacancy collapse mechanism. 
III. Growth mistakes. 
IV. Dissociation of whole dislocations. 


These will now be considered in turn. 

Matthews supposes that adjoining nuclei of silver are of different heights, 
grow with {111} facets, and impinge on a {111} plane with a mis-match in 
the stacking sequence. By suitable elaboration, the formation of both 
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simple and complex faults can be explained. If this mechanism were 
important, facetting would be evident at the edges of the holes through the 
present films. This is not observed. A further objection to Matthews 
mechanism is that growth of nuclei in the case of evaporated gold (Bassett 
and Pashley 1958-59) and presumably in silver, seems in the early stages to 
occur by the growth of some nuclei at the expense of others rather than by 
impingement and coalescense, yet stacking faults appear to form at a very 
early stage in the growth before impingement becomes common. 

A vacancy collapse mechanism will now be considered. In a fresh area 
of film, such as that shown in fig. 4, the density of major faults is about 
1x 10!%cem-2. Assuming an average width of 0-1 which is fairly typical 
of a fresh area and assuming that the faults formed by vacancy aggregation 
and collapse, this density is equivalent to an initial vacancy concentration 
of about 3x 10-3. This seems unreasonably high. Furthermore, most 
of the partial dislocations appear to be glissile in the plane of the associated 
fault. These considerations appear to rule out a vacancy collapse 
mechanism. 

Growth mistakes could occur at an early stage during growth, possibly 
in laying down the first few {111} layers, which could account for the 
presence of faults running from surface to surface of the films. The low 
stacking fault energy in silver would increase the probability of such 
mistakes occuring. More evidence is needed in order to properly assess 
the importance of this mechanism. 

The most important mechanism of formation of stacking faults is thought 
to be the dislocation dissociation one. Dislocations might be nucleated 
at silver patch edges during growth since it is probable that the silver is 
in a state of stress. Secondly, they could form by the impingement of 
adjacent patches on a {110} plane during the later stages of growth. 
Thirdly, it is also likely that dislocations grow from the substrate. Thus 
edge dislocations with a (110) Burgers vector emerging on a {100} cleavage 
surface of the rocksalt could grow into the silver during the first stage of 
deposition, and as growth continued would tend to dissociate into the 
appropriate Shockley partials. Since there is a marked tendency for the 
partials to lie in (110) directions, it would be expected that faults formed 
in these ways would commonly tend to be triangular in shape (type A’). 
The driving force for movement of the partials into (110) directions pre- 
sumably comes from a reduction of core energy in this orientation. Since 
the stacking fault energy is low, only a small stress is then required to cause 
further broadening of type A’ faults to form type C’. Considerable 
stresses arise during cooling of the film on the substrate due to differential 
contraction and also during stripping and drying due to surface tension 
effects. If both partials lay in the same (110) direction a fault of type B’ 
would be expected and could readily widen under stress. Faults such as 
type AG’ could result from the nucleation of a single Shockley partial at a 
patch edge. As demonstrated, further dissociation and interaction to form 
complex faults can readily occur. New faults can also be nucleated, 
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ABSTRACT 


The solar diurnal variation of the cosmic ray meson intensity caused by 
the periodic solar heating of the atmosphere is derived from a comparison of 
the ionization chamber and neutron monitor data from Huaneayo. The 
amplitude and the time of maximum are found to be 0:11% and 05-30 hr 
respectively. Recent meteorological data, free from appreciable radiation 
errors, are found to be in rough agreement with this result. Particularly good 
agreement is obtained in the phases which are both later than the value given 
by Dorman by about 4hr. The inclusion of this temperature correction 
enables the directional measurements at Mawson to be reconciled with the 
idea of a primary anisotropy. 


§ 1. INTRODUCTION 


TuE solar daily variation of cosmic ray intensity has been the subject of 
extensive investigation, chiefly because of the important information it 
ean yield about the anisotropic distributions of the primary radiation 
and the modulating mechanisms responsible for their production. Con- 
tinuous monitoring of the intensity, however, has been feasible only with 
the secondary radiation at sea-level or mountain altitudes where measure- 
ments have been made with ionization chambers, counter telescopes and 
neutron monitors. The intensity variations of the secondary radiation, 
however, are governed not only by the primary intensity variations but 
also by the variable conditions of the overlying atmosphere, and adequate 
corrections have to be applied for such atmospheric influences in order to 
relate the observed variations to those in the primary intensity. 

The neutron-producing component of the secondary radiation is sensitive 
to changes of barometric pressure but is not affected to any appreciable 
extent by atmospheric temperature. Since effects of barometric changes 
can be taken adequately into account, observations of the neutron inten- 
sity provide reliable measures of the solar diurnal variations of the primary 
radiation entering the earth’s atmosphere. In contrast, the total ionizing 
and y-meson components measured by ionization chambers and counter 
telescopes are influenced both by pressure variations and variations of 
temperature throughout the atmosphere. As will be seen later, difficulties 
are encountered in attempting to correct for the atmospheric temperature 


effects. 
+ Communicated by the Authors. 
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In spite of these difficulties over meteorological effects, counter telescopes 
possess an advantage over neutron monitors in their ability to select the 
radiation incident from particular directions and because of this, are of 
particular importance in a determination of the true direction and magnitude 
of the anisotropy existing outside the earth’s magnetic field. For instance, 
Brunberg and Dattner (1954) have shown how, by making some simple 
assumptions, the mean energy of the primary radiation giving rise to the 
solar daily variations and the direction of the anisotropy relative to the 
sun-earth line can be deduced from the directional measurements made 
with telescopes inclined at 30° to the vertical and pointing in the north 
and south directions. 

Measurements of this kind made by Elliot and Dolbear (1951) and 
Malmfors (1949) have shown that the observed solar daily variations were 
of about the same magnitudes as those originating in atmospheric tempera- 
ture variations. Under these circumstances, it is necessary to know the 
meteorological corrections with considerable accuracy in order to obtain 
any reliable information concerning the anisotropy. 

Now, although many elegant theoretical investigations have already 
contributed to an understanding of atmospheric effects on the ~-meson 
intensity near sea-level (Olbert 1953, Maeda and Wada 1954, Trefall 
1955, Dorman 1957), corrections for meteorological influences on the 
solar daily variations of meson intensity have not hitherto been evaluated 
to a high degree of accuracy. The chief difficulties encountered have been - 
a dearth of aerological information and where radio-sonde measurements 
of atmospheric temperature were available, the uncertainty as to the 
extent to which the measurements were contaminated by insolation and 
lag errors. 

There is general agreement that the radio-sonde measurements of the 
diurnal variation at high altitudes do not reflect the actual diurnal course 
of the air temperature (Kay 1951, Dorman 1957). However, Dorman has 
estimated the atmospheric contribution to the diurnal variation of meson 
intensity on the assumption that the radio-sonde measurements were 
reliable up to the 300 mb pressure level. The diurnal amplitude in the 
region above this isobar was estimated from cosmic ray measurements 
made underground at London to be <0-5°c while the upper air soundings 
gave amplitudes of over 2°c. When errors of over 300° occur just above 
300 mb, the limit of integration, it is conceivable that substantial radiation 
errors are also present in the measurements made below 300 mb and to 
this extent, the estimate of Dorman appears to be unreliable. An alter- 
native estimate of the atmospheric contribution is therefore desirable. 

In the following, we have attempted to deduce the 24 hr wave in meson 
intensity due to atmospheric temperature using only cosmic ray measure- 
ments and compare this result with the most recent meteorological data. 
Finally, it is shown that when this temperature correction is included, 
the directional measurements made at Mawson can be explained in terms 
of a primary anisotropy without any inconsistencies. 
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§ 2. Comparison oF NEuTRON Montror anv IonIzATION CHAMBER 
Data FROM HuaNncayo 


In this section the diurnal variation of meson intensity arising from 
changes of atmospheric temperature and hereafter referred to as the 
temperature wave, will be derived from a comparison of the ionization 
chamber and neutron monitor data from the equatorial station Huancayo. 
Tn this particular case where the response of the two detectors to primary 
intensity variations is likely to be similar, the basic cause of the difference 
between the two sets of results, after correction for barometric changes, 
is the residual atmospheric contribution left in the ion-chamber data. 

We shall first examine the basis for the assumption of similarity in 
detector response. Rough similarities are evident in (i) Directional 
Sensitivity and (ii) Energy Response as defined by the differential response 
curves giving the number of detected secondary particles produced, as a 
function of the primary magnetic rigidity. These will now be described 
briefly under (i) and (ii) 

(i) Both instruments are omni-directional in the detection of the 
secondary radiation and consequently the directional sensitivities to 
the primary radiation are determined solely by atmospheric 
absorption. , 


(ii) It is evident from the differential response curves given by Webber 
and Quenby (1959) that at rigidities higher than the equatorial 
cut-off value, the curves for an ion-chamber and a neutron monitor 
are similar at the atmospheric depth of Huancayo. 


The following experimental result lends further support to the belief 
that the instrumental responses are similar. It is known that during 
certain periods, particularly after magnetic and cosmic-ray storms, the 
amplitude of the daily variation may be as great as 0-5°%, whereas the 
temperature wave in the ion-chamber data is likely to be only 0-1% in 
amplitude. Ifthe response of the two instruments to primary variations 
is similar, it follows that the difference between the daily variations 
measured by neutron monitors and meson recorders during these periods 
should be small. 

Now, neutron monitors and wide-angled counter telescopes have 
concurrently been operated during 1957 and 1958 at the low-latitude 
stations Lae and Tokyo. The averaged results from several selected 
days of enhanced daily variation which followed Forbush decreases were 
subject to Fourier analysis and the 24 hr components are shown on harmonic 
dials in fig. 1. It is apparent from the diagram that the measured ampli- 
tudes are ~ 0-6°% while the vectorial difference between the two detectors 
at any one station is <~0-1%. Although the data from a neutron 
monitor and a counter telescope have been compared here, this result 
should also apply to the data from a neutron monitor and ion chamber at 
a low-latitude station because the differences between the differential 
response curves of an ionization chamber and a counter telescope are 


2Q2 
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negligible (Webber and Quenby 1959). It is concluded therefore, that 
the two types of detector at Huancayo respond in similar ways to primary 
intensity variations. The evidence cited above appears to be sufficient 
as only a qualitative similarity and not an exact correspondence in instru- 
mental response is required for our purpose. 

Let the mean solar 24hr waves during any given month originating 
from primary intensity changes be represented vectorially by N and I for 
the measurements made by the neutron monitor and the ion chamber 


Fig. 1 


Ob6 he’ LT. OG br. eC 


(2)\he LT, I2hr. L.T. 


Harmonic dials, showing the 24 hr components of the daily variation measured 
during days of enhanced variation. Simultaneous measurements made 
by a neutron monitor (V) and a wide-angled meson telescope (J) at 
Lae and Tokyo are shown. 


respectively. It will be assumed that the similarity of instrumental 
response leads to the following relationships : 

IN|=rll|, 

N .1=6|N]. |I| 
where r and @ are constants even though the rigidity dependence and the 
direction in space of the primary anisotropy may change in time. These 


conditions will be satisfied best if, as implied earlier, r and @ are both 
close to unity. 
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Since the data have already been corrected for barometric changes, the 
diurnal waves given by the measurements should be N for the neutron 
monitor and 1+ T where T is the atmospheric contribution to the meson 
data. Moreover, for the ion-chamber, the seasonal variations of the 
temperature wave are likely to be small at Huancayo which is near the 
geographical equator. Thus for practical purposes T can be assumed 
constant. 

If we now select a period of time during which the mean diurnal variation 
measured by the neutron monitor is zero, it can be shown that the average 
variation measured by the ionization chamber over this period is the 
temperature wave T; i.e. if $}N,;=0 where 7 denotes the monthly period, 


then > 1,=0, provided that r and @ are constant. This is easily seen when 


we consider a number of vectors N, selected so that on addition they give 
a closed polygon. Even when the amplitudes of the vectors are changed 
uniformly by a factor r and each vector is rotated through the same angle 
cos~'@, the polygon will be closed. Therefore the mean 24hr wave in the 
ion-chamber data is given by 


Ms 


(l,+7)=T. 


Sle 


a 


I 


1 
The condition }N,;=0 can be realized using the data from the period 


1953-1955. The neutron monitor data of Conforto and Simpson (1957) 
have been used in our analysis together with the results of Venkatesan 
and Dattner (1959) after allowing for some errors in the latter analysis. 
The definition adopted for the diurnal vector in the former publication, 
though different from the conventional one, is claimed by the authors to 
give equivalent results. 

The selected groups of months during which the mean neutron 
diurnal effect was nearly zero are shown in table 1 together with the actual 
mean values and the corresponding temperature waves from the ion- 
chamber data. The coefficients a, and 0, are given, where the diurnal 
wave is represented by 


7 t+b,sin 57. 

The atmospheric vectors are plotted on a harmonic dial in fig. 2 where © 
the standard deviation circles are also drawn. The vectors for each group 
are numbered 1 to 6 and the overall mean is marked M. The individual 
vectors will be seen to agree within the limits of statistical accuracy but 
we note that they are not entirely independent of one another because 
some months are common to several groups. The mean temperature 
wave obtained in this way has an amplitude of 0-11% and a time of 
maximum around 05-40 hr local time. 

Alternatively, r and @ can be determined from the data for two separate 
periods, again assuming 7, 6 and T to be constant. In order to achieve 
accuracy in this calculation, the measured diurnal waves during the two 
periods should show substantial differences. The mean diurnal waves 


a, COs 
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Table 1 


Mean neutron Atmospheric 
Group Selected months vector (%) vector (%) 
ay b, ay b, 


1 October, 53 Feb., March, April, 
May, June, July 1954 — 0-002 0-01 0-027 — 0-14 
2 March, April, September, Octo- 
ber, November, December 54 | —0-008 —0O-004 0:00 0-12 
3 September, October, November, 

December 53; January, Feb- 

ruary, March, April, August, 


September 54 —0-001  —0-005 0-01 © 0-12 
4 April, May, June, July, 53, 
August, September 54 =—U-O1 +0-01 0-02. 0-10 


5 March, April, May, June, August 
September, October, Novem- 
ber, December, 54, January, 


February, March, 55. —0-004 +0-002 0-00 O-11 
6 January 54 to March 55 
inclusive 0-00 0-00 0-02 0-10 


(+0-01, +0-11) 


OOfhr. L.T. 


0.10% 


0.05%. 


Solar time dial showing the estimated 24 hr component of the daily variation 
due to atmospheric temperature in the ion-chamber data from 
Huancayo. Vectors 1 to 6 (as in table 1) and the overall mean M were 
found by the null method. The vectors (54-54) and A were derived by 


alternative methods (vide text). The circles denote standard deviations 
of the estimates. 
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for the years 1953 and 1954 satisfy this criterion and these data yield 
values r=0-96 and 6=0-94, consistent with our earlier assumption r, O~ 1, 
Using these values of r and 6, the temperature wave in the Huancayo 
ion-chamber is found to have harmonic coefficients a,=0-026% and 
b, =0-102%. 

A further estimate of the temperature effect can be obtained by assuming 
that the primary variations produce 24 hr waves of identical magnitude 
and phase in the neutron monitor and the ionization chamber, or r=4=1. 
Here, the vector difference between the ionization chamber and neutron 
monitor data gives the temperature effect and the measurements made 
over the years 1953 to 1955 give values a,=0-035% and b,=0-106% for 
this vector. The estimates given by the last two methods are also shown 
in fig. 2 labelled (53-54) and A respectively. Although there is substantial 
agreement between the estimates given by the three methods, probably 
the most reliable one is obtained from the null method described first. 


§ 3. MreTroroLoaicaL Data 


We shall now examine the evidence for a daily temperature variation 
in the atmosphere. The difficulties of measuring the true air temperature 
during daylight hours have already been indicated but in recent years 
temperature sensing devices with a special coating of very low absorptive 
power for the solar radiation spectrum have been developed (Brasefield 
1948) and adopted for general use in routine upper air soundings in the 
U.S.A. An analysis of the upper air data obtained during the summer 
months of 1956 and 1957 at Washington, D.C. has been made by Harris 
(1959) who gives the diurnal and semi-diurnal components of the tempera- 
ture variation in the atmosphere at altitude intervals of 1km between 
0-5km and 9-5 km inclusive. 

The atmospheric contribution to the solar daily variation of cosmic 
ray meson intensity was calculated by Dorman (1957) from the upper 
air data compiled by Kaminer and it is therefore interesting to compare 
these data from the U.S.S.R. with the more recent results of Harris. 
The data for the same season as Harris’s results were analysed and differ- 
ences between the two sets of results were found to exist in both amplitude 
and phase. The differences can best be seen by resolving each 24hr wave 
of temperature into two components of the same periodicity, of which 
one has its time of maximum fixed at mid-day and the other one at 18 hr. 
Only the amplitudes of these components need be compared as the phases 
are constant. The advantage of this resolution is that the amplitude of 
only one of the components, viz. the one with a maximum at noon, will 
be sensitive to radiation errors in the measurements whilst the other will 
not be affected to any appreciable extent. This is because the radiation 
errors at times placed symmetrically with respect to noon should be 
roughly the same. 

We have calculated the amplitudes of these components at various 
isobaric levels from the results of Harris and Kaminer and plotted them 
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in fig. 3 as functions of the atmospheric pressure. It can be seen from 
the diagram that the two sets of results give roughly the same amplitudes 
at various pressure levels for the component with an 18hr maximum. 
Marked discrepancies are evident, however, in the other component, with 
the results of Harris showing a roughly constant amplitude of ~0-1°c 
between 750mb and 350mb while the Russian data give systematically 
larger amplitudes which increase with altitude from 0-5°c to about 1-2°c 
over the same range of the atmosphere. 


Fig. 3 


(a) COMPONENT WAVE WITH A MAXIMUM (b) COMPONENT WAVE WITH A MAXIMUM 
AT 2 Hr. AT 18 Hr. 


— Harris’s Data 


~Harrls’s Data 


— Kaminer's Data 


— Kaminer's Data 


AMPLITUDE OF TEMPERATURE VARIATION — °C 
AMPLITUDE OF TEMPERATURE VARIATION — °C 


ro) 200 400 600 800 [000 oO 200 400 600 600 [000 
ATMOSPHERIC PRESSURE-mb. ATMOSPHERIC PRESSURE — mb. 


The diurnal temperature variation at various isobaric surfaces in the atmos- 
phere according to the data of Harris (solid curve) and Kaminer (dotted 
curve). The amplitudes of the resolved components with a diurnal 
maximum at (a) 12 noon and (b) 18 hr L.T. are shown. 


In determining the reasons for this divergence between the two results, 
factors such as the wide differences in the epochs of measurement and the 
geographical locations of the observing stations have to be taken into 
account. But there can be little doubt that at least a part of it is due to 
the differences in the relative amounts of radiation errors in the temperature 
measurements. The rough agreement between the two sets of data 
regarding the component which is not affected by radiation errors and the 
growth of the difference in the other component with increasing altitude 
appear to support this view. 
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We shall next examine the effects of the differences in aerological data 
on the estimates of the atmospheric contribution to the meson intensity. 
The temperature data of Harris together with the densities of temperature 
coefficients calculated by Dorman have been used to compute the percen- 
tage changes in the counting rate of a cubical meson telescope at sea-level 
originating in the portion of the atmosphere between sea-level and 250 mb, 
the limits being set by the availability of temperature data. The calcula- 
tions give a first harmonic with an amplitude of 0-12°% and a time of 
maximum at 05-30hr, and a second harmonic with corresponding values 
of 0:04% and 06-30hr respectively, The daily variation curve of meson 
intensity, which includes the effects of both harmonics, is shown in fig. 4 
and a similar curve obtained by Dorman is also included for comparison 
purposes. The latter curve refers to the same season and the same range 
of the atmosphere, but is mainly diurnal with an amplitude of 0-22% 

‘and a maximum at 02hr. As indicated above, the larger amplitude and 
the advance in phase from 05-30hr to 02hr must partly be due to the 
larger radiation errors in Kaminer’s results. 


Fig. 4 


Harris's data 


Kaminer’s data 


pe 
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N 
a) 
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The solar daily variation due to atmospheric temperature in the counting rate 
of a cubical meson telescope at sea-level, calculated from the temperature 
data of Harris (solid curve). An earlier estimate of Dorman (dotted 
curve) is shown for comparison. 


The contribution from that portion of the atmosphere above 250mb 
cannot be estimated with any comparable accuracy but Kay (1951) has 
suggested that the temperature variation in the 150-80 mb region has an 
amplitude of <0-5°c, reaching a maximum value in the late afternoon, 
and if this is accepted as correct, our diurnal amplitude of 0-12% should 
be increased by about 0-03 to 0:05%. 
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As a result of the rapid attenuation of the diurnal amplitude with 
height evident in Harris’s data, about two-thirds of the 0-12°% amplitude 
originates in the atmosphere below the 750mb isobar. It is well known, 
however, that the daily temperature variations in the lower troposphere 
show a marked seasonal dependence and vary considerably with geo- 
graphical location. It can be expected, therefore, that the meson intensity 
variations will exhibit correspondingly large fluctuations from season to 
season at temperate latitudes and that the annual mean of the variation 
will have a smaller amplitude than the summer value of about 0-16% in 
agreement with the suggestion of Baliga and Thambyahpillai (1959). 
Moreover, the magnitudes of the temperature effects may vary substan- 
tially from one cosmic ray station to another. 

We conclude, however, that the temperature wave deduced in § 2 from 
cosmic ray data is not inconsistent with that given by aerological measure- 
ments. Of course, Huancayo is at an altitude of 3350m above sea-level 
and because of the lower mean energy of the »-mesons recorded there, 
the densities of temperature coefficients should be somewhat larger than 
at sea-level, but this cannot lead to more than a fractional change in the 
amplitude of the expected variation. This should, however, be offset by 
the integration over a smaller range of the atmosphere (0mb to ~ 650 mb) 
and the smaller amplitudes of the diurnal temperature variation in this 
range. These factors might explain the somewhat smaller amplitude of 
0-11°% deduced for the mountain station from cosmic ray measurements 
as compared with an amplitude of about 0-16°% at sea-level. Moreover, 
the time of maximum of 05-40 hr obtained in § 2 appears to be borne out 
by the temperature data of Harris. The remarkable agreement between 
the two independent derivations with regard to the phase leaves little 
room for doubt that the diurnal maximum of the meson intensity variations 
due to atmospheric temperature occurs around 05-30hr and not in the 
early hours of the morning, as has hitherto been believed (Dorman 1957). 


§ 4. DrrEcTIONAL MEASUREMENTS MADE AT Mawson 


In this section we shall consider the effects of correcting solar daily 
variation results for atmospheric influences using the temperature vector 
derived above. For this purpose we have selected the measurements 
made during the period May 1955 to February 1956 at Mawson, Antarctica 
(67° 36'S, 62°53’ E geographic) with a vertical telescope and another 
inclined at 45° to the vertical and rotated at regular intervals so as to 
look, in turn, due geographic north, east, south and west (Parsons 1957). 
The data, as published, were corrected only for changes of barometric 
pressure. These results were believed to be inconsistent with the view 
that the daily variation is the result of an anisotropic distribution of the 
primary flux, existing outside the earth’s magnetic field (Parsons) but 
this conclusion should be re-examined after taking into account the 
contribution from atmospheric temperature changes. Moreover, the 
dipole approximation to the earth’s magnetic field used by Parsons in his 


Atmospheric Temperature Effects 595 


discussion may not be sufficiently accurate for the purpose. Quenby and 
Webber (1959) have already pointed out the need for including the 
non-dipole part of the field in the calculation of cosmic ray cut-off rigidities 
and these terms are also of importance in calculating the asymptotic 
directions of all low energy cosmic rays. 

The diurnal components of the solar daily variations measured by 
Parsons after the application of pressure corrections are shown in table 2 
which gives the amplitudes and times of maxima for the various directions 
in which the telescopes were pointed. The features of interest in the 
results are the small amplitude of 0-014 + 0-016% for the south direction as 
compared with amplitudes lying between 0-06 and 0-08% for the other 
directions and the appreciable spread in the times of maxima. 


Table 2. The solar diurnal variations of cosmic-ray intensity at Mawson 
measured by telescopes pointing (a) in the vertical direction and (6) 
at zenith angles of 45° in various directions. The results corrected 
for barometer changes are shown together with the values obtained 
when a temperature corrections (vide text) is also included 


eee Pressure corrected eee ea iy 
eles 
direction ; : 
Time of max. Time of max. 
Ampl. (%) (hr (L.T.)) Ampl. (%) (L.T.) 
North 0:081+0-016 14-50+ 00-45 | 0-:179+0-025 16-28 + 00-32 
Kast 0:060+0-016 13-20+01-00 | 0-:147+0-025 16-124 00-40 
South 0:014+0-016 13-50+04-20 | 0-118+0-025 17-12 + 00-45 
West 0-:064+0-016 16-10+00-55 | 0-172+0-025 17-08 + 00-36 
Vertical 0-079 +0-006 13-45+00-15 | 0-:166+0-020 16-04 + 00-30 


In the absence of reliable information about the temperature changes over 
Mawson and the differences in temperature effects for vertical and inclined 
telescopes, we have used the temperature vector derived for the Huancayo 
ion chamber (amplitude 0-11°%: phase 05-40 hr) to correct all five of the 
24hr waves obtained by Mawson. This procedure can give only a rough | 
approximation to the actual state of affairs; but the uncertainty in the 
amplitude of the temperature correction has no marked effect on the 
conclusions to be drawn from the results. The diurnal variations obtained 
after the temperature correction are also given in table 2 and are shown on 
the harmonic dial in fig. 5. As a result of the correction, the ratio of the 
amplitudes in the south and north directions has been increased to about 
two-thirds and the spread in the phases for the various directions has been 
reduced. 

In discussing his results in terms of a primary anisotropy, Parsons has 
adopted the suggestion of Brunberg and Dattner (1954) that the deflection 
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of the primary particles in the earth’s magnetic field can be adequately 
taken into account by assuming a suitable mean value, between 20GV 
and 40GV, for the magnetic rigidity of the primary particles which are 
responsible for the daily variations. It is then possible to determine the 
asymptotic direction from which particles of this rigidity entered the 
earth’s field before being deflected so as to be recorded by the appropriate 
counter telescope. This asymptotic direction is characterized by a latitude 
angle, measured from the geographic equatorial plane, and a longitude 
angle measured eastwards from the geographic meridian plane through 
Mawson. Using the data of Brunberg and Dattner (1953), based on the 
dipole approximation to the geomagnetic field, Parsons calculated these 
angles defining the asymptotic directions for the five telescope directions 


0.10%, 


12 hr. L.T. 


Harmonic dial representation of the cosmic ray diurnal variations at Mawson 
in the vertical direction (V) and at zenith angles of 45° in the north (N), 
east (E), south (S) and west (W) directions after correction for the 


effects of atmospheric pressure and temperature. Standard deviation 
circles are also shown. 


used at Mawson and for particle rigidities of 15, 20 and 40GV. The 
calculated angles of latitude and longitude are shown diagrammatically 
in fig. 6. Parsons has pointed out the following difficulties in reconciling 
his results (uncorrected for temperature effects) with an anisotropic 
distribution of the primaries which exists outside the earth’s magnetic 
field. (a) Although the asymptotic latitude for the south pointing telescope 
was intermediate between those for the other telescope directions, the 
south telescope measured a much smaller amplitude than the other tele: 
scopes. This result cannot be explained easily because the strength of 
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the anisotropy can be expected to show a smooth variation with latitude. 
It will be seen from table 2 that this difficulty is largely removed by the 
application of a temperature correction. (6) The longitude angles are 


Fig. 6 


DIPOLE FIELD 


(¢) 
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40GV 
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The angles of latitude and longitude defining the asymptotic directions from 


which primary particles enter the earth’s magnetic field before being 
deflected so as to be recorded by the various telescopes at Mawson. 
Values for primary rigidities of 15, 20 and 40 GV are given. 


(a) and (6), asymptotic geographic latitudes; 

(c) and (d), asymptotic geographic longitudes (measured as in text); 

(a) and (c), take only the dipole part of the earth’s magnetic field 
into account; 

(b) and (d), take into account the dipole and non-dipole parts of the 
earth’s magnetic field. 
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different for the five telescopes and an excess of radiation entering the 
magnetic field from a particular direction should give rise to maxima in 
the measured intensity at different times. The sequence of the times of 
maxima should be in the order of decreasing longitude angle but the 
arrangement of measured times of maxima disagreed with the expected 
sequence. 

Now, the geomagnetic latitude of Mawson on the centred dipole approxi- 
mation is 73°, though the dip latitude, as defined by the angle of dip of the 
actual field lines, is 54°. Besides, the declinations of the horizontal com- 
ponents of the actual field and the field due to the centred dipole differ 
by 26°. Thus, the non-dipole part of the earth’s magnetic field is quite 
large at Mawson and has to be taken into account. To do this a new dipole 
was fitted by requiring that the field of the fitted dipole gave a best fit to 
the actual field strength along the radius vector through Mawson. The 
angles of latitude and longitude corrected for the non-dipole field effects 
are also shown in fig. 6., It will be seen that there is no significant change 
in the distribution of latitude angles for 15 GV and 20 GV particles. Since 
the measured amplitudes, after temperature correction, are approximately 
equal within the errors, there is still no major inconsistency with the idea 
of a primary anisotropy. 

The effect of including the non-dipole part of the earth’s magnetic field 
is to bring the asymptotic longitudes for 15GV and 20GV closer to one 
another and this gives better agreement with the measured diurnal varia- 
tions because the phases, after temperature correction, lie within a range 
of 20° while the standard deviations are about + 10° on each measurement. 

In order to determine the direction and the mean rigidity of the primary 
anisotropy, the angle 6 of the direction of maximum primary intensity 
measured eastwards from the sun-earth line was calculated from the 
equation : 

6=%,4 15°(tmax— 12) 

where ys, is the appropriate longitude angle corrected for non-dipole field 
effects and tmax is the time of maximum in hours of the diurnal variation 
when all atmospheric influences have been eliminated. 6 is plotted in 
fig. 7 as a function of the assumed mean rigidity of the primaries for the 
five telescope directions at Mawson. If the measured diurnal variations 
were due to a primary anisotropy, all curves should intersect at a point 
corresponding to the mean rigidity and the direction of the anisotropy. 
It will be seen from the diagram that all curves tend to converge at a 
primary rigidity of 10GV and the curves separate out at higher rigidities. 
Thus, within the limits of statistical accuracy of the measurements, the 
results appear to be consistent with a primary anisotropy located roughly 
70° due east of the sun—earth line, when the mean rigidity is taken to be 
about 10-20GV._ It is concluded, therefore, that the Mawson results can 
be reconciled with the idea of a primary anisotropy lying outside the 
earth’s magnetic field, when the atmospheric influences and the deflecting 
action of the non-dipole part of the geomagnetic field are taken into 
account. 
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§ 5. CONCLUSIONS 


Reasons have been given for believing that the corrections for atmos- 
pheric temperature effects on daily variation measurements of cosmic-ray 
intensity have not been evaluated accurately to date. Using purely 
cosmic-ray measurements, the atmospheric part of the solar diurnal varia- 
tion has been found for the ionization chamber data from Huancayo. 


LINE 


EAST OF SUN-EARTH 


Q@ DEGREES 


RIGIDITY (GV.) 


The longitude angle 6 between the sun—earth line and the direction of maximum 
primary intensity outside the geomagnetic field, as given by each 
telescope at Mawson, plotted against the assumed values for the mean 
rigidity of the primaries producing the diurnal variation. The standard 
deviations of 9, which are independent of rigidity, are also shown. 


This has an amplitude of 0-11°% and a time of maximum around 05-40 hr 
L.T. Also, using some recent meteorological data, considered to be largely | 
free from insolation errors, the atmospheric contribution has been calculated 
on the Dorman scheme for a cubical meson telescope operated at sea-level 
in temperate latitudes. The diurnal component here has an amplitude of 


600 On Atmospheric Temperature Effects 


about 0-16% and a time of maximum at 05-30 hr during the summer season. 
Whilst the small difference in the amplitudes given by these two independent 
derivations is to be expected, the remarkable agreement in phase shows 
that the diurnal maximum of temperature in the troposphere is reached 
near 18hr and not in the early afternoon as indicated by earlier radio- 
sonde measurements. When the latter calculation is compared with 
an estimate of Dorman for an ionization chamber at sea-level, it appears 
that the increased amplitude of 0-27% and the advance of phase to 02 hr 
found by Dorman are probably due to larger insolation errors in the 
aerological data used. 

As an application of the calculated temperature effects, the directional 
telescope measurements made at Mawson were corrected using the Huan- 
cayo temperature vector. These corrected results are consistent with the 
view that the daily variation arises from a primary anisotropy outside the 
earth’s magnetic field when the deflecting action of the geomagnetic field 
is taken into account. The results suggest that the primary radiation 
producing the daily variations possesses a mean rigidity in the range 
10-20GV and that the maximum intensity arrives from a direction 70° 
due east of the sun—earth line. A more comprehensive treatment of various 
aspects of the solar daily variation, taking into account the atmospheric 
effects discussed here, is now in preparation. 
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ABSTRACT 


It is shown that the general characteristics of the cosmic- -ray intensity 
variations can be explained in terms of a large-scale interplanetary magnetic 
field of predominantly dipole character but containing small-scale irregularities 
which act as scattering centres. It is suggested that the cosmic-ray data 
can be taken as evidence for the existence of such a field. The strength of the 
field, which is generated by current systems in the solar corona, is dependent 
on the level of solar activity but must in general be in the region of 10-% 
to 10-4 gauss at the earth’s orbit. 


§ 1. INTRODUCTION 


THE primary cosmic rays show certain characteristic variations in intensity 
which are rather closely correlated with solar and geomagnetic activity and 
it is generally believed that this relationship exists because of the modula- 
tion of the cosmic-ray flux by changing magnetic, and possibly electric, 
fields in interplanetary space. The central problem in the study of cosmic- 
ray intensity variations is that of constructing a detailed picture of these 
fields from the experimental data available and it is the purpose of this paper 
to discuss one particular model which appears to fit these data rather well. 
For any model to be acceptable, the configuration of fields chosen must be 
capable of explaining the eleven-year variation in cosmic-ray intensity 
over the sunspot cycle, the Forbush type decrease usually associated with 
great magnetic storms, the quasi-periodic 27-day variation, and the solar 
daily variation. Finally, it must be compatible with the propagation 
characteristics from sun to earth of particles accelerated to cosmic-ray 
energies in solar flares. 

During recent years several models of the interplanetary electric and 
magnetic fields have been proposed but none has so far proved completely 
satisfactory. Broadly speaking, these proposals may be divided into two 
categories; those which depend on randomly orientated fields of a kind 
generated by the turbulent motion of highly ionized gas (Morrison 1955, 
Meyer et al. 1956, Singer 1957, Parker 1958) and those which depend onrather 
highly organized and uniform magnetic fields associated with solar streams 
(e.g. Alfvén 1954, Brunberg and Dattner 1954, Cocconi e¢ al. 1957, Dorman 
1957). Since these magnetic fields vary with time it follows that there must 
also be electric fields associated with them and in some models these electric 
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fields play an important, even dominant, role. The possibility of an 
extensive interplanetary electric field capable of modulating the cosmic 
ray intensity has recently been discussed by Ehmert (1959). 

The modulating mechanism which will be discussed here has some 
features in common with these earlier models but differs from them in 
that it depends primarily on the existence of a general solar magnetic field 
which extends far out into interplanetary space beyond the earth’s orbit 
and which is essentially non-turbulent. The screening effect of a solar 
dipole field on the cosmic-ray intensity at the earth was first discussed by 
Janossy (1937) many years ago when it was believed that the sun possessed a 
dipole moment of ~ 1034gausscm?, corresponding to a polar field of about 
30 gauss. We now know that the solar dipole field, as measured in the 
photosphere, is much smaller than this ( ~ 1 gauss at the pole) and that it is 
variable (Babcock 1959). Such a field would not in itself be capable of 
producing any appreciable effect on the cosmic-ray intensity at the earth. 
If, however, large-scale current systems exist in the solar corona the effec- 
tive dipole moment of the sun may be much greater than would appear from 
measurements of Zeeman splitting of spectral lines originating in the 
photosphere. In the present model it is suggested that such currents do 
exist and that the field they produce is responsible for the cosmic-ray 
intensity variations. If we are to make a quantitative estimate of the 
effects of such a field on the cosmic rays, it is necessary to idealize and a 
close relation between the model and reality must not be expected at this 
stage. 

§ 2. PRoposED MopE&L oF THE SOLAR MAGNETIC FIELD 


In attempting to build up a picture of the interplanetary field from the 
cosmic-ray evidence we must begin with the simplest configuration capable 
of accounting for the qualitative behaviour of the cosmic rays and then 
elaborate it to the extent merited by the observational detail available. 
There are two simple models which lend themselves readily to quantitative . 
treatment. In the first of these it is supposed that the field changes 
rapidly in intensity and direction as a function of space coordinates. This 
corresponds physically to a highly turbulent field in which the motion of 
the cosmic-ray particles can be treated as a diffusion problem without the 
necessity of any knowledge about details of their trajectories. This 
model has already been extensively discussed in the references given 
above. 

At the other extreme, we can suppose that the current systems generating 
the field lie predominantly within the earth’s orbit and possess a high degree 
of axial symmetry about the sun generating a nearly dipole field at points 
outside the earth’s orbit. The motion of cosmic-ray particles in a dipole 
field is well understood and a quantitative treatment of this model is 
also possible. Close to the sun the field must be very disturbed and 
far from regular because of patchy photospheric fields, sunspot fields and 
the violent hydromagnetic motions of the solar atmosphere close to the 
surface. The occurrence of magnetic storms can be taken as evidence that 
these violent disturbances sometimes extend out as far as the earth’s orbit 
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and presumably beyond it. For the purpose of calculation we shall divide 
the solar field schematically into three regions as shown in fig. 1, an inner 
region in which the field is much disturbed and which is surrounded by 
another region in which the large scale axially symmetrical currents flow. 
These currents generate a dipole field which constitutes the third and outer 
region extending far beyond the earth’s orbit. This dipole field will be 
periodically disturbed by the outward moving solar streams responsible 
for magnetic disturbances and auroral activity. 


Fig. 1 


Schematic representation of the magnetic field generated by currents flowing 
in the solar corona. 


§ 3. Errect OF THE SOLAR FIELD ON THE Cosmic-Ray INTENSITY 


If the field outside the earth’s orbit was a completely unperturbed dipole 
field we should have the familiar allowed and forbidden regions given by the 
Stormer theory. That is to say, we should have some minimum magnetic 
rigidity, p,,;, below which no particle could reach the earth and a maximum 
Pax #bove which particles could arrive from all directions. In the inter- 
vening range betewen p,,,,, and p,,,,, the arrival of particles would be restric- 
ted to certain directions. The values of p,,,. and p,,, are given by 


P90 s ana 3p. =(4/2-1 20 


Pina ™ Tre ra InGGVs sxe (Ly 


where 7 is the radius of the earth’s orbit in centimetres (assumed circular), 
Mg is the dipole moment of the solar current system in gauss cm’, Z is the 
particle charge in units of the electron charge e. For a singly charged 
particle Z=1 and the magnetic rigidity in Gv is numerically equal to the 
particle momentum in Gev/c. 

2R2 
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In reality the situation is more complicated than this because of the 
perturbations introduced by outward-moving solar streams. If, as seems 
likely, these streams can produce gradients in the field of a sufficient 
magnitude, particles will be scattered into forbidden regions which would 
otherwise be empty according to the Stormer theory. The requirement 
for such scattering to be effective is that the field should change by a sub- 
stantial fraction in a distance corresponding to the radius of curvature 
of the particle trajectory in the unperturbed field. Alternatively the 
perturbing field, being time dependent, should change appreciably during 
the time the particle requires to cross the perturbing region. 


Fig. 2 


Scattering of cosmic-ray particle by irregularity in dipole field. 


§ 4. ScaTTERING BY Maqnetic IRREGULARITIES 


Near the equatorial plane of the solar dipole, where scattering by solar 
‘streams is most likely to occur, a cosmic-ray particle describes a nearly 
helical path with its axis perpendicular to the ecliptic plane. If such a 
particle encounters a sharp gradient in the field the radius of curvature of 
its trajectory at that point is changed and the particle is scattered into 
another orbit. As a result, the guiding centre is scattered through some 
distance A, The projection of the trajectory in the equatorial plane is 
shown in fig. 2. 

A detailed discussion of the scattering process is only possible if we have 
detailed information about the structure of the scattering region. In the 
absence of such information we can only make a rough schematic treatment 
as follows. 
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We suppose that the field strength in the irregularity is of the same order 
as that outside and that the physical dimensions of the irregularity are 
small compared with p, the radius of curvature in the unperturbed field. 
This means that on passing through the irregularity the particle suffers a 
small angular deflection 
_, 300 Bil 

Pp 
where 5; is the mean field strength along the particle trajectory in the 
perturbing field and / is its length. 

Then A,~ 6p ~(B;/B)l where B is the local value of the unperturbed 
field. If we suppose that the perturbing field B, is linearly related to B, 
so that B;/6=constant, it follows that A; is independent of particle rigidity 
and of position. We can then regard the motion of the cosmic-ray particle 
in the solar dipole field as a combination of that described by Stormer 
theory and diffussion with a constant transport mean free path given by A, 
above. 

In the absence of absorption, the effect of the diffusion would be to fill all 
the forbidden regions with particles which, by Liouville’s theorem, would 
eventually attain the same intensity as that existing in interstellar space. 
In fact, as we shall see, absorption by the sun reduces the intensity in the 
forbidden regions below that existing outside the solar system and an 
equilibrium intensity is reached which is dependent on the relative rates of 
diffusion and absorption. 


6 


§ 5. ABSORPTION OF Cosmic-RAy PARTICLES BY THE SUN 


In order to estimate the extent to which particles are absorbed by the 
sun during their diffusion inwards we have to take into account the complete 
field distribution and not merely the dipole part. The simplest assumption 
we can make about the inner part of the field is to suppose that it is generated 
by aspherical current shell of radius 7, in which the currents flow along lines 
of latitude. If the current strength is proportional to cosA, where A is the 
latitude, such a current system generates a dipole field outside the shell 
and a uniform field within it as in fig. 3. 

At each scattering collision there is a finite probability that the diffusing _ 
particle can be scattered into such a direction that it will travel directly 
along the field lines into the sun and be absorbed. Thus, in a sense, we may 
regard the solar absorber as spread throughout the solar system. If the 
particle is to reach the sun in this way it must be scattered into a small 
solid angle 2, as shown in fig. 3, such that ©, = 27 (1—cos%) where ¢ is 
given by the ‘ trapping condition’, sin*s=B/B,. That is to say, a particle 
with a pitch angle greater than [(7/2)—7] in the equatorial plane is trans- 
mitted at the surface of the sphere radius r¢ and may enter the turbulent 
region and finally strike the sun. A particle with a pitch angle less than 
[(z/2) —], on the other hand, will be reflected before reaching the surface 


ab T=Te- 
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Thus if the particle is scattered into the solid angle 

20 p= (l —4/ (1B) By) i 
it will enter the turbulent region of the field. This is a necessary, but not 
sufficient, condition for absorption by the sun since after entering the 
turbulent region the particle may either strike the sun or be scattered back 
into the ordered dipole part of the field. If P, is the probability that a 
particle entering the disturbed region will strike the sun, it follows that the 
absorption probability per scattering collision in the outer field is given by 


20.4 

Ps it 
The probability of absorption by the sun will, of course, depend on the 
structure of the field close to the sun and the simplest assumption we can 


P= 


Fig. 3 


Cosmic-ray particles scattered into the solid angle Q4 are no longer trapped by 
the dipole field and may be absorbed by the sun. 


make is that the probability of absorption is proportional to the solid angle 
subtended by the sun at the point where the particle enters the turbulent 
region. ‘Thus, in fig. 4, if the point of entry E is distant r’ from the centre 
of the sun we assume that the probability of absorption is 


P,(r')=4(1—cos8)P, where cosB={1—r,2/r’2}? 
and P¢ is the probability of a large angle scattering, i.e. 


, rs” 
Pal) ~4 Peas: 


In general the full area of the sun will not be accessible because of 
photospheric magnetic fields and because of this the accessible area will be 
to some extent dependent on particle rigidity. However, here we shall 
simply take the area to be a constant, As, so that 
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We must now relate this absorption probability to a particular line of 
force in the external dipole field. The dipole line of force which connects 
to EK cuts the equatorial plane at C at a distance 


2s re 
r= a = ot 
cos?7A or’? 
and therefore we can write 
Agr 
TAU U as = 3° 
c 


Combining this with the expression for Q,, we obtain for the probability 
of absorption per scattering collision 


20 Ag 
P(r) = * P,(r)=Pe— {1— V(1— B/By)} 

A : r.\3 (4) 
i.e. ei Aad ne on putting BiB,=(“) 

ae ; 


The geometry of the inner field region close to the sun. 


§ 6. QUANTITATIVE DISCUSSION OF THE SCREENING BY A SOLAR DipPoLe 
FIELD 

Our problem, then, is to calculate the cosmic ray intensity inside a 
diffusing medium with a transport mean free path A, and a probability of — 
absorption per collision given by P in equation (4). In order to do 80 we 
must first define the outer boundary to the diffusing region. If Mg is the 
magnetic moment of the solar current system and A the latitude angle 
measured from the equatorial plane, then according to Stormer theory, on 
the surface of revolution defined by 


300 MV, cos? A (5) 
SiN aT ccs in ee 
( p [1+ (1 —cos? A)? ] 


particles can arrive from all directions without restriction and on this 
surface cosmic rays of rigidity p or greater will have the full galactic 
intensity. 
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At first sight it might seem that this surface should be taken as the outer 
boundary of the diffusing region. In fact it is much more realistic to take 
instead the surface on which particles cannot arrive directly from any 
direction. Inside this surface particles can only arrive by diffusion whereas, 
as we shall see in §9, the flux outside this surface is very closely isotropic 
and of nearly full galactic intensity. 

This surface, S,, for which direct access is forbidden, is the surface of 
revolution generated by the curve. 


= 300 VM. cos? A 
us p " [1+ (1+ cos? A)t?] ° 


For A between + 30° we can take this to be part of the surface of a sphere 
of radius 7,, where r,, is the value of r corresponding toA=15°. Then 


axl (=S2) es 


and the diffusing region in the dipole field is bounded by the surfaces 
Sp ro loa): 


The boundaries of the region in which diffusion occurs. 


Provided the absorption probability is small the net flux crossing the 
surfaces S, and S,, fig. 5, is very nearly zero and consequently we can regard 
the system as spherically symmetrical so far as H is concerned. 

The problem is therefore simplified to that of calculating the intensity 
distribution inside a diffusing spherical medium of radius 1, = \/(45Ms/p) 
with a transport mean free path A, =B,l/B and a probability of absorption 
per collision P= P,(As/2r?). If N is the number of scattering collisions per 
second and n the number of cosmic rays per unit volume, the appropriate 
diffusion equation is | 


N(A,?/3) V2n(r, p, t)—NP(r) Nr, DP, j= —— < e 4 (7) 


We are interested in the stationary solution and we therefore set = =0 
t 


and put 4, the differential flux of particles of rigidity p, equal to nv. Then 
we have 


(A,?/3) V°9(r,p)—P(r)o(r,p)=0. . . . 2. (8) 


or V24(r,p)—K LEP) =0 where K= 3PcAs 
r De 
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The solution of this equation is 
b= Ar-HMl-VA+4E) 4 Bp + V(+4K)) 
but since ¢ must be finite at r=0 it follows that B=0. Finally using the 


other boundary condition, ¢=¢,, at r= ry, we obtain 
ie #{1—/(1+4K)} 
“(7 
uae Ms) M1 —/(1+4K)} 
=O ar | : ee eee) 


The intensity at the earth is ee given by 


5 M.)7#1-v +45} 
heeds | pea ees 2] 
where r,=1-5 x 10cm and r,2p < 45 Ms. 

In the calculation so far, we have taken A, to be constant but this cannot 
be so at low values of p when the scattering angle, 6 approaches 7/2. Under 
these circumstances, that is to say for large 0, A,>p. We can take this 
into account <a by writing 


(10) 


r=, 


Bil 
h & nl =u = “| 
—_ Winer es and p,=300 Bi. 
This gives A’, constant and equal to B,/B at large values of p and X’,=p at 
p<. Consequently in (10) we can replace K by 


Kt = 3PcAs (Po+P)” _ p[ Pot? |? 
Pas 2 rep? Pp 
giving 
4/ (45 Mg) |e - VG + 4 BU (D0 +9) 1)} 


bn=Pa (11) 


re pie 


§ 7. THe 11-yEAR VARIATION IN INTENSITY 


Equation (11) gives the steady-state cosmic-ray flux at the earth in terms 
of the galactic flux ¢,,, as a function of magnetic rigidity p. It seems likely 
that WM, will vary over the solar cycle and so will the probability of capture 
by the sun, P, because of its dependence on ©, (eqn. (3)). In deducing - 
the value of 2, used above we have assumed a strictly dipole field and 
therefore a particular relationship between B and By, (eqn. (2)). At the 
time of a major solar disturbance, however, this field must be seriously 
distorted by current systems associated with the outward moving solar 
material, indeed it has been suggested by Cocconi et al. (1957) that at such 
a time there may be a radial magnetic field connecting the earth and sun. 
In this extreme case, or for any increase in field strength in the region of 
solar emission between A= + 30°, the ratio B/B, is increased and so is the 
value of Q,. It follows from (3) that P is also increased and consequently 
at sunspot maximum, when the frequency of large disturbances is greatest, 
the average value of K must be greater than that appropriate to a dipole field 
becoming say wK, where « is a numerical factor greater than 1. According 
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to the present model it is this increase in both K and M, s with increasing 
solar activity which is responsible for the reduction in intensity of low- 
energy cosmic rays. The most convenient way of relating the predictions 
of the model to experimental observations is to compare the predicted 
variation of the ratio 4,,/¢.. with rigidity to that actually observed. 

The differential rigidity spectrum at the earth, $u(P), can be measured 
at any particular time during the solar cycle but since we do not know 
¢.(p), the galactic spectrum, we cannot deduce the radio by(P)/Po(P) 


Fig. 6 


THEORETICAL MODULATION FUNCTION 
M,=3x 10> :«K=0-5 2 p=lS 


05 


fi --~| EXPERIMENTAL MODULATION -]— 
4 FUNCTIONS 
f 


Comparison between calculated rigidity dependence of the 1l-year variation 
and that deduced from observation. 


directly from observation. There are however, two alternative plausible 
assumptions we can make which enable us to make a rough comparison 
between theory and observation. We can suppose that at sunspot 
minimum  ~ 0in which case the galactic spectrum is simply the differential 
spectrum observed at the earth atthattime. Alternatively, we can suppose 
that the galactic spectrum follows the same power law as that observed at 


the earth for high rigidity particles which are not affected by the solar 
modulation process. 
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The function [1—¢,/¢.,] deduced from experimental observations for 
each of these two assumptions (Elliot et al. 1959) is shown in fig. 6 together 
with that calculated for M,=3 x 108, vK =0-5 and p»=1-5av. It will be 
seen that there is reasonable agreement between theory and observation 
for these values of the three parameters. Using these values of «kK and M 
and py we can, with the aid of the differential response functions obtained 
by Webber and Quenby (1959), deduce the ratio of the change in intensity 
over the solar cycle for the nucleonic component and for the ,-meson 
intensity recorded at sea level. The calculated ratio of 3 is in rough 
agreement with the observed value of ~ 4 at Ottawa (Fenton et al. 1958). 


Fig. 7 
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‘Comparison between calculated and observed rigidity dependence of the 11-year 
variation for nucleons at mountain altitudes. 


A more detailed comparison between theory and observation is shown in 
fig. 7, where the calculated dependence of the long-term intensity variation 
on station cut-off rigidity is shown together with the observed changes at 
mountain stations for three different epochs of the solar cyle. It will be 
seen that the theoretical curves fit the data quite well provided we assume 
a systematic decrease in Ms, «K and py with decreasing solar activity. In 
this calculation it has been assumed that the galactic spectrum is the same 
as that observed at the earth at sunspot minimum. 


§ 8. ForBUSH-TYPE INTENSITY DECREASES 


At the time of a great magnetic storm it frequently happens that the 
cosmic-ray intensity at sea level drops sharply by from 5 to 10% and then 
slowly recovers, over a period of several days, back to the pre-storm value. 
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According to the present model these Forbush decreases are closely con- 
nected with the 11-year cycle in intensity and are observed when the earth 
lies directly in the path of the outward-streaming material emitted by an 
active region on the sun. If we suppose that this solar material to some 
extent carries the field with it then the trapped cosmic-ray particles 
spiralling round the field lines must be convected outwards at the same time. 
It follows that on such occasions an observer at the earth will sample a sub- 
normal cosmic-ray flux characteristic of that existing under normal 
conditions, at some point lying nearer to the sun. Alternatively, if we 
regard the outward motion of the solar gas as leading to an increased field 
strength in the region through which it passes illustrated in fig. 8, it follows 
that because grad (B) now has a tangential component, the cosmic-ray 
particles trapped in the field will tend to drift round the region of higher 
field, producing the same result, a decrease in intensity. 


Fig. 8 


BY OUTWARD MOVING 
SOLAR MATERIAL 


The outward convection of trapped cosmic radiation at the time of a Forbush 
decrease. 


This will not be the only effect however. An increase in the field strength 
near the ecliptic plane will also lead to a greater probability of particles 
being absorbed by the sun. This arises from the increase in the ratio B /B. 
and therefore in Q,, compared with that for the unperturbed dipole fell 
As pointed out in the previous section, the superposition of many such events 
will be contributory towards the reduction in cosmic-ray intensity at sun- 
spot maximum. Because of the time scales involved, however, this change 
in absorption probability is not very important in discussing the Forbush 
decrease and will be ignored in the following treatment. 
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We can obtain a rough idea of the form of the primary rigidity spectrum 
during a Forbush event in the following way. According to fig. 8, during the 
period of reduced intensity the cosmic rays received at the earth in the 
rigidity range p < (45.Ms/r,2) will be those normally trapped in orbits which 
do not extend out beyond some distance re’ from the sun. This distance 
re will be related to 7 in some rather complicated way, which we shall not 
attempt to deduce here, by the magnitude of the perturbing field. 


Fig. 9 
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Comparison between observed and calculated change in intensity of the 
nucleonic component at sea level during the Forbush decrease of March 
25/26, 1958. The intensity changes have been normalized to the value 
1-0 at high latitudes. 


It is sufficient for the present, to point out that at the time of the Forbush 
decrease we should observe a spectrum of the form 


j 4/ (45 Ms) 3{1—/(1+4aK[(po’ +H)/p}?)} 
$,'(0)=$o(p)| MS” 


where re and p, have been replaced by re’ and py’, and p,’ is greater than py 
because of the greater value of the field at re’. In fact p9' =p) (re/re')*, and 
if we take re/re’ = 1-22, a value which gives a decrease in intensity of 5-8% 
in the nucleonic component at sea level at high latitudes, we find 
Po =3:74av. Using these values for re’ and po’ and taking ae =0-3 and 
M,=10*°, we can calculate ¢,'(p) and the function [1—¢, (p)]- Using 
this, together with the appropriate differential response functions, we can 
calculate the corresponding change in the »-meson component which turns 
out to be 2:65%. The ratio of 2-2 between the changes in the nucleonic 
and meson components agrees well with the mean value of about 2-3 


(12) - 
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obtained by Fenton et al. (1958) for the Resolute and Ottawa stations in 
1957. The change in intensity of the nucleonic component at sea level has 
also been calculated as a function of station magnetic cut-off rigidity using 
eqns. (12). The calculated dependence is shown in fig. 9 together with the 
observed latitude dependence for the Forbush decrease of March 25/26, 
1958. Again the agreement is satisfactory. 


§ 9. THE SoutaR DaiLy VARIATION 


In the discussion so far we have assumed the cosmic-ray flux at the earth 
to be isotropic at all rigidities, but in fact a large-scale solar dipole field of 
the kind discussed above leads to a small anisotropy which would appear as 
a solar daily variation in intensity. This anisotropy arises in the following 
way. 


Fig. 10 


Anisotropy due to the gradient in the cosmic-ray flux. 


In a uniform magnetic field in which there exists a gradient in the particle 
number density the flux at a point # will be anisotropic as shown in fig. 10. 
The fractional magnitude of the anisotropy is given by 5 = 2p (grad d)/¢ and 
an observer rotating with the earth at # and receiving radiation from direc- 
tions close to the ecliptic plane would measure a sinusoidol variation in 
intensity with amplitude 5/2=p (grad ¢)/¢. 

From (12) we have 


: ere > PotPp)? 
a= — sea} af (44a reel 


Taking representative values for the three parameters M,=2 x 10%, 
ak =0-25 and py=1, we can calculate 8/2 as a function of p and this is 
shown in fig. 11 for l1<p<20@v. This calculation is valid in so far as the 
motion of the cosmic ray particles can be treated as a diffusion problem 
but it certainly breaks down for p > 30 av where p becomes comparable with 
the barrier thickness. For values of p> 48 Gv, however, we can obtain a 
rough estimate of the degree of anisotropy in the following way. 
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If the solar field was strictly dipolar in character, particles with rigidity 
less than 
, 300 M. 

ZT 97 
could only reach the earth as a result of scattering, whereas those with 
rigidity greater than 


Pmin.= (/2 1) 


could arrive freely from all directions. In the range between p,,,, and 
Pmax Singly charged particles could only arrive from directions lying outside 
a cone of semi-angle y given by 


peaers 300 MV, _ 300 Ms ath A heati a} 
Y pret pre 


and with its axis in the plane of the ecliptic and perpendicular to the sun— 
earth line (fig. 12). 


Fig. 11 
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The calculated rigidity dependence of the solar daily variation. 


For M,= 2:0 x 10°, pmin= 4°6 x 101° and pmax = 2°7 x 10" so that for this 
value of M, we might expect particles in this rigidity range to be aniso- 
tropically incident on the earth. Maximum contribution to the daily 
intensity variation would be expected from those particles for which y lies 
between about 45° and 135°, that is to say from particles with 
48av<p<1l2av. Taking the direction of the solar dipole to be in the 
same sense as that of the earth and neglecting the deflection in the earth’s 
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magnetic field, maximum intensity would occur at 18 hours local time and 
minimum intensity at 06 hours. Primary particles in the rigidity range 
48 to 112 av produce about 20% of the »-meson flux at sea level and if the 
solar field was entirely dipolar in character we should expect a large daily 
variation with an amplitude ~10%. Because of scattering, however, the 
forbidden directions will to a large extent be filled with trapped particles 
and the magnitude of the anisotropy much reduced (Alfvén (1947), Kane 
et al. (1949), Treiman (1954)). The intensity in the forbidden directions 
will be determined by the probability of absorption by the sun and can be 
estimated very roughly as follows. 


Fig. 12 


The Stormer forbidden cone at the earth due to the dipole field. 


In order that a particle in the momentum range we are considering shall 
be transferred from an allowed intoa forbidden direction it has to be scattered 
through an angle comparable to y given by (13). This angular deflection 
will be made up for a number n, of small deviations = 300B,l/p such that 

ae { 9 eh (= Ms\ _ 300 Ms\ _ 12300 Bil 
Hite? sf Va ADPae ad Dt 
From eqn. (4) and §7 the probability of absorption per collision, P is given 
by P=Pe¢ (Aga/2re”), and taking 
, 8PcAsaf/B \? aye 
aK = oR (5) == (Zo 
as indicated by the 11-year variation, we find 


We AN ETE 
pore 
(3) 12r¢? 


The intensity in the forbidden direction at the earth, ¢,, is therefore given 
by 


2 
a te p ae 300 Ms\ _ 300Ms)}? 
br bo( 12(300Bre)? & PC pre? pre : ue 


Putting in the value of B at the earth’s orbit appropriate to M,=2 x 1035 
and taking p=5~x 10! for instance, we obtain (bo —by)/b = 0-016, 
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indicating that for this rigidity the intensity in the forbidden directions 
is only 1-6% less than that in the allowed directions. At p>67«@v the 
forbidden cone has a half angle, y, less than 7/2 and in consequence can be 
filled to the full galactic intensity because particles can be scattered into it 
with a negligible probability that they will strike the sun in the process. It 
follows that the magnitude of the anisotropy, 5/2, in primary intensity must 
fall from ~0-8% at 50 av to zero at p= 67 av as shown in fig. 11. 

In the intervening region between 20Gv and 48«@v the extent of the 
anisotropy is very uncertain but it seems not unreasonable to join the two 
branches of the curve smoothly together by the dashed line in fig. 11. Any 
attempt to calculate the magnitude of 6 for this part of the spectrum re- 
quires a more detailed specification of the field than has been given so far 
and it seems likely that it is just this dependence on detail which is 
responsible for the great variability of the daily variation. 

If, however, for the present, we accept the curve in fig. 11 as giving an 
approximate estimate of the degree of anisotropy for particle trajectories 
which lie in or near the plane of the ecliptic, we can use it to calculate the 
expected diurnal variation in intensity of the primary particles. In doing 
so we must take account of the deflection of the primary particles by the 
earth’s magnetic field (Brunberg and Dattner 1953)). Finally, using the 
appropriate differential response curves (Webber and Quenby 1959), the 
variation of primary particle intensity can be converted into the corres- 
ponding variation for a particular component of the secondary radiation. 

A calculation of this kind leads to an expected amplitude of ~ 0-43°% for 
the nucleonic component at Huancayo with maximum intensity occurring 
at 13.30 hours local time. These values are in reasonable agreement with 
the experimentally observed amplitude of the nucleonic variation at 
Huancayo of 0-3% and the time of maximum at 12.15 hours (de Moura et 
al. 1959, private communication). This calculation of the daily variation 
to be expected from the solar field is in the nature of a preliminary estimate 
but it serves to show that there is no gross inconsistency between the 
observed daily variation and that expected on the basis of the present 
model. 


§10. THe Propacation TO THE Eartu OF PARTICLES ACCELERATED 
IN SoLAR FLARES 


It is well known that occasionally particles with energies of a few Gev can 
reach the earth from a solar flare and that in the early stages they do so in 
_ well defined impact zones. These impact zones soon disappear, however, 
leaving a slowly decaying flux of solar particles which is isotropic to a 
high degree. Cocconi ef al. (1957) have shown that it is possible to account 
for this sequence of events if the photospheric magnetic fields can be 
extended by outward moving solar material into the form of long tongues 
which trap the solar particles and guide them to the earth. 

In § 8 we have seen that distortion of the general field by outward-moving 
solar material can account for the Forbush decrease. In an extreme case 


P.M. 28 
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the general field may well be dragged out into a tongue with the same capa- 
bilities of guidance and storage as those envisaged in the model proposed by 
Cocconi et al. If we accept their picture we must conclude that the present 
model is in no way incompatible with the propagation characteristics for 
solar flare particles. 


§ 11. SuMMaRY 


It appears from the calculations given above that the general character- 
istics of the cosmic-ray intensity modulation process can be explained in 
terms of a large-scale solar field of predominantly dipole character which 
contains small-scale irregularities acting as scattering centres and which is 
occasionally very seriously distorted by outward-moving solar material. 
In particular this model appears to be capable of giving a semi-quantitative 
account of the 11-year variation, the Forbush-type decrease and the solar 
daily variation. The 27-day intensity variation appears to be basically 
the same phenomenon as the Forbush decrease (Van Heerden and 
Thambyahpillai 1955) and consequently does not require separate 
treatment. 

In fitting the field model to the experimental data on the 11-year cycle 
three adjustable parameters have been used. The first of these Ws, 
determines the strength of the field, assumed to be predominantly dipole, 
whilst the second, «K, is a measure of the probability of absorption by the 
sun of particles diffusing through the dipole field. The values of Ms lie 
in the range 1 to 3 x 10° gauss cm? and those of aK in the range 0-15 to 0-30. 
The third parameter, p, is related to the size and field strength of the scat- 
tering irregularities by the equation py = 300 Bl, and if we take B; ~ (IM g/r¢°) 
and py) ~ 1 Gv, as indicated by the cosmic-ray data, we find /~ 10" cm near 
theearth. No particular care has been taken in the present paper to fit the 
experimental data with great accuracy but it would appear that if we are to 
attempt a more detailed description of the interplanetary field on the basis 
of cosmic-ray data we shall need greater observational accuracy on the 
primary spectrum variations than at present exists. If however, as 
suggested here, there exists a dipole like field with a strength of 3 to 8 x 10-> 
gauss at the earth’s orbit its detection by space probe magnetometer 
measurements cannot be long delayed. 
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Electron-microscopic Observations on Radiation 
Damage in Graphite} 
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ABSTRACT 


Transmission electron microscopy has been applied to the study of lattice 
defects introduced by neutron irradiation of graphite. The dark-field 
technique has proved to be especially useful for this purpose. The technique 
of preparation is described and a tentative interpretation of the observations 
is given. 


§ 1. IyTRODUCTION 


A GREAT deal of research work has been done on radiation damage in 
graphite, and this has been reviewed up to 1956 in various papers presented 
to the “ First Conference on Peaceful Uses of Atomic Energy ’’ (1956). 
An important problem is the mechanism of energy storage. The research 
work presented here, which is a part of a more extended research 
programme, has been undertaken with the aim of studying, by means 
of transmission electron microscopy, lattice defects introduced by neutron 
irradiation. 


§ 2. TECHNIQUES 


For the preparation of the microscope specimens, graphite scrapings, 
suspended in a water—alcohol mixture, were collected with a normal 
200-mesh copper grid, which was covered by a collodion foil full of holes 
of about 1» diameter. Thus, thin graphite flakes could be placed over 
such holes without being backed by amorphous carbon from the collodion 
support. Another technique was similar to that used in taking a replica : 
a Triafol (acetobutyrate) foil wetted, and on its surface dissolved by — 
acetone, was pressed on to a piece of graphite. After the foil had dried 
on the graphite, small particles stuck to it. A small piece of the Triafol 
was placed on the copper-collodion grid in a Soxhlet where the Triafol 
was dissolved in acetone. By this method very thin specimens could be 
obtained. A ‘dry’ method was used by dusting the graphite on to a 
sticky copper grid (the sticky grid was prepared by dipping a normal 
grid into a dilute benzene solution of ‘Scotch-tape * cement). In this 
case the graphite flakes stuck out in all directions and very thin foils 
which are the most interesting ones, were found to be folded. The 
observations were made in a Siemens Elmiskop I at 80 kv. 


+ Communicated by the Author. 


622 W. Bollmann on 


The pictures are taken in dark-field illumination in the ‘light’ of 
1010 reflections selected by the aperture diaphragm, the illumination 
system being tilted to allow the selected beam to pass through the centre 
of the lenses, so that the full resolution of the microscope could be 
obtained. Dark-field, as compared with bright-field, has the advantage 
of much higher contrast because practically only those electrons which 
are diffracted contribute to the image, and the undeviated or diffusely 
scattered electrons are eliminated. Thus the influence on the image of 
the amorphous contamination caused by oil vapours is greatly reduced. 
In addition, only a small group from all the diffraction spots is selected, 
giving detailed information about a small component of the whole 
specimen. From a stack of layer-packets rotated with respect to each 
other, particular ones can be selected by this method. A comparison of 
fig. 4+ and fig. 5 shows the difference in information obtained from the 
same area by dark-field and by bright-field illumination. In the case 
of the bright-field picture, none of the fine details of radiation damage 
can be seen. It is thought that for other applications also, dark-field 
transmission electron microscopy will markedly increase the amount of 
information obtainable, compared with the usual bright-field method. 
It is noteworthy, however, that even if only a part of a stack of graphite 
layers becomes visible by diffraction, the remainder of the material can 
have a harmful influence because of inelastic scattering, which introduces 

a chromatic error, with the effect of fading out fine details. The difference 
in concentration of fine details between the centre and the top of fig. 1 
can be attributed to this effect. Normally, complicated moiré fringe 
systems are observed owing to relative rotation of superposed groups of 
graphite layers. The specimen has to be scanned in dark-field to choose 
areas with especially widely spaced moiré fringes, i.e. especially small 
rotation of superposed layers. The spacing of the fringe in fig. 1 
corresponds to a rotation angle of 16-7’. 


§ 3. RESULTS AND Discussion 


The graphite studied was grade A reactor graphite irradiated by 
107° n/em? at about 30°c. Figures 1 and 2 show the situation as it was 
observed in this material. Black and white dots with diameters ranging 
from 10 to 1004 are visible. In another specimen of the same material 
small black and white loops can also be seen (fig. 3). Unirradiated 
graphite showed none of these dots: in this case only smooth moiré 
fringe systems were visible. Graphite which had been irradiated by 
10% n/em? showed the same type of dots but at a much lower density. 

To study the effects of annealing, specimens were examined of irradiated 
graphite (107° n/em?), heated in vacuum for 3 hours to different tempera- 
tures, i.e. 200°, 400°, 1000°, 2200°c. The strongest annealing effect is 
to be expected around 200°c as the differential energy release there is 


{ All figures are shown as plates. 
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a maximum (Woods ef al. 1956). Figure 4 shows a specimen heated 
to 400°c. Both a certain coarsening of the loops and a decrease of the 
dimensions of the white and black dots are observed. In addition 
certain lines appear, which might be dislocations. Black and white 
dots were also observed in a specimen annealed at 1000°c but a specimen 
heated to 2200°c showed the same appearance as unirradiated material. 
The parallel annealing experiments of unirradiated graphite showed 
none of these features. 

To obtain complete understanding of the nature of these dots, the 
contrast due to electron diffraction effects will have to be studied in 
detail. Tilting experiments have shown that the appearance of the dots 
can change from white to black and vice versa. The loops are reminiscent 
of those observed by Hirsch e¢ al. (1958) in quenched aluminium, which 
are formed by vacancy clusters. Thus, as a tentative explanation, we 
suggest that the observed defects might be extrinsic and intrinsic stacking 
faults or groups of these corresponding to interstitial and vacancy 
clusters. 

The difficulty in estimating the density of the dots lies in the fact that 
the thickness of the layer is not known. To produce the picture, fig. 1, 
at least three to five monatomic layers of graphite are needed, but probably 
there are more of them. If we suppose about 10 layers, 107° n/cm?, and 
a mean free path of the neutrons of about 3 cm, we observe about 5 to 
10 times fewer dots than expected neutron impacts. Thus, the dots 
could originate from individual neutron impacts, if only the high-energy 
neutrons produce a visible effect. This does not mean that the dot should 
be at the place where the impact actually occurred, because it could mark 
the end of the path of an atom which has been hit, unless further 
rearrangement and annealing happened in the time between irradiation 
and observation. 

Some of the properties of irradiated graphite, e.g. x-ray line broadening 
and differential energy release at 200°C, show a significant change for 
a dose of about 4 to 6x 107n/em? (Woods e al. 1956). At higher 
doses the x-ray lines broaden very strongly and the energy release on 
annealing is more and more displaced to higher temperatures at the 
expense of the release at 200°c. Looking at fig. 1 and considering a - 
6-times higher dose it becomes obvious that a significant part of the 
graphite lattice must be disturbed, or even destroyed, which would be 
compatible with the x-ray line broadening. So, the change from the 
low dose annealing, happening essentially at around 200°c to the high 
dose annealing which is displaced to higher temperature might correspond 
to the change from repairing an existing crystal lattice to nucleating a 
new one. 
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ABSTRACT 


A solution of the non-linear differential equations of a domain boundary 
is reported for the case of a thin film in which the demagnetizing energy 
associated with the magnetization at right angles to the plane of the film is 
included. The solution leads to a proposal for the internal structure of the 
cross-tie wall which has been observed in thin films of permalloy. The effects 
of the magnetic poles in this model of a cross-tie wall are discussed and 
compared with observations. The energy density of the wall is suggested 
to be about 10 ergs/em® and the correlation between cross-tie spacing and 
length is derived. 


§ 1. LyTRODUCTION 


Iy order to understand coercivity effects in thin films when the reversal 
process is one of domain wall motion it is important to have detailed 
knowledge of the energies of the domain boundaries. Néel (1955) has 
described how the energy associated with a Bloch wall rises as the thickness 
of a film decreases because of the increasing magnetostatic contribution 
to the wall energy arising from the poles on the surface of the film. When 
this magnetostatic energy becomes sufficiently great it is energetically 
favourable for the magnetization M in the boundary region to rotate in the 
plane of the film. Such a wall is known as a Néel wall and in permalloy 
films of zero magnetocrystalline anisotropy and induced uniaxial aniso- 
tropy of 1500ergs/cm? it occurs below a thickness of about 300A and 
decreases in energy with decreasing film thickness. Further, the Néel wall 
is much thicker than a Bloch wall for the same film thickness because of the 
low value of the induced anisotropy which provides the force tending to 
restrict the wall thickness. 

Huber et al. (1958) have observed a ‘ cross-tie’ structure domain wall in 
thin films of permalloy and have proposed a model for the magnetization 
distribution within the wall. Figure 1 is reproduced from their paper and it 
shows how their model includes both Néel and Bloch like segments. 
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The object of this paper is to describe some calculations made to try to 
derive a cross-tie wall from a perfectly general set of equations for a domain 
wall. Fig. 1 
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Huber, Smith and Goodenough’s model for the cross-tie wall. 


§ 2. SOLUTION OF THE WALL EQuaATIONS 


By describing the behaviours of the amplitude-constant magnetization 
vector M with two angles (6, ¢) defined in fig. 2 and in the presence of a 


Fig. 2 


The definition of coordinate angles. 
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general anisotropy, K(6, 6) which can be chosen to include magnetostatic, 
magnetoelastic and magnetocrystalline terms Harrington (1958) has shown 
that the conditions to be satisfied by a static boundary in the absence of an 
applied field and in a material with exchange constant A are 


2A V?0=2A sin @cos 6(V¢)?+0K/00, . . . . . . (1) 
2A sin? 6V*6 = —4A sin 8 cos 6(VO. Vd) + 0K /d¢. i ae ea 


Harrington derived these equations in a semi-classical manner by 
formulating a Langrangian of the system and then using a variational 
treatment to obtain a set of Euler-Lagrange equations. It was from these 
EKuler—Lagrange equations that eqns. (1) and (2) were deduced. 

In a thin permalloy film the total anisotropy K can be represented by 
two terms: 

(i) A term K sin? 6 expressing the fact that such films show an induced 
uniaxial anisotropy and 

(ii) A term proportional to the square of the component of magnetization 
at right angles to the plane of the film which occurs because of the magneto- 
static energy associated with magnetic poles at the surface of the film. 
Thus the total anisotropy term can be written as 


K(6,46)=Ksin?6(1+bcos?f) . . . .. . (3) 


where 0 is a constant determined by the magnetostatic energy density of a 
thin film magnetized uniformly perpendicular to its plane. 

The solution of eqns. (1) and (2) for a given set of boundary conditions is 
not straightforward and in order to proceed from the K (4, 4) in (3) above 
let us choose to find a solution in which 6 varies normal to the wall and ¢ 
varies along the wall so that the derivatives of the two angles are orthogonal 


and V@. V¢=0. 
Thus the two equations become: 
Ad?6/dy?= A sin 6 cos 6(06/0x)* + Ksin 6 cos 6(1 +6 cos? 4), (4) 
PA Oe Ole =. — COB SIOD SE ee she Pe a ee Sei ety ge (0) 
Equation (5) may be directly solved giving 
L—Xy= (ez) log tan (f/2+7/4) . . . . . (6) 


inserting (@¢/@x)? in (4) and solving yields the result : 


A 
y= I (Rarer) ooo s Siassdalic Mage L) 


Equations (6) and (7) together show that the wall thickness (y direction) 
varies with the value of ¢ and hence with the value of x A plot of y at 
§=175° against 2 is shown in fig. 3 using the following values 

A=10-%ergs/em 

K =1500ergs/cm3 (i.e. H;, ~ 4 oersteds), 

b =(47M2/K) = 8-5 x 10°. 
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This value of 6 has been used because the energy associated with the 
demagnetizing field normal to the plane of the film is 4 M?. 


Fig. 3 
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The variation of wall thickness for @=5° and @=175°. 


It is clear that the solution represented by eqns. (6) and (7) is not cyclic in 
x so that it cannot be the solution for a cross-tie wall. However, the 
following points are to be noticed. 

(i) This solution contains a Bloch segment and two Néel segments and 
so is a compromise between pure Bloch and Néel walls in the same way as 
the original interpretation by Huber, Smith and Goodenough of the cross- 
tie wall. 

(ii) Using the constants above, the distance between the centre of the 
Bloch segment and the Néel segment is about 300 A which is a distance very 
short compared with the diameter of a typical film. The magnetostatic 
energy associated with the free poles in the Néel segments has not been 
taken into account in the calculation above and it is likely that the Néel 
segments will be reduced in length in order to minimize this energy. 
Further, there will be a magnetostatic interaction between the poles in 
adjacent Néel segments. This will be discussed later. 

(iii) The wall thickness in the Néel region is of the order of 15 0004 
which compares reasonably with the cross-tie width as shown in the paper 
by Huber e¢ al. and in a paper by Moon (1959). 
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(iv) With increasing film thickness the value of b falls and both the length 
and the width of the Bloch segment increase so tending to a simple Bloch 
wall. The wall thickness in the Néel segment decreases with increasing 
uniaxial anisotropy K as is expected. 

It is suggested that the solution above could in fact represent one cycle 
of a cross-tie wall and the reason that a cyclic solution has not been found is 
associated with the fact that the original formulation of the domain wall 
equations did not take into account the energetic consequences of a non- 
constant value of div M in the wall. These magnetostatic interactions will 
be dependent upon the film thickness and will be important in determining 
the width and spacing of the Néel segments. In fig. 4 is drawn a model for 


Fig. 4 
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Proposed model for the magnetization distribution within a cross-tie wall. 


the proposed structure of a cross-tie wall based upon the solution derived 
above. The principal difference between this model and that of Huber et 
al. is that the wide sections of the wall are Néel segments and do not contain 
a magnetization vector at right angles to the plane of the film at the centre 
of the wall. This model is still consistent with the shift observed in a 
powder pattern when a field is applied in a plane normal to the surface of the 
film. 

The magnetic colloid will reveal a cross-tie type of wall structure for this 
magnetization distribution because it will lie in the places where the mag- 
netic field at the surface of the film is high. Thus it will collect along the 
length of the wall because here the magnetization intersects the film 
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surface and there are magnetic poles, and along the cross-ties because here 
there is a stray field linking the poles in the two halves of Néel segment 
on each side of the wall. 


§3. Tue Errects or MaGNetTostatTic INTERACTIONS 


Before attempting to evaluate the effects of magnetostatic interactions in 
the cross-tie wall it is convenient to consider the thickness and energy of the 
simple Néel type of boundary. Ifthe energy density and width of a wall in 
the bulk material are written as o, and a, respectively then the energy 
density in a Néel wall of width a in a film of thickness D can be written 
(Néet 1955) as: 


e a(< % “) by 2n7aDM? 

2\a) @ a+D 
The equilibrium wall energy and width for a film of given thickness can be 
found from this equation by putting dH/da=0. The variation of # and a 
with film thickness is shown in fig. 5. The variation in the width of the 
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The thickness and energy density of Néel wall in permalloy film using 
o9=0-28 ergs/em?, ag=2x10-*cm, 4=10-* ergs/em and K =103 ergs/cm’. 


wall with changing film thickness is attributed to the fact that as the film 
thickness increases the area of free poles in the wall increases and as these 
are of opposite sign in each half of the wall they attract each other 
increasingly and the wall becomes thinner. 

Next it is instructive to consider the wall shown in fig. 4 when the free 
poles in the Néel segments are regarded as being at a point distance p from 
the axis of the wall and the cross ties are separated by a distance gq. Itis 
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easily shown by considering nearest neighbours and magnetostatic effects 
alone that the wall is stable when 


256p* = (q? + 4p?) 
1-0: q~1-5p. 

The dimension p is determined primarily by the free poles in the Néel 
segments alone, consequently as the film thickness decreases the dimensions 
p and q both increase. This is consistent with the observations of 
Rubinstein and Spain (1959) who report that the spacing between cross ties 


increases as the films thickness decreases and that the larger the spacing 
between cross ties the longer is the cross tie. 


§ 4. THE WaAxLL ENERGY 


The sum of the exchange energy and the anisotropy energy at a point 
(8, ) in the wall is 


f=2K sin? (1+ 2b cos? 4). 


Neglecting the magnetostatic contributions the wall energy per unit length 
can then be written 


a= [-_ fdy=[" flaylasyas. 
-o / 0 
Thus 
o,=4,/(AK(1+ 26 cos? )]. 
An approximate term for the average wall energy density o,, can then be 


derived by assuming that the solution above is one for a cyclic wall and 
averaging the energy o, over one cycle. 


Thus 


77/2, 
[7 ovtaelagyag 
es a A STC 
Ww 
Cv/(A/bK) 
where C is a constant which allows for the fact that the cycle length is in 
fact longer than \/(A/bK). It can be seen from fig. 3 (where the value of 
,/(e2A/bK) is marked) that a reasonable value for Cis about 2-5. Equation 


(8) can then be expressed as 


vo Mt EB) 


where a=1/2b. The integral is a complete elliptic integral of the second 
kind which is equal to unity because b> 1, therefore 


DMN) 2A). sii Ses tees etre tg. (9) 


For the values of M and A used above this means that o,, ~ 10 ergs/cm? and 
is independent of film thickness and uniaxial anisotropy K when the 
magnetostatic contributions are neglected. This latter result is reasonable 
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because in this case the energy balance is between the exchange A and the 

demagnetizing field 47M. 
Using Néel’s data for Bloch and Néel walls fig. 6 has been plotted showing 

the variation of wall energy with film thickness for the three types of wall. 


Fig. 6 
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Plots of wall energy density versus film thickness for Bloch, Néel and cross-tie 
walls. 


It is clear that a cross-tie wall of energy density 10 ergs/em* can only be 
expected in thin permalloy films in the thickness range 200-1000 A. This is 
in substantial agreement with the work of Methfessel e¢ al (1959) who report 
that they can observe cross-tie walls in thin permalloy films only in the 
thickness range 400-900 A. 


§ 5. CONCLUSIONS 


Although it has not been possible to derive a cross-tie wall directly from 
the general wall equations a solution describing a wall half-way between a 
Bloch and a Néel boundary has been found. This solution appears to be 
similar to one-cycle of a cross-tie wall and makes possible the presentation 
of a more plausible model for such a wall than that suggested by Huber, 
Smith and Goodenough. The energy balance in this wall is between the 
exchange energy and the demagnetizing energy contained in the component 
of the magnetization at right angles to the plane of the film. The spacing 
and lengths of the cross-ties are determined by the magnetostatic interac- 
tions between the cross-ties arising because they are Néel type regions and 
so contain magnetic poles. An approximate calculation of the energy of 
this wall reveals that permalloy films should only exhibit cross-tie walls 
when between 200 and 10004 thick. Thisis in fair agreement with observa- 
tion. 
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The Generation of Vacancies in Metals} 


By R. 8. Barnzs 
Atomic Energy Research Establishment, Harwell, Berkshire 
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ABSTRACT 


The regions where injected helium atoms first appear as gas bubbles are 
those where thermal vacancies originate. Microscopic examination of metals, 
injected with helium by using them as targets for energetic alpha-particles, 
reveals grain boundaries as the principal suppliers of vacancies. Detailed 
examination distinguishes those boundaries which generate vacancies from 
those which merely conduct them. The behaviour and nature of these 
boundaries and other vacancy sources are discussed. 


§ 1. [InTRODUCTION 


SINCE the concept of the vacant lattice site (or vacancy) was introduced by 
Frenkel, it has become increasingly clear that this point defect plays an 
important role in many physical processes of crystalline solids; particu- 
larly at high temperatures where the equilibrium vacancy concentration 
rises steeply, but even at low temperature when non-equilibrium vacancy 
concentrations have been introduced by quenching, deformation, or by 
atomic particle bombardment. A recently developed technique (Barnes 
et al. 1958) enables the sources of vacancies to be observed. These observa- 
tions and some of the deductions from them will be given in this paper. 
Mono-energetic alpha-particles from a cyclotron have been used to bom- 
bard a metal surface and inject a layer of helium atoms beneath the surface 
at a depth equal to the range of the alpha-particles in the metal (see fig. 1). 
The helium atoms in this layer, which is normally quite thin (~ 10~* cm), 
are unable to migrate once their initial momentum is spent unless the tem- 
perature of the metal is sufficiently high. Ifthe temperature is below this 
value the helium atoms remain dispersed in ‘solution’ in the metal, the 
concentration steadily increasing far beyond the solubility limit as the 
alpha-particles continue to be fired in. High concentrations of helium can: 
be retained in solution in this way. On heating such a bombarded metal 
the helium atoms become mobile at some temperature and precipitate from 
solution with the ultimate formation of gas bubbles. Gas bubbles have been 
osberved on heating copper (Barnes et al. 1958), aluminium (Barnes 1959), 
beryllium (Barnes and Redding 1959) and uranium which has been injected 
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with helium. Some study has also been made of gas inserted by other means 
into the metals (Ells and Perryman 1959, Mykura 1959, Rich et al. 1959). 
Figure 2+ is a photomicrograph of a sectioned copper single-crystal injected 
with helium by bombardment with 38 ev alpha-particles and then heated 
to 700°c for lhour. The bombarded surface is seen at the top of the figure, 


Fig. 1 
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The injection of helium atoms by bombardment with alpha-particles. 


and the layer containing the helium bubbles can be clearly seen at the 
bottom{. It is observed that in the first instant only those helium atoms 
accessible to vacancies precipitate ; only the periphery of the layer contain- 
ing bubbles. This is the region immediately accessible to vacancies from 
outside the layer. After shorter times (or lower temperatures) the bubbles 
are even more marginal, and as the heating progresses the vacancies pene- 
trate the layer further until the whole of the layer is filled with bubbles. 
These observations show that there are no vacancy sources within the 
helium layer in single crystals of copper. By heating for a time insuffic- 
ient to allow the complete precipitation of the helium it is possible to outline 
with bubbles, the sources of the vacancies necessary for the process, and 
so observe these sources. 


§ 2. Vacancy SOURCES 


2.1. Free Surface 


That the free surface of the metal is a source of vacancies is demonstrated 
by fig. 3. This shows that although the heat treatment has been insufficient 
to fill the helium layer with bubbles, near the free surface (lying on the left 
of the figure), the layer is filled. (The free surface has been covered with an 


i Figures 2-11 are shown as plates. 
{ In this sample the individual bubbles cannot be resolved, yet: the layer can 


be She mainly due to its lower reflectivity than the rest of the electro-polished 
surface. 
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electro-deposit of copper after the heat treatment to preserve flatness of the 
micrographed surface during polishing.) Vacancies have apparently 
formed at the free surface during heating, and have migrated inwards a 
distance demarked by the limit of the bubbles. The shape of this demar- 
cation line must be an index of the concentration of the helium in the layer. 
If a given number of vacancies per helium atom is necessary before bubbles 
can form, then the penetration of the vacancies is inversely proportional to 
the helium concentration. If this is so the shape suggests that the helium 
is assymmetrically distributed about the mean range of the alpha-particles, 
and this is borne out by the distribution of the bubbles once they have 
formed. 

Evidence that the free surface is a source of vacancies is provided also 
at each heat treatment, for it is only after a sufficiently thick layer has been 
removed from the surface that an unbubbled region is seen. For instance, 
fig. 4 shows the same area (a) after a thin surface layer has been removed 
and (b) after much more has been removed, when the vacancies from the 
original surface are not contributing to the helium bubbles. 


2.2. Grain Boundaries 


Vacancies are also generated at grain boundaries within the helium 
layer, Figure 4 shows this in the same piece of copper as that seen in fig. 3, 
but after heating at 825°c for 1 hour, and the penetration of vacancies into 
the layer is very similar, suggesting that the grain boundary supplies 
vacancies at the same rate as the free surface. Both figs. 3 and 4 were of 
copper which had been well annealed (heated to 1050°c for 48 hours in 
vacuo) before the helium was injected. Other samples of copper with a 
lower temperature of anneal had much smaller grain size. This material 
showed similar behaviour at the grain boundaries (see fig. 5), but due to the 
more numerous sources the vacanices were more plentiful. Photographs of 
samples sectioned in a plane through the centre of the band and parallel to 
the bombarded surface (parallel sections) are shown in figs. 5and 6. These 
show clearly that most of the grain boundaries are supplying vacancies 
to the helium in the surrounding metal. 


2.3. Twin Boundaries 


There are many coherent annealing twin boundaries shown in figs. 4, 5 
and 6. These are never surrounded by bubbles and thus are never 
vacancy sources. As such a coherent twin boundary is not surrounded by 
bubbles, even where it intersects a grain boundary which is, suggests that 
itis also unable to conduct vacancies. However, a non-coherent twin bound- 
ary, although normally free of bubbles, frequently is surrounded by them 
if it connects with a grain boundary, suggesting that such a boundary is 
conducting vacancies from the true vacancy source. Not all non-coherent 
boundaries do this, e.g. fig. 6 shows two non-coherent boundaries (NC) of 
similar nature, both intersecting bubbled regions but only one of these is 
bubbled itself. 
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There are certain twin boundaries which invariably contribute vacancies 
to their surroundings. These are boundaries separating two individual 
annealing twins. In every such case the boundary contributes vacancies 
equally as well as a normal grain boundary. Two examples of boundaries 
(one being very short) between the two twins T, and T, of the parent grain P 
are shown in figure 6. 

2.4. Point Sources 

In addition to the vacancy sources already described there are point 
sources which give rise to small circular regions of bubbles (see, for example, 
fig. 5). These point sources can in a few cases be identified as being due to 
voids within the metal. Figure 7 shows two photomicrographs of a void in 
aluminium. The first (a) is of a section almost through the centre of the 
void and (0) is of a subsequent section cutting the void below its centre. 
Such a void would be able to contribute vacancies to its surroundings until 
it was completely filled in by the aluminium atoms deposited there. The 
bubble free zone immediately surrounding the void is presumably due to 
these deposited atoms. 

Inclusions which at the heating temperature are surrounded by a gap due 
to differential expansion should behave as a void, but due to the smallness 
of the free volume available the number of vacancies generated there should 
be limited. Figure 8 illustrates this behaviour at an inclusion (1) in copper. 
The inclusion is surrounded by a very narrow circle of extremely small 
bubbles (~ 300A in diameter) which can only just be discerned in the 
figure but have been revealed by electron micrographs of replicas from this 
region. In comparison the voidt, which lies very close to the inclusion in 
fig. 8, has supplied sufficient vacancies for bubbles to have formed a large 
distance from it. 

Other point sources, which do not contain a helium free zone (see the 
electron micrograph of a replica in fig. 9) are less well understood. As can 
be seen, there are very few in material which has been well annealed (figs. 2, 
3 and 4) but they are very numerous in material heated sufficient only for 
recrystallisation before injecting the helium (fig. 5). A piece of copper 
(that shown in fig. 6) which had been well annealed was cold-wrought and 
then recrystallized before injecting it with helium. The material was then 
found to contain point sources, but these were fewer in number and more 
irregularly sited than those seen in the rolled material treated identically, 
but not annealed at a high temperature (cf. fig. 6 with 5). It is unlikely 
that these point sources are due to voids, as it is difficult to believe that 
voids could be reintroduced merely by working the material. 


2.5. Dislocation Lines 
On no occasion in any of the material studied have dislocation lines them - 
selves been identified as supplying vacancies to the helium atoms. The 
Se Sa Sener eT 


+ The void has permitted a greater penetration of those alpha-particles which 
have passed through it, as is shown by bubbles lying beneath the void which 


have prevented locally the advance of the grain boundary running along the 
base of the gas layer. 


Generation of Vacancies in Metals 639 


well-annealed material contains ~10* dislocation lines em-2 and the 
recrystallized material many more, but in neither case are bubbled regions 
observed on a scale as fine as this. Samples have been purposely bent or 
cold-worked in other ways, but although slip occurred in the region con- 
taining the helium atoms, the subsequent precipitation of the helium 
followed the normal pattern for the undeformed material, despite the 
inevitable introduction of dislocations into the helium layer. 


§ 3. Deposition or ATOMS 


One consequence of the generation of vacancies is that there must be a 
deposition of atoms at the source. This deposition arises from the manner 
in which vacancies migrate and should be observable. A vacancy migrates 
by exchanging places with a neighbouring atom, so that although each 
atom only moves one jump, there is in effect always a counter flow of atoms 
in equal number to the flow of vacancies. Thus the generation of a vacancy 
at a source and its subsequent migration to a helium bubble is equivalent to 
the migration of an atom from the bubble to the source, where it is deposited. 
Thus for each vacancy generated at a source an atom of the matrix is 
deposited there. 

The helium bubbles, which are inert markers in the crystal, will move 
away relative to the source as these atoms are deposited. Thus at at each 
source there should be a region free of bubbles containing as many atoms of 
the matrix as there are vacancies generated. This phenomena can be 
observed for instance in fig. 7 (a) and (b) which shows a zone free of bubbles, 
at the surface of the void. The width of the zone free of bubbles is rela- 
tively much wider in the lower section (fig. 7 6) as would be expected. Event- 
ually the void can be completely filled and will then cease to be a source of 
vacancies as occurs at the inclusion shown in fig. 8. The majority of point 
sources do not contain such a region free of bubbles suggesting that they 
operate in a different way. 

Figure 10(a) is a shadowed two-stage carbon replica of the electro- 
polished surface of a region in the helium filled layer in copper after heating 
for two hours at 800°c. The individual bubbles are easily resolved and 
each casts a shadow. The diameter of the bubbles can be estimated from 
the width of these shadows, but these may over-estimate due to enlarge- 
ment by the etchant. Figure 10 (b) shows a thin film of copper (~ 30004 
thick) viewed in transmission in the electron microscope. The bubble size 
is here truely represented as the majority of them are completely embedded 
in the copper (diameter 3 x 10-> cm). Both these photographs show grain 
boundaries running across them and at each there is a region free from 
bubbles, similar to that observed at the void in fig. 7. 

If the explanation for this free zone is correct then the volume of the 
bubblesformed per unit area of grain boundaryshould be equal to the volume 
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+ No such clear zone can be seen at the free surface shown in fig. 3. This may 
have been etched away during the preparation for the electro-deposition, or 
smeared out by surface diffusion. 
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of the material deposited upon that area of the boundary. Rough estimates 
from fig. 10, where the average bubble radius is taken to be 1-5 x 10-°cm 
and the number of bubbles ~ 10!2cm~* show that the volume of the deposit 
is about 10-4cm3 and the equivalent bubble volume 1-3 x 10-+cm*; agree- 
ment well within the accuracy. The number of vacancies generated at 
each atom site of the grain boundary is thus 3 x 10° during the two hours 
heating at 800°c and the deposited atoms produce a local volume increase 
around the boundary which in the above example is about 11%. 

Not all grain boundaries show this effect, suggesting that only some of the 
grain boundaries (those with helium-free zones) are really acting as sources 
the rest are merely conducting vacancies into the crystal from elsewhere. 
Some local diffusion of the helium to those boundaries surrounded by bubbles 
occurs as some small bubbles can always be seen on the boundaries them - 
selves. The helium free region cannot be caused by diffusion of the helium 
however, as this would occur at all boundaries equally. 

In general the grain boundaries which are the most copious vacancy 
sources are those which lie in, or nearly in, the plane of the helium layer. 
This behaviourcan be understood when itis realised that if a grain boundary 
is to be a true source of vacancies it must have atoms deposited upon it, and 
the region must expand laterally. A grain boundary lying perpendicular to 
the helium layer will be heavily restrained by the metal above and below the 
layer and such expansion will be difficult. Grain boundaries lying in the 
plane of the layer will have no such restraint. 

Boundaries which do not have this helium free region frequently contri- 
bute fewer vacancies into the surrounding crystal as the distance from the 
real source of vacancies increases (e.g. some boundaries shown in fig. 5). 
This adds to the evidence that they are only conducting vacancies from this 
source. 


§ 4. Discussion 


The observations show that the first bubbles to form surround both the 
true vacancy sources and also those features which are merely conducting 
vacancies from these sources. A source can be distinguished from a vacancy 
conductor by the region of precipitated metal atoms, devoid of helium 
bubbles, adjacent to it. The vacancy sources and conductors will be dis- 
cussed separately. 


4.1. Vacancy Sources 


The free surface is the most straightforward source of vacancies. Ona 
clean surface there will be no restraint to reduce the normal atomic fluctua- 
tions which occur. A surface film, e.g. an oxide layer, may retard such 
a vacancy source, but if it is not a coherent one the effect should be slight. 

The situation at a grain boundary is similar to that at a free surface with 
the added complication caused by the superimposed crystal of a different 
orientation. The inherent structure of a boundary is expected to determine 
its ability as a vacancy source. If there is a space in the boundary of 
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approximately an atomic size then there would be little hinderence to the 
generation of a vacancy. Boundaries having little free space, ie. nearly 
coherent boundaries, would not be good sources. However, not only must 
there be spaces in the boundary large enough for extra atoms to be forced 
in but spaces must then form elsewhere in the boundary for the process to 
continue. 

A hydrostatic pressure (P) upon a sample generating vacancies would 
necessitate an extra amount of work Pa’ (where a? is the volume of an atom) 
to be performed in generating a vacancy. A pressure of more than 104 
atmospheres must be applied before this extra work (~7x 10-ev) is 
comparable with that needed to create a vacancy at a free surface (~ lev). 
If such a pressure was applied then the number of vacancies emanating 
from each and every source would be limited. Such high pressures cannot 
be maintained within these samples, so that local pressure, e.g. those 
caused by unequal vacancy production at different boundaries, cannot 
influence the emanation of vacancies from a boundary. 

However, internal pressures, although unable to prevent the formation 
of bubbles, can change the behaviour of a boundary from that of a source 
to that of a conductor by forcing the atoms which are normally deposited 
on the boundary to flow along it to some other boundary (not subjected to 
this pressure) where they are deposited. A laternal pressure of about 40 
atmospheres (600 p.s.i.) is estimated to be sufficient to cause this flow of 
atoms out of a boundary in a typical experiment. Such pressures can be 
attained and maintained, and explain why boundaries lying in the plane 
of the helium layer contain the zone free of bubbles near to them and those 
lying perpendicularly do not. Such pressures around a cylindrical or 
spherical source would also be expected to force the atoms to escape by any 
possible path. If no such path existed then the deposited atoms would 
build up a pressure sufficient to cause the surrounding material to yield. 

Thus external pressure, or pressures imposed by the material itself, can 
induce a vacancy source to act as a conductor. But this cannot explain 
why a coherent twin boundary is never a source of vacancies. This must 
be due to the nature of the boundary itself. 


4.2. Vacancy Conductors 


Certain crystal defects are unable to generate vacancies but can conduct 
them7+. Such conductors are distinguished from true sources by the absence 
of an intervening region of deposited metal atoms between them and the 
helium bubbles. There is no change in volume in such a region as the 
vacancies are conducted in from their true source where the actual volume 
increase occurs. ‘The conduction process is thus unaffected by restraining 
pressures except in so far as the source itself is affected. 


+ Strictly, the concept of a vacancy in a boundary is meaningless because of 
the absence of a simple crystallographic arrangement there. But the end effect 
of the diffusion of atoms (not necessarily a particular atom) along the boundary 
is equivalent to the diffusion of a vacancy along the boundary. 
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A grain boundary which appears to be a good source is sometimes con- 
verted into a poor conductor where one of the adjacent grains is twinned 
(fig. 11) and occasionally a change in boundary direction prevents a 
boundary from even conducting vacancies. These effects cannot be 
caused by local pressure changes and clearly show the importance of the 
nature of a boundary upon its behaviour. 

As a conductor can only conduct a limited number of vacancies from the 
true source, and as the vacancies are ‘pumped’ away into the helium 
bubbles in the surrounding crystal, the supply becomes exhausted at a 
distance from the real source which is dependent upon the relative 
diffusion rates of vacancies along the conductor and through the surrounding 
crystal. In small grain material the penetration of vacancies from the 
source boundaries into the surrounding crystal is the same for all, but in 
large grain material the distance is very variable. This is because in the 
latter case the conduction paths are much longer so that exhaustion of the 
supply occurs more frequently. 

Occasionally a conducting boundary is seen in a parallel section, where 
the bubbled region tapers to a point, the vacancy supply being exhausted 
there (see fig. 5). In some of these cases the penetration of vacancies from 
the adjacent vacancy source is more than the width of the unbubbled 
region of the helium layer, i.e. more than the penetration in a direction 
perpendicular to the photograph. Thus in these particular boundaries 
the conduction of vacancies is anisotropic. Vacancies conducted along 
the boundary from outside the layer are unable to travel as far as those 
conducted in the plane of layer. In other cases the exhaustion occurs 
in less than the width of the unbubbled region and here the boundary is 
probably a poor conductor in all directions. 

The circular regions of bubbles observed in figs. 5 and 6 suggest that 
vacancies are emanating from a point. They cannot arise from voids 
which were originally present but are now completely filled as there is no 
central region devoid of helium bubbles (see fig. 9). The complete absence 
of such central regions also demonstrated that these spheres of bubbles 
are not caused by point sources of vacancies. Vacancies must be supplied 
to the region and then radiated from the central point. As these regions 
do not in general lie on any observable boundary the vacancies must be 
conducted in via some unseen defect, presumably dislocation lines. 
Certainly the number of circular regions is very much smaller and their 
size more non-uniform in well annealed material (in fact none can be seen 
in figs. 2, 3 and 4 where the number of dislocation lines is ~10%em-?), In 
this material the number of circular regions is only ~104cm-2 
(r~2x 10cm), i.e. the number of spherical regions is ~ 2x 107cem-3, 
whereas the number of dislocation nodes must be ~ 10° em-%, almost two 
orders of magnitude greater. Thus if the dislocation nodes are the 
originating points, only a few of them are operating. The regions might 
originate upon some special dislocation configuration which can climb and 
thus supply vacancies itself. 
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4.3. Dislocation Lines 


It is surprising that the majority of the dislocation lines which are known 
to exist in the metal are unable to climb and thus supply vacancies to the 
helium atoms and permit precipitation as bubbles. Dislocations of both 
edge and screw character are known to be sinks for supersaturated point 
defects and it seems reasonable that they should behave reversibly when 
the vacancy concentration is depleted. The explanation may be that the 
dislocations are effectively locked by the helium atoms lying upon them 
and are thus unable to climb. Certainly the material in the layer con- 
taining helium is much harder (the microhardness is approximately doubled) 
presumably because the helium atoms impede slip of dislocation lines. 


§ 5. STRUCTURE OF BOUNDARIES 


A coherent twin boundary, which is never either a source or a conductor 
of vacancies, is a particularly simple boundary, lying ona(111) plane. The 
diagram (fig. 12) shows the atom positions of a (111) plane as dots. The 
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Diagram of the atomic arrangement at a non-coherent boundary NC in a face- 
centred-cubic crystal. 


plane of atoms above this is represented. as crosses. The right hand side 
can be seen to have suffered a twinning shear, i.e. the passage of a partial 
dislocation with Burgers vector b on each plane. If the successive (111) 
planes in a perfect crystal are labelled ABCABC then the fault at the twin 
plane can be represented by the sequence ABACBA. This is a perfectly 
fitting boundary of very low energy and it sould be surprising if it was able 
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to generate or conduct vacancies, for the formation energy would be almost 
as high as than in perfect crystal (i.e. ~ 5ev). 

A non-coherent twin boundary is that which separates the twinned from 
the untwinned material and yet does not contain the twinning vector b. 
Such a boundary (NC) is represented in the diagram (fig. 12) and must 
contain an appreciably lower density of atoms than the perfect crystal. 
Thus a boundary of this nature might be expected to conduct vacancies 
particularly if it lay normal tob. As every third (111) plane (the B plane) 
is common to both the parent and twin, the non-coherent twin boundary 
will consist of well-defined channels, lying parallel to the (111) twin plane. 
It should be relatively easy for atoms (and thus effectively vacancies) to 
flow along these channels, but most difficult for them to flow in any other 
direction. The boundaries would be expected to conduct vacancies very 
anisotropically, as appears to be the case. Due to the tie-rod action of 
every third layer, it would be difficult for extra atoms to be accommodated 
in this boundary, and it would thus not be a source of vacancies, in 
the sense we have defined. 

The number of possible configurations of a grain boundary is enormous. 
Some simple boundaries may be considered. A simple tilt boundary 
between two grains with a small angle of misorientation may be treated as 
a collection of edge dislocations and a simple twist boundary as a collection 
of screw dislocations. Such boundaries have low energy and appear to be 
unable to create vacancies, but may conduct them in the same way as it is 
thought that dislocation lines may. The conduction of such boundaries 
should be highly directional. High-energy boundaries have considerable 
atomic misfit along them and such boundaries should generally be good 
sources of vacancies. Itis a general observation that low energy boundaries 
are rarely sources of vacancies but high energy ones frequently are. 


§ 6. Tor Errect of DispLACEMENT VACANCIES 


A large number of vacancies is produced during the bombardment of the 
target with alpha-particles, and it should be explained why these are not 
able to permit helium bubbles to form before vacancies from grain boun- 
daries, etc., arrive. One can conclude only that the number of vacancies 
involved is too small. The following simple model illustrates the situation. 

Consider a unit volume of crystal with a simple cubic structure containing 
no vacancy or interstitial atom sinks which is bombarded with alpha- 
particles and has attained an equilibrium situation with NV vacancies and 
an equal number of interstitial atoms. Consider also that the tem- 
perature is such that the interstitial atoms have a high mobility (v atomic 
jumps/sec) and that the vacancies are unable to change their position by 
thermal activation—this is likely to be the situation in copper at room 
temperature. For this equilibrium state to be maintained the number of 
interstitial atoms combining with vacancies in time St must be balanced by 
the number (nét) being produced in that time. In time 8¢ each interstitial 
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atom explores vét interstitial sites and as each of these enables access to at 
least four new vacancy positions the number of vacancies annihilated is 
4(N/N’)vdt where N’ is the number of atomic positions in the unit volume. 
As there are N interstitial atoms, which we shall assume behave inde- 
pendently, 

4(N/N’)vdtN = ndt 


ie. N=4y/(nN'|v). 


Substituting the values for copper at 600°K bombarded in two hours with 
alpha-particles sufficient to produce an atomic concentration of 10-3 helium 
atoms in the layer beneath the surface, gives a value for the equilibrium 
atomic concentration of vacancies to be ~10~°, i.e. only one vacancy per 
10° helium atoms. Thus it is reasonable to assume that the annihilation 
of the vacancies by the interstitial atoms produced at the same time, 
removes the vast majority of the vacancies. The remainder are so few 
that they are unable to assist in the formation of the helium bubbles of the 
size observed. Only if the number of sinks for the interstitial atoms are 
so numerous that more than one hundredth of the interstitial atoms are 
annihilated at positions other than at vacancies would these vacancies be 
expected to affect the observations appreciably. 


§ 7. CONCLUSIONS 


1. The vacancies introduced during the alpha-particle bombardment 
are mainly annihilated by the interstitial atoms also produced, so 
that those remaining are too few to permit the formation of the 
helium bubbles observed. 

2. The necessary vacancies for bubbles to form are drawn from vacancy 
sources in the material by thermal activation. 

3. The nature of a defect boundary in the crystal determines whether 
it is inherently a source, or merely a conductor of vacancies, or 
neither. 

4. A boundary can be restrained from acting as a source, but not as a 
conductor, by the application of a lateral pressure. 

5. Extremely high pressures are needed to prevent the supply of 
vacancies. 

6. A boundary which conducts vacancies can do so at different rates in 
different directions. 

7. Line defects in the crystal are unable to act as vacancy sources in the 
conditions of these experiments, but may conduct vacancies. 
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REVIEWS OF BOOKS 


Solid State Physics in Blectronics and Communications. Vol. 1, Part 1, 
Semiconductors. Edited by M. Dmstrant and J. L. Micurats. 
and New York: Academic Press, 1960.) [Pp. 638.] £6 8s. 6d. 


TuIs is not a book : it is the Report of a Conference held in Brussels in June 
1958. It has all the disadvantages of such a document : (a) eighteen months 
is far too long a delay ; if the papers were worth publishing, they were worth 
publishing quickly, and three more volumes are yet to appear ; (b) many of 
the papers are mere interim accounts of work in progress; they may have 
been worth presenting verbally at the conference, but are not worth printing ; 
(c) there is no reason to suppose that any of these papers were refereed ; they 
do not, therefore, carry the weight of work published in reputable scientific 
journals ; (d) it is more lavishly and expensively printed than such ephemeral 
work deserves ; no private person would afford to buy it for the occasional 
reference ; (e) although libraries may be persuaded into purchasing it, the 
problem of finding papers in it by reference will be much more complicated 
than in the ordinary scientific journals ; even the title is too long for a neat 
footnote, and who can guess how it will be indexed in the book stacks. 

For these reasons I cannot commend the publication of such works, however 
admirable may be some of the actual research that is here reported (and some 
of it is quite interesting). Perhaps the publishers or the organizers of the 
conference can explain the rationale of their procedure, and tell us its advantages 
compared with the publication of ordinary scientific papers through the usual 
channels in the recognized journals. el ae Z: 


(London 


A New Method in the Theory of Superconductivity. By N. N. BocgouruBoy, 
V. V. Totmacuey and D. V. Surrkov. (New York: Consultants Bureau, 
Ine.; London: Chapman & Hall, Ltd.) [Pp.121.] Price $5.75, about £2 10s. 


THE theory of superconductivity has undergone a rapid development in the 
last three years. Indeed it would be fair to say that physicists now have a 
qualitative understanding of many of the main features of the superconducting 
state. The initial impetus for this development came from the now celebrated 
paper of Bardeen, Cooper and Schrieffer. Immediately after this Bogoliubov, 
Tolmachev and Shirkov developed an alternative approach to the theory. 
The volume under review is essentially a translation of this and several other 
Russian papers on the subject. As such it will be of the first importance to 
those theoretical physicists whose Russian is inadequate to deal with the - 
original papers. It is stated on the jacket that the book contains a complete 
solution to the theoretical problems of superconductivity. This, at least in 
the opinion of the reviewer, is not so and it is unlikely that the authors themselves 
would support such a claim. The book does, however, contain a clear account 
of the Russian contributions to the subject. The main criticism that might be 
made of it is that the presentation is rather abstract and that very little physical 
discussion of the theoretical results is given. Nevertheless the serious student 
of the subject will find it particularly valuable for the clear account of the 
quasi particle transformation and the principle of compensation of dangerous 
graphs, for the discussion of the role of the coulomb interaction and for the 
detailed investigation of the collective modes of motion in the superconducting 
state. There are unfortunately several obscurities in the translation, rather 
too many considering the price of the book. G. VC. 
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The Conceptual Foundations of the Statistical Approach in Mechanics. By 
P. and T. Exrenrest. Translated by M. J. Moravesik. (Ithaca, N.Y.: 
Cornell University Press.) [Pp. xiit+114.] Price £1 4s. 

THis is a most welcome translation of one of the few classics in the field of 
statistical mechanics. It has been relatively difficult to get hold of the original 
(at any rate until the appearance last year of the Collected Works of P. 
Ehrenfest) and, as is pointed out by Kac and Uhlenbeck in a short preface, 
the original is couched in rather difficult German. Anybody who is interested 
in the development of statistical mechanics, and especially in its foundations, 
will want to possess this booklet, and luckily its price does not have to be a 
hindrance, as it so often is these days. It is rather unfortunate that the 
translation is far from smooth and reads even less easily than the original. 
A number of German terms were deemed to be untranslatable, but I fail to see 
why recurrence paradox, reversibility paradox, and equivalent would be incor- 
rect renderings of Wiederkehreinwand, Umkehreinwand and gleichberechtigt, to 
name a few instances. A personal element is introduced by a special preface 
for this translation (also printed in the Collected Works) by Mrs. Ehrenfest 
consisting of a few remarks which were omitted in the original on her instigation, 
but which she feels should now be included. D. TER Haar. 


BOOK NOTICES 


Confluent Hypergeometric Functions. By L. J. Suarer. (Cambridge: University 
Press, 1960.) [Pp. 247.] Price £3 5s. 

A Simple Approach to Electronic Computers. By E. H.W. HeErsrEr. (London: 
Blackie & Son.) [Pp. 104.] Price 12s. 6d. 

Lectures on Fourier Integrals. By 8. Bocuner. Translated by M. Tenenbaum 
and H. Pollard. Annals of Mathematics Studies, No. 42. (Princeton: 
University Press ; London: Oxford University Press.) Price £2. 
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VAS PHILLIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 57. 


Structure of fresh area of monocrystalline film of silver made by evaporation 
onto heated rock salt. Note simple stacking faults AG’, D’, X’, etc., 
dislocation loops A, and parallel line features B,_,;. Holes appear in 
white contrast. Foil thickness about 6504. x 60 000. 


V. A. PHILLIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 58. 


Region showing line or point defects C, whole dislocations D, stacking faults 
a, 8, stacking-fault tetrahedra E, and triangular loops F which are 
probably Frank sessile dislocations. The simple stacking faults a and 
complex faults 8B show typical extinction bands. Foil thickness 
500-650 A, x LOO. 000. 


VeeAc PHILLIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 59. 


Fig. 4 (a) 


(a) Dark field electron micrograph of fresh area of film showing stacking faults 
in white contrast. Feature I is probably a twin. The interpretation 
of features J, and J, is in doubt (see text). The remaining features are 
stacking faults and may be classified according to shape as described 
in the text. Note groups of parallel faults H, and Hj. x 60000. 


V. A. PHILLIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 60. 


Fig. 4 (6) (continued) 


(6) Light field micrograph of same area as (a) taken immediately following. 
In the time interval between exposures, pairs of simple faults at L and M 
have interacted and become single-bend faults. Faults now appear in 
dark contrast or show finely spaced extinction bands. x 60 000. 


Ve A> PHILLIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 61. 


Typical appearance of film after long thermal stressing in the beam. Contrast 
figs. 2 and 4. The faults here have widened and interacted to produce 
complex faults. Holes, dislocations and other faults act as barriers to 
further dissociation. Holes and thin regions act as pinning points for 
partial and stair-rod dislocations. x 60 000. 


V. A. PHILLIPS Mag. Ser. 8, Vol. 5, Pl. 62. 


Fig. 6 


(2) 

Sequence of micrographs showing effect of thermal stressing on structure. 
Note pile-up of dislocations Q,, Q., R,—-R, against a group of stacking 
faults ; dissociation is least at the head of the pile-up. Q, and Q, 
are on a different but probably parallel {111} plane to that of R,-R,. 
The hole edge W appears to be a source of dislocations, since a second 
fault on a different plane has appeared between exposures (a) and (b) 
forming a single-bend fault, while in (c) a third fault on a {111} plane 
parallel to that occupied by the fault in (a) has appeared. The 
single-bend fault S has widened progressively between (a) and (b) 
and between (b) and (c), being blocked by a hole in (c). Complex 
fault UT has changed to a simple fault T in (c) permitting dislocation 
V in (b) to move to a new position in (c). There appear to be two parallel 
faults at X which have extended between (b) and (c). Dislocations 
at Y and Z in (a) are held up at faults. Z has moved away in (c). 
Note the pinning of partials P and W at holes, in (a). Note the 
tendency of partials associated with faults Q;, Q, and W in (c) to be 
parallel to ¢110) directions in the fault plane. Stair-rod dislocations 
as at the bend in fault S in (a) are parallel to (110) directions in the 
fault plane, x 140 000, 
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Fig. 6 (continued) 
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V.7A. PHILCIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 64. 


Fig. 7 


Sequence of ciné photographs showing combination of two simple faults to 
form a single-bend complex fault. In (a) the two faults SI and 82 are 
visible. In (b) fault S2 has extended but the two partials (retouched 
for reproduction) indicated by arrows are separately visible. In (c) the 
two partials have combined giving a stair-rod dislocation at B (not 
visible) and forming a single-bend fault ABC (angle between faulted 
planes acute). Magnification about x 80 000 parallel to screen marker 
and about x 67 000 perpendicular to this. 


VeCAgerEILLIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 65. 


Sequence of micrographs showing conversion of single-bend fault DEF in (a) 
to a complex double-bend fault DEHGFI in (6). x 140 000. 


| V.-3A. PHIEEIPS Phil. Mag. Ser. 8, Vol. 5, Pl. 66. 


Sequence of micrographs showing changes in fault configuration. Between 
exposures (a) and (b): the double-bend-fault JKLM has become a 
single-bend fault KLM, a simple fault NO has formed between a hole 
at N and another complex fault at O, and a simple fault RS has formed 
between a hole at R and a fault PQ. Between exposures (b) and (c) : 
the single-bend fault KLM has become a simple fault ML, fault PSQ 


has contracted to SQ and its partial appears to be held up at the partials 
of fault SR. x 140 000. 


W. BOLLMANN Phil. Mag. Ser. 8, Vol. 5, Pl. 67. 


Graphite irradiated by 10° n/em? at 30°C, 320 000 x. 


W. BOLLMANN Phil. Mag. Ser. 8, Vol. 5, PI. 68. 


Same type of graphite as fig, showing loops, 320 000 x . 


W. BOLLMANN Phil. Mag. Ser. 8, Vol. 5, Pl. 69. 


Same graphite as fig. 1 annealed at 400°c, 320 000 x . 


W. BOLLMANN Phil. Mag. Ser. 8, Vol. 5, PI. 70. 


Selected area of the specimen fig. 4 in bright-field illumination, 320 000 x. 


R. S. BARNES Phil. Mag. Ser. 8, Vol. 5, Pl. 71. 


The incomplete precipitation of injected helium atoms in a single crystal of 
Johnson Matthey spectroscopically standardized copper during heating 
at 700°c for 1 hour. x 570. 


Fig. 3 


Bubbles at a free surface in copper after heating to 800°c for | hour. 
x 750. 


R. S. BARNES Phil. Mag. Ser. 8, Vol. 5, Pl. 72. 
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The surface exposed during heat treatment of 1 hour at 825°c (a) after a thin 
layer has been removed by mechanical and electro-polishing and (6) after 
the similar removal of a further layer. x 750. 


Fig. 5 


A parallel section of a similar piece of fine grain copper after an anneal of 
1 hour at 750°c. x 570. 


R. S. BARNES Phil. Mag. Ser. 8, Vol. 5, Pl. 73 
Fig. 6. 


x 570. 


Fig. 7 


(2) 


Microphotographs of a void in Johnson Matthey spectroscopically standardized 
aluminium (a) a section almost through the centre and (b) one below 
the centre. x 570. 


R. S. BARNES Phil. Mag. Ser. 8, Vol. 5, Pl. 74. 


Photomicrograph of an inclusion (1) and a void in copper showing their relative 
ability to supply vacancies. x 1000. 
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R. S. BARNES Phil. Mag. Ser. 8, Vol. 5, Pl. 75. 
Fig. 10 


(0) 


Electron micrograph of a grain boundary showing a zone with no bubbles 
adjacent to the boundary (a) shadowed carbon two-stage replica x 7000 
and (6) transmission through a thin film of the metal extracted from the 
helium layer. x 12000. 
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Fig. 11 


A grain boundary which is either a good vacanc i 
y source (central region) or a 
poor conductor, depending upon whether th joini i 
twinned or not. x 3000. aa 0 ae - 


